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Abstract 

We review in a pedagogical way the present status of the impact of non-Gaussianity (NG) on the Cosmic 
Microwave Background (CMB) anisotropics. We first show how to set the initial conditions at second- 
order for the (gauge invariant) CMB anisotropics when some primordial NG is present. However, there 
are many sources of NG in CMB anisotropics, beyond the primordial one, which can contaminate the 
primordial signal. We mainly focus on the NG generated from the post-inflationary evolution of the CMB 
anisotropies at second-order in perturbation theory at large and small angular scales, such as the ones 
generated at the recombination epoch. We show how to derive the equations to study the second-order 
CMB anisotropies and provide analytical computations to evaluate their contamination to primordial 
NG (complemented with numerical examples). We also offer a brief summary of other secondary effects. 
This review requires basic knowledge of the theory of cosmological perturbations at the linear level. 
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1 Introduction 

Cosmic Microwave Background (CMB) anisotropies play a special role in cosmology, as they allow an 
accurate determination of cosmological parameters and may provide a unique probe of the physics of the 
early universe and in particular of the processes that gave origin to the primordial perturbations. 

Cosmological inflation [I] is nowadays considered the dominant paradigm for the generation of the 
initial seeds for structure formation. In the inflationary picture, the primordial cosmological pertur- 
bations are created from quantum fluctuations "redshifted" out of the horizon during an early period 
of accelerated expansion of the universe, where they remain "frozen". They are observable through 
CMB temperature anisotropies (and polarization) and the large-scale clustering properties of the matter 
distribution in the Universe. 

This picture has recently received further spectacular confirmations from the results of the Wilkinson 
Microwave Anisotropy Probe (WMAP) five year set of data [2]. Since the observed cosmological pertur- 
bations are of the order of 10~ 5 , one might think that first-order perturbation theory will be adequate for 
all comparisons with observations. This might not be the case, though. Present [21 [3] and future experi- 
ments dj may be sensitive to the non-linearities of the cosmological perturbations at the level of second- 
or higher-order perturbation theory. The detection of these non-linearities through the non-Gaussianity 
(NG) in the CMB [5] has become one of the primary experimental targets. 

There is one fundamental reason why a positive detection of NG is so relevant: it might help in dis- 
criminating among the various mechanisms for the generation of the cosmological perturbations. Indeed, 
various models of inflation, firmly rooted in modern particle physics theory, predict a significant amount 
of primordial NG generated either during or immediately after inflation when the comoving curvature 
perturbation becomes constant on super-horizon scales [5]- While standard single-field models of slow- 
roll inflation [6] and - in general - two (multi)-field [7] models of inflation predict a tiny level of NG, 
"curvaton"-type models, in which a significant contribution to the curvature perturbation is generated 
after the end of slow-roll inflation by the perturbation in a field which has a negligible effect on inflation, 
may predict a high level of NG [8]. Alternatives to the curvaton model are models where a curvature 
perturbation mode is generated by an inhomogeneity in the decay rate [9l [10] , the mass [11] or the in- 
teraction rate [12] of the particles responsible for the reheating after inflation. Other opportunities for 
generating the curvature perturbations occur at the end of inflation [13], during preheating p3], and at 
a phase-transition producing cosmic strings [15] • Also, within single-field models of inflation, a high level 
of NG can be generated breaking the standard conditions of canonical kinetic terms and initially vaccum 
states: e.g. this is the case of Dirac-Born-Infeld (DBI) models of inflation [16], and initially excited states 
respectively [T7] . For every scenario there exists a well defined prediction for the strength of NG and its 
shape [TH1 US] as a function of the parameters. 

Statistics like the bispectrum and the trispectrum of the CMB can then be used to assess the level 
of primordial NG (and possibly its shape) on various cosmological scales and to discriminate it from 
the one induced by secondary anisotropies and systematic effects [20"1 [2"Tj [22] . A positive detection 
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of a primordial NG in the CMB at some level might therefore confirm and/or rule out a whole class of 
mechanisms by which the cosmological perturbations have been generated. 

Despite the importance of evaluating the impact of primordial NG in a crucial observable like the CMB 
anisotropy, the vast majority of the literature has been devoted to the computation of the bispectrum 
of either the comovig curvature perturbation or the gravitational potential on large scales within given 
inflationary models. These, however, are not the physical quantities which are observed. One should 
instead provide a full prediction for the second-order radiation transfer function. A preliminary step 
towards this goal has been taken in Ref . [23] (see also [24l ESI EH E3 EEJ EH] ) where the full second-order 
radiation transfer function for the CMB anisotropies on large angular scales in a flat universe filled with 
matter and cosmological constant was computed, including the second-order generalization of the Sachs- 
Wolfe effect, both the early and late Integrated Sachs- Wolfe (ISW) effects and the contribution of the 
second-order tensor modesH 

There are many sources of NG in CMB anisotropies, beyond the primordial one. The most relevant 
sources are the so-called secondary anisotropies, which arise after the last scattering epoch. These 
anisotropies can be divided into two categories: scattering secondaries, when the CMB photons scatter 
with electrons along the line of sight, and gravitational secondaries when effects are mediated by gravity 
[31j . Among the scattering secondaries we may list the thermal Sunyaev-Zeldovich effect, where hot 
electrons in clusters transfer energy to the CMB photons, the kinetic Sunyaev-Zeldovich effect produced 
by the bulk motion of the electrons in clusters, the Ostriker-Vishniac effect, produced by bulk motions 
modulated by linear density perturbations, and effects due to reionization processes. The scattering 
secondaries are most significant on small angular scales as density inhomogeneities, bulk and thermal 
motions grow and become sizeable on small length-scales when structure formation proceeds. 

Gravitational secondaries arise from the change in energy of photons when the gravitational potential 
is time-dependent, the ISW effect, and gravitational lensing. At late times, when the Universe becomes 
dominated by the dark energy, the gravitational potential on linear scales starts to decay, causing the 
ISW effect mainly on large angular scales. Other secondaries that result from a time dependent potential 
are the Rees-Sciama effect, produced during the matter-dominated epoch by the time evolution of the 
potential on non- linear scales. 

The fact that the potential never grows appreciably means that most second order effects created by 
gravitational secondaries are generically small compared to those created by scattering ones. However, 
when a photon propagates from the last scattering to us, its path may be deflected because of the gravi- 
tational lensing. This effect does not create anisotropies, but only modifies existing ones. Since photons 
with large wavenumbers k are lensed over many regions (~ k/H, where H is the Hubble rate) along the 
line of sight, the corresponding second-order effect may be sizeable. The three-point function arising from 
the correlation of the gravitational lensing and ISW effects generated by the matter distribution along 
the line of sight [32l [33] and the Sunyaev-Zeldovich effect [M] are large and detectable by Planck [35jl36]. 
A crucial issue is the level of contamination to the extraction of the primordial NG the seconndary effects 
can produce. In Sec. [8] we briefly summarize some recent results about the level of CMB NG generated 
by some of these secondary effects. 

Another relevant source of NG comes from the physics operating at the recombination. A naive esti- 

1 A similar computation of the CMB anisotropies up to third-order from gravitational perturbations has been per- 
formed in Ref. [30] , which is particularly relevant to provide a complete theoretical prcdition for cubic non-linearities 
characterizing the level of NG in the CMB through the connected four-point correlation function (trispectrum) [3U] . 
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mate would tell that these non-linearities are tiny being suppressed by an extra power of the gravitational 
potential. However, the dynamics at recombination is quite involved because all the non-linearities in the 
evolution of the baryon-photon fluid at recombination and the ones coming from general relativity should 
be accounted for. This complicated dynamics might lead to unexpected suppressions or enhancements of 
the NG at recombination. A step towards the evaluation of the three-point correlation function has been 
taken in Ref. [37] where some effects were taken into account in the so-called squeezed triangle limit, 
corresponding to the case when one wavenumber is much smaller than the other two and was outside 
the horizon at recombination. Refs. [38|, 139] (see also [40]). present the computation of the full system 
of Boltzmann equations, describing the evolution of the photon, baryon and Cold Dark Matter (CDM) 
fluids, at second order and neglecting polarization, These equations allow to follow the time evolution of 
the CMB anisotropics at second order on all angular scales from the early epochs, when the cosmological 
perturbations were generated, to the present time, through the recombination era. These calculations 
set the stage for the computation of the full second-order radiation transfer function at all scales and 
for a a generic set of initial conditions specifying the level of primordial non-Gaussianity. Of course, 
for specific effects on small angular scales like Sunyaev-Zel'dovich, gravitational lensing, etc., fully non- 
linear calculations would provide a more accurate estimate of the resulting CMB anisotropy, however, as 
long as the leading contribution to second-order statistics like the bispectrum is concerned, second-order 
perturbation theory suffices. 

The goal of this review is to summarize in a pedagogical form the present status of the evaluation of 
the impact of NG on the CMB anisotropics. This implies first of all determining how to set the initial 
conditions at second-order for the (gauge- invariant) CMB anisotropy when some source of primordial NG 
is present. The second step will be determining how primordial NG flows on small angular scales. In this 
review we will focus ourselves on the study of the second-order effects appearing at the recombination era 
when the CMB anisotropy is left imprinted. We will show how to derive the equations to evaluate CMB 
anisotropies, by computing the Boltzmann equations describing the evolution of the baryon-photon fluid 
up to second order. This permits to follow the time evolution of CMB anisotropies (up to second order) 
on all scales, from the early epoch, when the cosmological perturbations were generated, to the present 
time, through the recombination era. We will also provide the reader with some simplified analytical 
computation to evaluate the contamination of the recombination secondary effects onto the detection 
of primordial NG. The formalism for a more refined numerical analysis is also displayed and results for 
some worked examples will be also reported. The review is mainly based on a series of papers written by 
the authors along the past years on the subject (with various updates) and, as such, follows both a logic 
and a chronological order. It requires knowledge of the theory of cosmological perturbation at the linear 
level (which however we summarize in the Appendices). We have tried to write the different sections in 
a self-contained way. Nevertheless, we alert the reader that the level of complexity increases with the 
number of the sections. 

The review is organized as follows. In Section 2 we provide a simple, but illuminating example to show 
why we do expect some NG present in the CMB anisotropy regardless if there is or not some primordial 
NG. In Section 3 we provide the reader with the necessary tools to study the dynamics at second-order in 
the gravity sector. In Section 4 we show how to set the initial conditions for the primordial NG, while in 
Section 5 we provide a gauge-invariant way to define the CMB temperature anisotropy at second-order on 
large scales. In Section 6 we go to small scales and present the full procedure to compute the Boltzmann 
equations necessary to follow the evolution of the non-linearities from the recombination epoch down 
to the present epoch. Section 7 presents some analytical solutions of the Boltzmann equations in the 



6 



tight coupling limit, along the same lines of what is done at the linear level. The issue of contamination 
is addressed in Section 8, while in Section 9 we offer the reader with an analytical estimate of such 
a contamination. A more refined numerical work is presented in Section 10. Finally, in Section 11 
some conclusions are given. This review has also some hopefully useful appendices: in Appendix A the 
reader can find the energy-momentm tensors at second-order, Appendix B gives the solutions of Eistein 
equations for the perturbed fluids up to second-order, while Appendix C offers the analytical solutions 
of the linearized Boltzmann equations in the tight coupling limit. 

2 Why do we expect NG in the cosmological perturba- 
tions? 

Before tackling the problem of interest - the computation of the cosmological perturbations at second- 
order after the inflationary era- we first provide a simple, but insightful computation, derived from 
Ref. [28], which illustrates why we expect that the cosmological perturbations develop some NG even 
if the latter is not present at some primordial epoch. This example will help the reader to understand 
why the cosmological perturbations are inevitably affected by nonlinearities, beyond those arising at some 
primordial epoch. The reason is clear: gravity is nonlinear and it feeds this property into the cosmological 
perturbations during the post-inflationary evolution of the universe. As gravity talks to all fluids, this 
transmission is inevitable. To be specific, we focus on the CMB anisotropies. We will adopt the Poisson 
gauge which eliminates one scalar degree of freedom from the goi component of the metric and one scalar 
and two vector degrees of freedom from gij . We will use a metric of the form 

ds 2 = -e^dt 2 + 2a(t)oj i dx i dt + a 2 (i)( e - 2 %- + X ij)dx i dx j , (1) 

where a(t) is the scale factor as a function of the cosmic time t, and U{ and Xij are t ne vector and tensor 
peturbation modes respectively. Each metric perturbation can be expanded into a linear (first-order) 
and a second-order part, as for example, the gravitational potential = <J >( - 1 ) + <J>( 2 )/2. However in the 
metric (pQ) the choice of the exponentials greatly helps in computing the relevant expressions, and thus we 
will always keep them where it is convenient. From Eq. ([1]) one recovers at linear order the well-known 
longitudinal gauge while at second order, one finds $^ 2 ^ = (p^ — 2(4>^) 2 and = V^ 2 ^ + 2(ip^) 2 where 
f^ 1 ) and (j)( 2 \ (with cfP-> = and tp^ 1 ' = VE^ 1 ') are the first and second-order gravitational 
potentials in the longitudinal (Poisson) gauge adopted in Refs. [HJ [5] as far as scalar perturbations are 
concerned. 

We now consider the long wavelength modes of the CMB anisotropies, i.e. we focus on scales larger 
than the horizon at last-scattering. We can therefore neglect vector and tensor perturbation modes in the 
metric. For the vector perturbations the reason is that we are they contain gradient terms being produced 
as non-linear combination of scalar-modes and thus they will be more important on small scales (linear 
vector modes are not generated in standard mechanisms for cosmological perturbations, as inflation). 
The tensor contribution can be neglected for two reasons. First, the tensor perturbations produced from 
inflation on large scales give a negligible contribution to the higher-order statistics of the Sachs- Wolfe 
effect being of the order of (powers of) the slow-roll parameters during inflation (this holds for linear 
tensor modes as well as for tensor modes generated by the non-linear evolution of scalar perturbations 
during inflation). 
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Since we are interested in the cosmological perturbations on large scales, that is in perturbations 
whose wavelength is larger than the Hubble radius at last scattering, a local observer would see them in 
the form of a classical - possibly time-dependent - (nearly zero-momentum) homogeneous and isotropic 
background. Therefore, it should be possible to perform a change of coordinates in such a way as to 
absorb the super-Hubble modes and work with a metric of an homogeneous and isotropic Universe (plus, 
of course, cosmological perturbations on scale smaller than the horizon). We split the gravitational 
potential $ as 

$ = + (2) 

where <3?^ stands for the part of the gravitational potential receiving contributions only from the super- 
Hubble modes; & s receives contributions only from the sub-horizon modes 



d?k 
d 3 k 



*' = I -£kW-*B) (3) 



where H is the Hubble rate computed with respect to the cosmic time, H = a/ a, and 6{x) is the step 
function. Analogous definitions hold for the other gravitational potential ^f. 

By construction <E>£ and are a collection of Fourier modes whose wavelengths are larger than the 
horizon length and we may safely neglect their spatial gradients. Therefore &e and ^ are only functions 
of time. This amounts to saying that we can absorb the large-scale perturbations in the metric (JT]) by 
the following redefinitions 

dt = e® e dt, (4) 
a = oe"*'. (5) 

The new metric describes a homogeneous and isotropic Universe 

ds 2 = -dt 2 + a 2 5ij dx i dx j , (6) 

where for simplicity we have not included the sub-horizon modes. On super-horizon scales one can regard 
the Universe as a collection of regions of size of the Hubble radius evolving like unperturbed patches with 
metric © [33]. 

Let us now go back to the quantity we are interested in, namely the anisotropies of the CMB as 
measured today by an observer O. If she/he is interested in the CMB anisotropies at large scales, the 
effect of super-Hubble modes is encoded in the metric During their travel from the last scattering 
surface - to be considered as the emitter point £ - to the observer, the CMB photons suffer a redshift 
determined by the ratio of the emitted frequency us to the observed one wo 

To = T £ ^, (7) 

where To and Tg are the temperatures at the observer point and at the last scattering surface, respec- 
tively. 
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What is then the temperature anisotropy measured by the observer? The expression (JT]) shows that 
the measured large-scale anisotropies are made of two contributions: the intrinsic inhomogeneities in 
the temperature at the last scattering surface and the inhomogeneities in the scaling factor provided by 
the ratio of the frequencies of the photons at the departure and arrival points. Let us first consider the 
second contribution. As the frequency of the photon is the inverse of a time period, we get immediately 
the fully non-linear relation 

wo ojo 

As for the temperature anisotropies coming from the intrinsic temperature fluctuation at the emission 
point, it maybe worth to recall how to obtain this quantity in the longitudinal gauge at first order. By 
expanding the photon energy density p-y oc T*, the intrinsic temperature anisotropies at last scattering are 
given by 5\T £ /T £ = (l/4)5ip 7 //9 7 . One relates the photon energy density fluctuation to the gravitational 
perturbation first by implementing the adiabaticity condition Sip^/py = (4/3)#ip m /p m , where S±p m / p m 
is the relative fluctuation in the matter component, and then using the energy constraint of Einstein 
equations <f>W = -(l/2)5ip m / p m . The result is 5iT £ /T £ = -2<5> l£ /3. Summing this contribution to the 
anisotropies coming from the redshift factor (|8j) expanded at first order provides the standard (linear) 
Sachs-Wolfe effect 5\To/Tq = $ig/3. Following the same steps, we may easily obtain its full non-linear 
generalization. 

Let us first relate the photon energy density p 7 to the energy density of the non-relativistic matter 
~p m by using the adiabaticity conditon. Again here a bar indicates that we are considering quantities in 
the locally homogeneous Universe described by the metric ([6]). Using the energy continuity equation on 
large scales dp/di = —2>H(j) + P), where H = dlna/dt and P is the pressure of the fluid, one can easily 
show that there exists a conserved quantity in time at any order in perturbation theory 



3 -> (p' + P' 



The perturbation 5 J- is a gauge-invariant quantity representing the non-linear extension of the curvature 
perturbation £ on uniform energy density hypersurfaces on superhorizon scales for adiabatic fluids |44j . 
Indeed, expanding it at first and second order one gets the corresponding definition £i = —ipi — 5\p/p 
and the quantity £2 introduced in Ref. [55]. We will come back to these definitions later. At first order 
the adiabaticity condition corresponds to set (j 7 = Cim f° r the curvature perturbations relative to each 
component. At the non-linear level the adiabaticity condition generalizes to 

1 f dp m 1 f dp 7 



(10) 

3 J Pm 4 J p^ 

or 

lnp m = lnp 7 / 4 . (11) 

To make contact with the standard second-order result, we may expand in Eq. (jlip the photon energy 
density perturbations as <5p 7 //0 7 = 5\py/p^ + ^2P 7 //0 7 , and similarly for the matter component. We 
immediately recover the adiabaticity condition 

8 2 p 1 _ 4 5 2 p m | 4 / 5 1 p m \ 2 ^ 

P 7 3 p m 9 V Pm J 
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given in Ref. [5]. 

Next we need to relate the photon energy density to the gravitational potentials at the non-linear level. 
The energy constraint inferred from the (0-0) component of Einstein equations in the matter-dominated 
era with the "barred" metric Q is 

—2 8ttG N _ 

H = -^—Pm ■ (13) 

Using Eqs. (j3J) and ([5]) the Hubble parameter H reads 

-ff=~ = e-*'(H-* t ), (14) 
a at 

where H = dlna/dt is the Hubble parameter in the "unbarred" metric. Eq. thus yields an expression 
for the energy density of the non-relativistic matter which is fully nonlinear, being expressed in terms of 
the gravitational potential 

Pm = Pm.e' 2 ^ , (15) 

where we have dropped ^ which is negligible on large scales. By perturbing the expression (|15|) we are 
able to recover in a straightforward way the solutions of the (0-0) component of Einstein equations for 
a matter-dominated Universe in the large-scale limit obtained at second-order in perturbation theory. 
Indeed, recalling that is perturbatively related to the quantity cp = (pi 1 ) + (pi 2 ) /2 used in Ref. [5] by 
(pi 1 ) = (pi 1 ) and (pi 2 ) = (pi 2 ) — 2{<pi 1 )) 2 , one immediately obtains 



Pm 
1 hPm 



-0(2) +4 (0W) 2 . (16) 

The expression for the intrinsic temperature of the photons at the last scattering surface oc p\ 



2 p ri 

c ten 

follows from Eqs. (jlip and (115[> 



T e = T e e" 2 *^ 3 . (17) 

Plugging Eqs. (jHJ) and (JT7J) into the expression ((7J) we are finally able to provide the expression for the 
CMB temperature which is fully nonlinear and takes into account both the gravitational redshift of the 
photons due to the metric perturbations at last scattering and the intrinsic temperature anisotropies 

'^W e *i/3 

From Eq. (|18|1 we read the non-perturbative anisotropy corresponding to the Sachs- Wolfe effect 



T ° = ^) T£e ' (18) 



= e $,/3 _ 1 (lg) 

To 

Eq. (|19p is one of the main results of this paper and represents at any order in perturbation theory the 
extension of the linear Sachs- Wolfe effect. At first order one gets 



T r 



(20) 
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and at second order 

This result shows that the CMB anisotropies is nonlinear on large scales and that a source of NG is 
inevitably sourced by gravity. 



3 Perturbing gravity 



In this Section we provide the necessary tools to deal with perturbed gravity, giving the expressions for 
the Einstein tensor perturbed up to second-order around a flat Friedmann- Robertson- Walker background, 
and the relevant Einstein equations. In the following we will adopt the Poisson gauge which eliminates 
one scalar degree of freedom from the goi component of the metric and one scalar and two vector degrees 
of freedom from gij . We rewrite the metric ([I]) as 

ds 2 = a 2 {rj) l-e^dr] 2 + 2u i dx i dr ] + (e~ 2 %- + Xij^dxi] , (22) 

where a{rj) is the scale factor as a function of the conformal time rj. As we previously mentioned, for 
the vector and tensor perturbations, we will neglect linear vector modes since they are not produced in 
standard mechanisms for the generation of cosmological perturbations (as inflation), and we also neglect 
tensor modes at linear order, since they give a negligible contribution to second-order perturbations. 
Therefore we take Wj and Xij to be second-order vector and tensor perturbations of the metric. 

Let us now give our definitions for the connection coefficients and their expressions for the metric (|TJ) . 
The number of spatial dimensions is n = 3. Greek indices (a, /3, fj,, u, ....) run from to 3, while 
latin indices (a,b, ...,i,j,k, ....m,n...) run from 1 to 3. The total spacetime metric has signature 
(—,+,+,+). The connection coefficients are defined as 

The Riemann tensor is defined as 

jya -pa -pa . pa pA pa pA {<~>A\ 

n pnv — 1 0u,n ~ 1 0n,v 1 Xp L pv ~ 1 \u l 0n ■ V Z4 J 

The Ricci tensor is a contraction of the Riemann tensor 

RfjLU = R^au ) (25) 
and in terms of the connection coefficient it is given by 

Rpu = d a r ^ u — T ua + T aa F^ u — F au T^ a . (26) 
The Ricci scalar is given by contracting the Ricci tensor 

R = R^. (27) 

The Einstein tensor is defined as 

Gp,u — Rpu — —g^uR ■ (28) 
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3.1 The connection coefficients 

For the connection coefficients we find 

r° = H + 

r* ro = J + nJ + c 2*+2*|i , 

= («-*<) % + i*<- + I(||-§|), 

(29) 



3.2 Einstein equations 

The Einstein equations are written as = k 2 T^ u , so that k 2 = 8-7tGn, where Gn is the usual Newtonian 
gravitational constant. They read 



G° 



a" 
=2* r 



3H 2 - m^' + 3(^') 2 - e 2 ® +2 * (diV&V - 2V 2 ^ 

1 



K T 



I 



V V + ( - 2 



2a 2 



a"\ w 



a J a 



2rpi 

2 ' 



-2$ 



H -2 - 3(*') + ZH (& + 2^') + 2^' 

a 



+ e 2 * (^.$8*$ + V z $ - V*tf 
e 2 * / • 

U 



% (ff Uj + - — 2 (d^ + ay) + - 2 [ u x ) + 5x5" - 5 V V 



k 2 TV 



(30) 
(31) 



(32) 



Taking the traceless part of Eq. (f32j) , we find 



(33) 



where Q is defined by V 2 Q = -P + 3iV, with P = P i h 



(34) 
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and V 2 N = d^P^. The trace of Eq. ((32J gives 



e 2 * / 

H (2<9 fc <I><9 fc <I> + d k $>d k ty - 2d k <S>d k ^ + 2(P - 3iV) 

^Vr*. (35) 



3 

From Eq. (|3ip . we may deduce an equation for w J 

I-72, ,i 



-VV + (4?T - 2— J u/ 1 (36) 

= - - ^Jj (2 (d»V + (ft - 9?'*) - STrCWe" 2 *^ 

Here T^, is the energy momentum tensor accounting for different components, photons, baryons, dark 
matter. We will give the expressions later for each component in terms of the distribution functions. 



4 Setting the initial conditions from the primordial NG 

In this review we are concerned with the second-order evolution of the cosmological perturbations. This 
requires that we well define the initial conditions of the cosmological perturbations at second order. These 
boundary conditions may or may not contain already some level of NG. It they do, we say that there 
exist some primordial NG. The latter is usually defined in the epoch in which the comoving curvature 
perturbation remains constant on large superhorizon scales. In the standard single-field inflationary 
model, the first seeds of density fluctuations are generated on super-horizon scales from the fluctuations 
of a scalar field, the inflaton pp. 

In order to follow the evolution on super-horizon scales of the density fluctuations coming from the 
various mechanisms, we use the curvature perturbation of uniform density hypersurfaces C = + 
£( 2 )/2 + where = —ip^ — HSp^/p' and the expression for £( 2 ) is given by |45] 

C (2) = _^(2)_^^ + AC ( 2)) (37) 



with 



AC (2) = 2% — j-—r 
P P 

+ 2^V 1 )' + 2^ 
P 

/ 5 (i) r \ 2 / j, 



P' I \ P 



P) ( n ^- n '-2U 2 )+2^ 2 . (38) 
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The crucial point is that the gauge-invariant curvature perturbation £ remains constant on super-horizon 
scales after it has been generated during a primordial epoch and possible isocurvature perturbations are 
no longer present. Therefore, we may set the initial conditions at the time when £ becomes constant. 
In particular, f^ 2 -* provides the necessary information about the "primordial" level of non-Gaussianity 
generated either during inflation, as in the standard scenario, or immediately after it, as in the curvaton 
scenario. Different scenarios are characterized by different values of <^ 2 ). For example, in the standard 

single-field inflationary model = 2 (C^) + O (e, 77) [6], where e and rj are the standard slow-roll 
parameters [I]. In general, we may parametrize the primordial non-Gaussianity level in terms of the 
conserved curvature perturbation as in Ref. [26j 



C (2) = 2a NL (c«) 2 , (39) 



where the parameter onl depends on the physics of a given scenario. For example in the standard scenario 
onl — 1> while in the curvaton case <znl = (3/4r) — r/2, where r (/v/p)d is the relative curvaton 
contribution to the total energy density at curvaton decay [Ml [5] . Alternatives to the curvaton model are 
those models characterized by the curvature perturbation being generated by an inhomogeneity in the 
decay rate j9[ [10] or the mass [TT] or of the particles responsible for the reheating after inflation. Other 
opportunities for generating the curvature perturbation occur at the end of inflation [13] and during 
preheating |14] . All these models generate a level of NG which is local as the NG part of the primordial 
curvature perturbation is a local function of the Gaussian part, being generated on superhorizon scales. 
In momentum space, the three point function, or bispectrum, arising from the local NG is dominated by 
the so-called "squeezed" configuration, where one of the momenta is much smaller than the other two. 
Other models, such as DBI inflation [16] and ghost inflation [JTj, predict a different kind of primordial 
NG, called "equilateral", because the three-point function for this kind of NG is peaked on equilateral 
configurations, in which the lenghts of the three wavevectors forming a triangle in Fourier space are equal 

HE]. 

One of the best tools to detect or constrain the primordial large-scale non-Gaussianity is through the 
analysis of the CMB anisotropies, for example by studying the bispectrum [5]. In that case the standard 
procedure is to introduce the primordial non-linearity parameter /nl characterizing the primordial non- 
Gaussianity via the curvature perturbation [5^1 

C = C (1) + ^/nl(c (1) *C (1) ) , (40) 

where the coefficient 3/5 arises from the first-order relation connecting the comoving curvature pertur- 
bation and the gravitational potential, £W = —5(f>^/3, and the ^-product reminds us that one has to 
perform a convolution product in momentum space and that /nl is indeed momentum-dependent. To 
give the feeling of the resulting size of /nl when |onl| S> 1> /nl — 5onl/3 (see Refs. [5J [26]). Present 
limits on NG are summarized by -4 < < 80 [48] and -125 < fZT a < 435 [49] at 95% CL, where 



/NL ^ ou P^l miu ^ JNL 

/j^l and /nl"' stand for the non-linear parameter in the case in which the squeezed and the equilateral 
configurations dominate, respectively. 



2 The sign covention is such that it follows the WMAP covention for /nl, see Sec. [5] 



14 



5 CMB anisotropics at second-order on large scales 



In this section, we provide the exact expression for large-scale CMB temperature fluctuations at second 
order in perturbation theory. What this Section contains, therefore, should be considered as a more 
technical elaboration of what the reader can find in Section 2. The final expression we will find has 
various virtues. First, it is gauge-invariant. Second, from it one can unambiguously extract the exact 
definition of the nonlinearity parameter /nl which is used by the experimental collaborations to pin 
down the level of NG in the temperature fluctuations. Third, it contains a "primordial" term encoding 
all the information about the NG generated in primordial epochs, namely during or immediately after 
inflation, and depends upon the various fluctuation generation mechanisms. As such, the expression 
neatly disentangles the primordial contribution to the NG from that arising after inflation. Finally, the 
expression applies to all scenarios for the generation of cosmological perturbations. 

In order to obtain our gauge-independent formula for the temperature anisotropies, we again perturb 
the spatially flat Robertson- Walker background. We expand the metric perturbations (I22D in a first and 
a second-order part as 



ffoo = -a 

.2 



9ij 



a~ 



2 (l + 2^ + + 2 (0«) 2 ) , 90i = a 2 ( w « + \JP 
(1 - - + 2 (V«) 2 )^ + (xff + \xf) 



(41) 

Again, the functions cj)^ , uj[' , ip W and Xij\ where (r) = (1,2), stand for the rth-order perturbations 
of the metric. We can split oj^ = c^uA-* + where uW is the scalar part and uj^ is a transverse 
vector, i.e. d^cj^' = 0. The metric perturbations will transform according to an infinitesimal change of 
coordinates. From now on we limit ourselves to a second-order time shift 

V V ~ a W + ^{a' {1) a {1) - a {2) ) , (42) 

where a prime denotes differentiation w.r.t. conformal time. In general a gauge corresponds to a choice 
of coordinates defining a slicing of spacetime into hypersurfaces (at fixed time rj) and a threading into 
lines (corresponding to fixed spatial coordinates x), but in this Section only the former is relevant so that 
gauge-invariant can be taken to mean independent of the slicing. For example, under the time shift, the 
first-order spatial curvature perturbation ip^ transforms as ip^ — > ip^ — % am (here % = a' /a), while 

<fi^ — > (jft-} + a',^ + HcS l \ u)j^ — > uf' — diCt( l \ and the traceless part of the spatial metric x\j turns 
out to be gauge-invariant. At second order in the perturbations we just give some useful examples like 



the transformation of the energy density and the curvature perturbation [41] . 

S&) p 5® p + p ' a{2)+a{l) (p" a(1) +p'a{ 1) + 2S® p') (43) 

and 



^(2) ^ ^(2) + 2a (1) ^ (1)/ - Ua 2 (1) - Hama'm 



2u\ 1) -af r] )af-Ua {2) . (44) 
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We now construct in a gauge-invariant way temperature anisotropics at second order. Temperature 
anisotropies beyond the linear regime have been calculated in Refs. [241 141] . following the photons path 
from last-scattering to the observer in terms of perturbed geodesies. The linear temperature anisotropies 
read [Ml E] 

^ = 4>-4» ei+ TP- f'dXAW, (45) 



where = ifi^ + <p^ + u^e 1 — \xtj' 1 e J ' , the subscript £ indicates that quantities are evaluated 
at last-scattering, e l is a spatial unit vector specifying the direction of observation and the integral is 

evaluated along the line-of-sight parametrized by the affine parameter A. Eq. (|45p includes the intrinsic 

(l)i 

fractional temperature fluctuation at emission T£, the Doppler effect due to emitter's velocity v £ and 
the gravitational redshift of photons, including the Integrated Sachs- Wolfe (ISW) effect. We omitted 
monopoles due to the observer O {e.g. the gravitational potential i^q evaluated at the event of obser- 
vation), which, being independent of the angular coordinate, can be always recast into the definition of 
temperature anisotropies. Notice however that the physical meaning of each contribution in Eq. (j45|) is 
not gauge-invariant, as the different terms are gauge-dependent. However, it is easy to show that the 
whole expression ([45 p is gauge- invariant. Since the temperature T is a scalar, the intrinsic temperature 
fluctuation transforms as Tg — > + (T' ' /T)at\) = t^P — T-La^, having used the fact that the temper- 
ature scales as T a a" 1 . Notice, instead, that the velocity Vg does not change. Therefore, using the 
transformations of metric perturbations we find 



ATW ATW , fvs da' {1) AT( 1} 



(46) 



where we have used the fact that the integral is evaluated along the line-of-sight which can be parametrized 
by the background geodesies = (A, (Xq — Ag)e i ) (with dX/d-q = 1), and the decomposition for the 
total derivative along the path for a generic function /(A, x l (X)), f' = jj\ = §\ + ®if e Eq. ([^6]) shows 
that the expression (j45p for first-order temperature anisotropies is indeed gauge-invariant (up to monopole 
terms related to the observer O). Temperature anisotropies can be easily written in terms of particular 
combinations of perturbations which are manifestly gauge- invariant. For the gravitational potentials we 
consider the gauge-invariant definitions tpQ^ = ip^ — Hu)^ and 4>qj = + Hlj^ + a/ 1 )'. For the 
(0 — i) component of the metric and the traceless part of the spatial metric we define = w| 

and x[^ G1 = Xi^ ■ For the matter variables we use a gauge-invariant intrinsic temperature fluctuation 

Tqj = r^ 1 ) — Huj^, while the velocity itself is gauge-invariant Vq^ = v^ 1 under time shifts. Following 
the same steps leading to Eq. (j46p one gets the linear temperature anisotropies in Eq. (|45p m terms of 
these gauge-invariant quantities 



where = <f>^ + ip^ + u^ Gl ei — ^Xi]^ Glele ^ an d we omitted the subscript £. For the primordial 
fluctuations we are interested in the large-scale modes set by the curvature perturbation i^ 1 ). Defining 
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a gauge-invariant density perturbation S^pqi = S^'p + p'uj^ l \ we write the curvature perturbation 

Since for adiabatic perturbations in the radiation (7) and matter 
(to) eras (1/4) (JW p 7 /p 7 ) = (l/3)(5^ p m /p m ), one can write the intrinsic temperature fluctuation as 

r^ 1 ) = (l/4)(5^ 1 )p 7 /p 7 ) = —%{5^p/p') and a gauge-invariant definition is T<~jj = —'H{5 <yl ^pQi/p'). In the 

large-scale limit, from Einstein equations, in the matter era </>qj = V>gi = ~~ |Cgi • Thus we obtain the 

AT^^ (I s ) (11 (1) 

large-scale limit of temperature anisotropies (|4"7j) — = 2iPq{ + Cqi = V'gi/^ i- e - * ne usua l Sachs- Wolfe 
effect. 

At second order, the procedure is similar to the one described so long, though more lengthy and 
cumbersome. We only provide the reader with the main steps to get the final expression. The second- 
order temperature fluctuations in terms of metric perturbations read \2A\ |4"T] 

AT( 2 ) 1 (2) 1 (2)< 1 (2) 

—f— = 2^£ 2 £ i + 2 £ ^ £ > 

+ (li(Ae) + v^et) - rP + + )) 



+ + (4 1)J + x«V) - t,ge* + r« 

2 V £i v £ T £ + Q d i a v e e ^e 

(!) ( , ,(1)* 



+ ; -PiM) ■ (48) 

Here / 2 is the second-order ISW 

/ 2 (A,) = f X£ d\£AW - (cW' - X g ) V)(JfcW < + e'fcCDO) + 2fc«°,4«' 

JA 2 

+ 2^ (1) 'A( 1 ) + X W°AW" + xW^'] , (49) 

where = ^ (2) + ^ {2) + u^e* - ^xJjW, while fcW°(A) = -2^ - wW'ej + Ii(A) and fcW^A) = 
— 2</>( 1 ^e* — a/ 1 )* + x^ l ^ e j ~ I{(X) are the photon wave vectors, withli(A) given by the integral in Eq. (|45p 
and /{(A) is obtained from the same integral replacing the time derivative with a spatial gradient. Finally 
in Eq. (@8D 

X W°(A) = / A dA' [-2<-/> (1) - wfV + (A - A')^ (1) 'l (50) 



and 



x (1)l (A) = - / A dA' [2^ (1) e* + - x (1)ii ej + (A - AOAW'*] (51) 

are the geodesies at first order, and S- 1 ^ 1 = e l — is the direction of the photon emission. As usual 

we have omitted the monopole terms due to the observer. Using the transformation rules of Ref. |41j . 
it is possible to check that the expression ([48 j) is gauge-invariant. We can express the second-order 
anisotropies in terms of explicitly gauge-invariant quantities, whose definition proceeds as for the linear 
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case, by choosing the shifts such that u)^ r > = 0. For example, we consider the gauge-invariant 
gravitational potential [25] 



(2) 



(1) 



a J a 
1 



+ («' + 2^ 2 ) W W] + 2u,«' foW + + - (aa^) 



(52) 



where a( 2 ) = uj^ + w^W 1 )' + V~ 2 <9*[— Aip^ diUi^ — 2u^' diU^]. Expressing the second-order tem- 
perature anisotropies (|48|) in terms of our gauge-invariant quantities and taking the large-scale limit we 
find 



AT§/T = (l/2)*g> - (1/2) + (1/2)4? + , 



(53) 



(having dropped the subscript £), and the gauge-invariant intrinsic temperature fluctuation at emission 
is 



(54) 



We have dropped those terms which represent integrated contributions and other second-order small-scale 
effects that can be distinguished from the large-scale part through their peculiar scale dependence. At 
this point we make use of Einstein's equations. We take the expression for C (2) in Eqs. (|57jt and (15511. 
and we use the (0 — 0) component and the traceless part of the (i — j) Einstein's equations (|30|) and (|32l) 
after having appropriately expanded the exponentials. Thus, on large scales we find that the temperature 
anisotropies are given by 



AT, 



(2) 



GI 
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K_ 
10 



1 

To 



^GI 



where we have defined a kernel 



K = WV- A did j (d l ^dj^) - V~ 2 (— d^d^) . 



(55) 



(56) 



Eq. (|55p is the main result of this Section. It clearly shows that there are two contributions to the final 
nonlinearity in the large-scale temperature anisotropies. The contribution, — 2(£qj) 2 ], comes from 
the "primordial" conditions set during or after inflation. They are encoded in the curvature perturbation 
C which remains constant once it has been generated. The remaining part of Eq. (|55p describes the post- 
inflation processing of the primordial non-Gaussian signal due to the nonlinear gravitational dynamics, 
including also second-order corrections at last scattering to the Sachs- Wolfe effect [Ml HI]- Thus, the 
expression in Eq. (155p allows to neatly disentangle the primordial contribution to NG from that coming 
from that arising after inflation and from it we can extract the nonlinearity parameter /nl, see the 
expression (I40p which is usually adopted to phenomenologically parametrize the NG level of cosmological 
perturbations and has become the standard quantity to be observationally The definition of /nl adopted 
in the analyses performed in Refs. [35] goes through the conventional Sachs- Wolfe formula AT/T = — <3>/3 
where $ is Bardeen's potential, which is conventionally expanded as (up to a constant offset, which only 
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affects the temperature monopole) $ = $l + /nl (^l) * (^l)- Therefore, using = — fV'Gi during 
matter domination, from Eq. (I55D we read the nonlinearity parameter in momentum space 



/ NL (ki,k 2 ) 

where 



5 n \ j. 1 3 r 
3 (1 - *nl) + 6 " lo* 



(57) 



/C = 10 (ki • k 3 )(k 2 • k 3 )/£; 4 - H kl • k 2 /A: 2 , (58) 



with k 3 + ki + k2 = and fe = | k 3 1 . To obtain Eq. (j57|) we have made use of the expression (|39j) to set the 
initial conditions. In particular within the standard scenario where cosmological perturbations are due 
to the inflaton the primordial contribution to NG is given by onl = 1 — i( n C — 1) [S] , where the spectral 
index is expressed in terms of the usual slow-roll parameters as — 1 = — 6e + 2r\ p]. The nonlinearity 
parameter from inflation now reads 

A£ = -^(»c-i)-| + ^a:. m 

Therefore the main NG contribution comes from the post-inflation evolution of the second-order per- 
turbations which give rise to order-one coefficients, while the primordial contribution is proportional to 
|n£ — 1| <C 1. This is true even in the "squeezed" limit first discussed by Maldacena in [6], where one of 
the wavenumbers is much smaller than the other two, e.g. k\ <C &2,3 and /C — > 0. 



6 CMB anisotropics at second-order at all scales 

As we already mentioned in the Introduction, despite the importance of NG in CMB anisotropics, little 
effort has ben made so far to provide accurate theoretical predictions of it. On the contrary, the vast 
majority of the literature has been devoted to the computation of the bispectrum of either the comovig 
curvature perturbation or the gravitational potential on large scales within given inflationary models. 
These, however, are not the physical quantities which are observed. One should instead provide a full 
prediction for the second-order radiation transfer function. A preliminary step towards this goal has been 
taken in Refs. (23j [24l EE] (see also [29]) where the full second-order radiation transfer function for the 
CMB anisotropies on large angular scales in a flat universe filled with matter and cosmological constant 
was computed, including the second-order generalization of the Sachs- Wolfe effect, both the early and 
late Integrated Sachs- Wolfe (ISW) effects and the contribution of the second-order tensor modeso We 
have partly reported about these works in the previous Sections. 

In this Section we wish to offer a summary of some of the second-order effects in the CMB anisotropies 
on small scales. There are many sources of NG in CMB anisotropies, beyond the primordial one. The 
most relevant sources are the so-called secondary anisotropies, which arise after the last scattering epoch. 
These anisotropies can be divided into two categories: scattering secondaries, when the CMB photons 
scatter with electrons along the line of sight, and gravitational secondaries when effects are mediated 
by gravity [31]. Among the scattering secondaries we may list the thermal Sunyaev-Zeldovich effect, 
where hot electrons in clusters transfer energy to the CMB photons, the kinetic Sunyaev-Zeldovich effect 

3 For some recent papers focusing on the generation and evolution of tensor perturbations at second-order see, 
e.g., Refs. [501 EU [52] . 
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produced by the bulk motion of the electrons in clusters, the Ostriker-Vishniac effect, produced by 
bulk motions modulated by linear density perturbations, and effects due to reionization processes. The 
scattering secondaries are most significant on small angular scales as density inhomogeneities, bulk and 
thermal motions grow and become sizeable on small length-scales when structure formation proceeds. 

Gravitational secondaries arise from the change in energy of photons when the gravitational potential 
is time-dependent, the ISW effect, and gravitational lensing. At late times, when the Universe becomes 
dominated by the dark energy, the gravitational potential on linear scales starts to decay, causing the 
ISW effect mainly on large angular scales. Other secondaries that result from a time dependent potential 
are the Rees-Sciama effect, produced during the matter-dominated epoch by the time evolution of the 
potential on non-linear scales. 

The fact that the potential never grows appreciably means that most second order effects created by 
gravitational secondaries are generically small compared to those created by scattering ones. However, 
when a photon propagates from the last scattering to us, its path may be deflected because of the gravi- 
tational lensing. This effect does not create anisotropies, but only modifies existing ones. Since photons 
with large wavenumbers k are lensed over many regions (~ k/H, where H is the Hubble rate) along the 
line of sight, the corresponding second-order effect may be sizeable. The three-point function arising from 
the correlation of the gravitational lensing and ISW effects generated by the matter distribution along 
the line of sight [32j[33] and the Sunyaev-Zeldovich effect [M] are large and detectable by Planck [35] . 

Another relevant source of NG comes from the physics operating at the recombination. A naive esti- 
mate would tell that these non-linearities are tiny being suppressed by an extra power of the gravitational 
potential. However, the dynamics at recombination is quite involved because all the non-linearities in the 
evolution of the baryon-photon fluid at recombination and the ones coming from general relativity should 
be accounted for. This complicated dynamics might lead to unexpected suppressions or enhancements of 
the NG at recombination. A step towards the evaluation of the three-point correlation function has been 
taken in Ref. [37] where some effects were taken into account in the in so-called squeezed triangle limit, 
corresponding to the case when one wavenumber is much smaller than the other two and was outside the 
horizon at recombination. 

This Section, which is based on Refs. [381 139] . present the computation of the full system of Boltz- 
mann equations, describing the evolution of the photon, baryon and Cold Dark Matter (CDM) fluids, 
at second order and neglecting polarization, These equations allow to follow the time evolution of the 
CMB anisotropies at second order on all angular scales from the early epochs, when the cosmological 
perturbations were generated, to the present time, through the recombination era. These calculations 
set the stage for the computation of the full second-order radiation transfer function at all scales and 
for a a generic set of initial conditions specifying the level of primordial non-Gaussianity. Of course on 
small angular scales, fully non-linear calculations of specific effects like Sunyaev-Zel'dovich, gravitational 
lensing, etc. would provide a more accurate estimate of the resulting CMB anisotropy, however, as 
long as the leading contribution to second-order statistics like the bispectrum is concerned, second-order 
perturbation theory suffices. 

6.1 The collisionless Boltzmann equation for photons 

We are interested in the anisotropies in the cosmic distribution of photons and inhomogeneities in the 
matter. Photons are affected by gravity and by Compton scattering with free electrons. The latter are 
tightly coupled to protons. Both are, of course, affected by gravity. The metric which determines the 
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gravitational forces is influenced by all these components plus CDM (and neutrinos). Our plan is to write 
down Boltzmann equations for the phase-space distributions of each species in the Universe. 

The phase-space distribution of particles g(x l , P^, rj) is a function of spatial coordinates x l , conformal 
time rj, and momentum of the particle = dx^/dX where A parametrizes the particle path. Through 
the constraint P 2 = g flu P >1 P 1 ' = —m 2 , where m is the mass of the particle one can eliminate P° and 
g{x l ,P 3 ,rj) gives the number of particles in the differential phase-space volume dx 1 dx 2 dx 3 dP 1 dP 2 dP 3 . 
In the following we will denote the distribution function for photons with /. 

The photons' distribution evolves according to the Boltzmann equation 



(60) 



where the collision term is due to the scattering of photons off free electrons. In the following we will 
derive the left-hand side of Eq. (|60|) while in the next section we will compute the collision term. 

For photons we can impose P 2 = g^ w P^P u = and using the metric ([1]) in the conformal time rj we 
find 



e 2<& ( pU ) , + £_ + 2 ^pUp 4 



:>0\2 



P 



0, 



where we define 



From the constraint ([6T 



P 



9ij 



iiP l P J • 



1/2 



(61) 
(62) 

(63) 



Notice that we immediately recover the usual zero and first-order relations P° = p/a and P° = p(l — 
$«)/«. 

The components P l are proportional to pn l , where ri 1 is a unit vector with n l rii = Sijtfn 3 = 1. We 
can write P l 



Cn l , where C is determined by 

g^pi = C 2 a 2 (e- 2 * + Xy nV) = p 2 



so that 



P' 



P 



I, \ 



'1/2 



P i * 

—ne 
a 



1 



(64) 



(65) 



where the last equality holds up to second order in the perturbations. Again we recover the zero and 
first-order relations P l = pn l /a and P % = pn l (l + VP^J/a respectively. Thus up to second order we can 
write 

30 _ -*P 



a 



1 + uji n' 



(66) 



Eq. ()65p and (j66|) allow us to replace P° and P % in terms of the variables p and ri 1 . Therefore, as it is 
standard in the literature, from now on we will consider the phase-space distribution / as a function of 
the momentum p = pri 1 with magnitude p and angular direction ri 1 , f = f(x l ,p,n l ,r]). 

Thus, in terms of these variables, the total time derivative of the distribution function reads 



df df df dx % df dp df dn l 
dr/ drj dx l dr) dp dr\ dri dr\ 



(67) 
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In the following we will compute dx 1 /drj, dp/drj and dn l /dr]. 



a) dx 1 jdr\: 
From 



and from Eq. ([65]) and ([66 



dx 1 dx 1 drj dx 1 ^ 
dX drj dX drj 



(68) 



dx^_ 
drj 



1 



fc„m 



(69) 



b) dp/ drj: 

For dp/drj we make use of the time component of the geodesic equation dP°/dX = — r^P a P^, where 
d/dX = {drj/dX) d/drf = P° d/dq, and 



dP° 



pa p/3 



Using the metric ([!]) we find 



2T° ap P a pP 



,0i> 



drj Q/3 P° 

. dg ua daf: 
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papfi 



2(ft + ( P° ) " + 4$ ;P°P l + muj i P°P i 



+ 2e 



-2* 



1 , 



(70) 



(71) 



On the other hand the expression (166p of P in terms of p and n* gives 



dP° 
drj 



pd&_ 

a drj 



(l + Win*) + e"* (l + n<) 



dip /a) 
drj 



P^ djuij rt) 
a drj 



Thus Eq. ([70]) allows us express dp/drj as 



p ar/ 2 J 



(72) 



(73) 



where in Eq. (]7ip we have replaced P° and P* by Eqs. (|66|) and (]65p . Notice that in order to obtain 
Eq. (j73p we have used the following expressions for the total time derivatives of the metric perturbations 



drj 



<9$ d$ dx' 
drj dx 1 drj 

8$ 



~dr~j + (Jx* 71 * 6 ' 1 ~~ Uj ^ ~ 7:Xkmn n 



(74) 
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and 



dfarf) =n J dun | dujdxA = 3uH ni | dui n i n j (75) 
drj \ drj dxi drj J drj dx 3 

where we have taken into account that cjj is already a second-order perturbation so that we can neglect 
dn % jdrj which is at least a first order quantity, and we can take the zero-order expression in Eq. (I69p . 
dx l /drj = n\ In fact there is also an alternative expression for dp/ drj which turns out to be useful later 
and which can be obtained by applying once more Eq. f)T4f) 



1 dp „ , d& T , , . , i 1 , , ; „• , s 

= -U - ~r + & + *' - w X - -x . (76) 
p dr/ drj 2 J 



c) dn % /drj: 

We can proceed in a similar way to compute drt /drj. Notice that since in Eq. (|67p it multiplies df/dn l 
which is first order, we need only the first order perturbation of dn l /drj. We use the spatial components 
of the geodesic equations dP l /dX = —T l a gP a P^ written as 



d pi papfi 

For the right-hand side we find, up to second order, 



(77) 



9 



dg a y dgfjy _ dg_ 
dxP dx a 



dx u 



4(H - V)P l P u + 2 X k + ~ P P 



pap): 



5O r>k 



ox 1 



(PO) 



2 

- 4 



dx k 



pi pk 



8^ u 

-2—5 km P k P m 
ox 1 



2HJ5 k -(^ + ^-^ 
3 \ dx k dxi dxi 



pjpk 



(78) 



while the expression (|65[) of P % in terms of our variables p and ri 1 in the left-hand side of Eq. (|77[) brings 



dP { 
drj 



P ^ 

-e 
a 



d^ 

di] 



,■ ad(p/a) i 
n % H K -^—^n l + 



dn l 
dr] 



k^m 



1 - -jXkmn n 



P . 
—n e 

a 2 



9 1 d (xk m n k n m 



drj 



(79) 



Thus, using the expression for P % and (JBHj) for P° in Eq. (JTSJ), together with the previous result (|75|) . 
the geodesic equation (|77p gives the following expression drt /drj (valid up to first order) 



— = ($,i + $t)nV 
drj 



1" 



(80) 
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To proceed further we now expand the distribution function for photons around the zero-order value 
jf (°) which is that of a Bose-Einstein distribution 



/(°W = 2 * , (81) 

exp l5¥j ~ 1 

where T{rj) is the average (zero-order) temperature and the factor 2 comes from the spin degrees of 
photons. The perturbed distribution of photons will depend also on x 1 and on the propagation direction 
n l so as to account for inhomogeneities and anisotropies 

f(x\p, n\ n) = (p, rj) + /« {x\ p, n\ V ) + ±/< 2 > (x\p, n\ V ), (82) 

where we split the perturbation of the distribution function into a first and a second-order part. The 
Boltzmann equation up to second order can be written in a straightforward way by recalling that the 
total time derivative of a given i-th. perturbation, as e.g. df^ 1 /drj is at least a quantity of the i-th. order. 
Thus it is easy to realize, looking at Eq. (|67p . that the left-hand side of Boltzmann equation can be 
written up to second order as 

* = ^ + I^_ p ^A(*(i) + 
drj drj 2 drj dp drj \ 2 

dfW_d_ t (1) (1) i (2) i (2) 

^ P dp dr,\ ^ ^2 ^2 

df(°) du { { 1 d Xij t j 

- P^r~ -7T n ~nP^ -^-nn 3 , 83) 

op orj 2 op or] 

where for simplicity in Eq. (|83p we have already used the background Boltzmann equation (df /drj)\^ = 0. 
In Eq. (I83|) there are all the terms which will give rise to the integrated Sachs- Wolfe effects (corresponding 
to the terms which explicitly depend on the gravitational perturbations), while other effects, such as the 
gravitational lensing, are still hidden in the (second-order part) of the first term. In fact in order to 
obtain Eq. (j83j) we just need for the time being to know the expression for dp/ drj, Eq. ([76]) . 



6.2 Collision term 
6.2.1 The Collision Integral 

In this section we focus on the collision term due to Compton scattering 

e(q)7(p) <— > e(q / )7(p) , (84) 

where we have indicated the momentum of the photons and electrons involved in the collisions. The 
collision term will be important for small scale anisotropies and spectral distortions. The important 
point to compute the collision term is that for the epoch of interest very little energy is transferred. 
Therefore one can proceed by expanding the right hand side of Eq. (|60p both in the small perturbation, 
Eq. (I82j) . and in the small energy transfer. 
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The collision term is given (up to second order) by 

C(p) = C(p)ae* , (85) 

where a is the scale factor and0 

r( v J_ f dq dq!_ dp' 

' [P) E(p)J (2vr) 3 2£(q) (2vr) 3 2£(q') (2vr) 3 2£(p') 

x (2irY8\q + p-q' -p')\M\ 2 

x { 5 (q')/( P ')[l + /(P)] " S(q)/(P)[1 + f(p')}} (86) 

where -E(q) = (q 2 +m 2 ) 1 / 2 , M is the amplitude of the scattering process, 5 A {q+p — q'—p') = <5 3 (q+p — q' — 
p')5(E(q) + p — E(q') —p') ensures the energy-momentum conservation and g is the distribution function 
for electrons. The Pauli suppression factors (1 — g) have been dropped since for the epoch of interest the 
density of electrons n e is low. The electrons are kept in thermal equilibrium by Coulomb interactions 
with protons and they are non-relativistic, thus we can take a Maxwell-Boltzmann distribution around 
some bulk velocity v 

/ 2vr \ 3 / 2 f (q-m e v) 2 l , . 

g(q)= n J—-) exp - iq o ° \ (87) 



\m e T e J 2m e T e 

By using the three dimensional delta function the energy transfer is given by E(q) — i£(q + p — p') and 
it turns out to be small compared to the typical thermal energies 

E(q) - E(q + p - p') ~ (P ~ ^ q = 0(Tq/m e ) . (88) 

m e 

In Eq. (|88p we have used E(q) = m e + q 2 /2m e and the fact that, since the scattering is almost elastic 
(p ~ p'), (p — p') is of order p ~ T, with q much bigger than (p — p'). In general, the electron momentum 
has two contributions, the bulk velocity (q = m e v) and the thermal motion (q ~ {m e T) 1 / 2 ) and thus 
the parameter expansion q/m e includes the small bulk velocity v and the ratio (T /m e ) l l 2 which is small 
because the electrons are non-relativistic. 

The expansion of all the quantities entering the collision term in the energy transfer parameter and 
the integration over the momenta q and q' is described in details in Ref. [53]. It is easy to realize that 
we just need the scattering amplitude up to first order since at zero order g(q') = g(q + p — p') = <?(q) 
and 5(E(q) + p — E(q') — p') = 5(p — p') so that all the zero-order quantities multiplying \M\ 2 vanish. 
To first order 

\M\ 2 = 6vro- T m 2 [(l + cos 2 (9) - 2cos6»(l - cos6»)q • (p + p')/m e ] , (89) 

where cos6 = n • n' is the scattering angle and gt the Thompson cross-section. The resulting collision 
term up to second order is given by [53j 

Ap 

+ c { H(p,p') + c^(p,p')+c {2 \p,p') 



C(P) = -^J dpP ^ 



c W( P ,p') + cg ) (p,p , ) + 4 2) (p,p / ) 

(90) 



4 The reason why we write the collision term as in Eq. (j55j) is that the starting point of the Boltzmann equation 
requires differentiation with respect to an affine parameter A, df /d\ = C'. In moving to the conformal time rj one 
rewrites the Boltzmann equation as df/dij = C(P°) _1 , with P° — drj/dX given by Eq. (|66|) . Taking into account 
that the collision term is at least of first order, Eq. (|85|) then follows. 
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where we arrange the different contributions following Ref. [53]. The first order term reads 



c«(p,p') = (1 + cos 2 



S(p- P ')(f (1) (p')-f (1) (p)) 



+ (/ (0) (P')-/ (0) (P))(P 



p j • V- 



dp' 



(91) 



while the second-order terms have been separated into four parts. There is the so-called anisotropy 
suppression term 

4 2) (p> pO = \ (i + cos 2 e) s(p - p')(/ (2) (p) - / (2) (p)) ; (92) 

a term which depends on the second-order velocity perturbation defined by the expansion of the bulk 
flow as v = vW + v( 2 ) /2 



4 2) (P, P') = \{l + cos 2 0)(/(°)(p') - /W(p))(p - p') • v( 2 ) 



a set of terms coupling the photon perturbation to the velocity 

4W) = (/ (1 V)-/ (1) (p) 

d5(p — p') 



1 + cos 2 6>) (p - p') • v 



dp' 



2cos0(l - cos6)5(p - p')(n + n') • v 



and a set of source terms quadratic in the velocity 

<$(p,p') = (/ ( V)-/ (0) (p)) (p-pO-v 



1 + cos 2 6 



(p - p') " v d 2 5(p - p') 



2 dp' 2 
2cos#(l — cos#)(n + n') • v 



d5(p — p') 
dp' 



The last contribution are the Kompaneets terms describing spectral distortions to the CMB 



4W) 



(l + cosV (P " P,)2 



2m f 



[p')-f {0 \p))T, 



x 



d2 %/' ] ~ (f i0) (p') + f i0) (p) + 2/ (0) (p')/ (0) (p) 



<9<5(p — p') 
dp' 



+ 



2(p — p')cos#(l — cos^ 



x f^(p')(l + f^(p))(f {0) (p')-f i0) (p) 



5{p-p') 

d5(p — p') 
dp' 



(93) 



(94) 



(95) 
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Let us make a couple of comments about the various contributions to the collision term. First, notice the 
term c\, (p, p ) due to second-order perturbations in the velocity of electrons which is absent in Ref. [53]. 
In standard cosmological scenarios (like inflation) vector perturbations are not generated at linear order, 
so that linear velocities are irrotational v^ 1 = d t v^ 1 \ However at second order vector perturbations are 
generated after horizon crossing as non-linear combinations of primordial scalar modes. Thus we must 
take into account also a transverse (divergence-free) component, v^ 1 = d l v^ +Vj^ % with diV^ 1 = 0. As 
we will see such vector perturbations will break azimuthal symmetry of the collision term with respect 
to a given mode k, which instead usually holds at linear order. Secondly, notice that the number density 
of electrons appearing in Eq. (|90|) must be expanded as n e = n e (l + 5 e ) and then 

6McW(p,p') (96) 

gives rise to second-order contributions in addition to the list above, where we split <5 e = 5^+6^/2 
into a first- and second-order part. In particular the combination with the term proportional to v in 
c^(P'P') gives r i se to the so-called Vishniac effect, as discussed in Ref. |53j . 

6.2.2 Computation of different contributions to the collision term 

In the integral (|90p over the momentum p' the first-order term gives the usual collision term 



C«(p) 



n e a T 



/o (1) (p) + ^/ 2 (1) ^(v-n)-/«-p^v.n 



(97) 



where one uses the decomposition in Legendre polynomials 

/ (1) (x,p,n) =^(2£ + l)/, (1) (p)P,(cosi?), (98) 

where •& is the polar angle of n, cost? = n • v. 

In the following we compute the second-order collision term separately for the different contributions, 
using the notation C(p) = CW(p) + C^(p)/2. We have not reported the details of the calculation of the 
first-order term because for its second-order analog, (p, p') + c v (p, p') 5 the procedure is the same. 
The important difference is that the second-order velocity term includes a vector part, and this leads to 
a generic angular decomposition of the distribution function (for simplicity drop the time dependence) 



ft = HrV^ 1 / d^Yl m {n) . (100) 



/«(x, P ,n) = e E /i2(*.p)(-o' v^r y ^( n ) ■ (") 

l m=—£ ' 

such that 

47T 

Such a decomposition holds also in Fourier space. The notation at this stage is a bit confusing, so let us 
restate it: superscripts denote the order of the perturbation; the subscripts refer to the moments of the 
distribution. Indeed at first order one can drop the dependence on m setting m = using the fact that 
the distribution function does not depend on the azimuthal angle (p. In this case the relation with is 

fi% = H)- t (U + l)6 ma f< L) . (101) 
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a) c^(p,p'): 

The integral over p' yields 

c%\p) = 



3n e aT 



/<i' (P , P i 



3n e <TT 
Ap 



dp'p'5(p — p') 



/^[l + (n.n0 2 ][/( 2 )( P 0-/ (2) (p)] 



(102) 



To perform the angular integral we write the angular dependence on the scattering angle cos# = n • n' in 
terms of the Legendre polynomials 



[1 + (n • n'Y] 



l + ip 2 (n-n') 



1 + J E W«n" 



2£ + l 

m=—2 

where in the last step we used the addition theorem for spherical harmonics 

47T 



(103) 



P i = E ^(n)y/ m (n') 



(104) 



m=-2 



Using the decomposition ()100p and the orthonormality of the spherical harmonics we find 



C"a ) (p) = ""-e^T 



/S } (p)-/ (2) (p)-^ E ||/&)Wn 



m=-2 



(105) 



It is easy to recover the result for the corresponding first-order contribution in Eq. (|226p by using 
Eq. (p02|) . 

b)cf(p,p'): 

Let us fix for simplicity our coordinates such that the polar angle of n' is defined by // = v^ 2 ) • n' 
with (/)' the corresponding azimuthal angle. The contribution of c\; (p, p') to the collision term is then 

C ( 2)(p) = ^,( 2 )/ W[/ (o) (y) _ /( o )(p)] ^_P0 



4p 



We can use Eq. (1103P which in our coordinate system reads 



2 m 4l 2 (2 + m)! 21 ; 2 1 ; 



so that 



J ^P 2 (n • n') = P 2 (n • V ( 2 ))P 2 (n' • v( 2 )) = P 2 (/,)P 2 (//) 



(106) 



(107) 



(108) 
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By using the ortho-normality of the Legendre polynomials and integrating by parts over p' we find 

C® (p) = -n e cjt P^vM . n . (109) 

As it is clear by the presence of the scalar product v^ 2 ) p the final result is independent of the coordinates 

chosen. 

c) <S(p,p'): 

Let us consider the contribution from the first term 

4Vp' p') = i 1 + cos2 *) (/ (1 V) - / (1) (p)) (p - p') • • 

where the velocity has to be considered at first order. In the integral (|90p it brings 

3n e aTV 



1 r (2) 
2°Av(J) 



4p 



j dp'/-^ 1 f_ x d 4 [/ (1) (pO-/ (1) (p)] 



x (W-PV) r^(l + cos 2 0) 

JO 



(110) 



The procedure to do the integral is the same as above. We use the same relations as in Eqs. (|107[) 
and (|108p where now the angles are those taken with respect to the first-order velocity. This eliminates 
the integral over (f>', and integrating by parts over p' yields 



1 



(2) 



2 ^Av(I) 1 



3n e aTV f 1 dfj,' 
4p~~ 7_i T~ 



fill) 



pfr - 2/i')(/ (1 W) ~ + p\» ~ /* / )- 5/(1)(P ' / ' /) 



dp 

We now use the decomposition (|98[) and the orthonormality of the Legendre polynomials to find 

*HL,,'fV)( n „>\P„(,A Ma*'../, 



E 



(112) 



where we have used //^(aO-^KpO = -sP\{jJi!) + |P3(//), with P\{p!) = //. Thus from Eq. (jllip we get 



V (2) ^ 



?T, e cj T < v • n 



. .^(p) 



/ w (p)-/r(p)-p- 



+ V 



2/ 1 (1) (p)+p 
r(l)/ 



^ + i ft (v.„)f2/« W 



(113) 
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In c( 2 )(p, p') there is a second term 



J 2 ) 

-Av(II) 



-2cos#(l - cosfl) (/ (1) (p') - / {1) (P)) S(p - p')(n + n') • v, 



whose contribution to the collision term is 

1^(2) / \ 3n e a T v 



dp'p'8(p-p') £^(/W(p') 



/«(p))0u + //) /" 2 ^cos0(l-cos0). 
J 27r 



(114) 



This integration proceeds through the same steps as for C/\^j^(p)- I n particular by noting that cos#(l 
CO s6>) = -1/3 + Pi(cos0) - 2P 3 (cos6>)/3, Eqs. (fTUTD and IflOBjt allows to compute 



| ^cos#(l - cos9) = ~ + Pi(/i)Pi(/i') - ^P 2 (/i)P 2 ( / u') , 
and using the decomposition (f98l) we arrive at 

■ VT {vn/ 2 (1) (p)(l-P 2 (v.n)) 
ip 2 (v-n) (s/f'W-Sfb) 



(115) 



-r (2) m 



+ U 



(116) 



We then obtain 



/ (1) (P)-/ (1) (P)-P 



9/ (1) (p) # (D 



+ ^P 2 (v-n)(/ 2 (1) (p)-p 



5/ 2 (1) (p) 



dp 



+ v 



y[ l \p) 



(117) 



As far as the remaining terms, these have already been computed in Ref. [53] (see also Ref. [5l]) and 
here we just report them 

d)c$(p,p'): 

The term proportional to the velocity squared yield a contribution to the collision term 



i<7«(p) = Vr{(vn) : 



' dfW_ ii 2 a 2 /(°) 

^ <9p 20^ Sp 2 



+ V 



_3_ 2 £> 2 /(°) 

<9p 20^ <9p 2 



(118) 
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e) c£W): 

The terms responsible for the spectral distortions give 



r ^r+/ (0) (l+/ (0) ) 



(119) 



Finally, we write also the part of the collision term coming from Eq. ([96 



#>c«(p,p') ^ 8 ( pCU{p)=n e o T 8V 



+ 2/«p 2 (v. n )-/ 



(i) 



/^(p) 

d/(° ) 
° dp 



-v • n 



(120) 



6.2.3 Final expression for the collision term 

Summing all the terms we find the final expression for the collision term (|90p up to second order 



with 



and 



C7W(p) 



n e a T 



/d 1) (p) + ^2 1) ^(v-n)-/( 1 ) 



P" 



Op 



-v • n 



^ (2) (P) = Wr{\$<J>)-\ £ 5/2(P)Wn) 

I rre=— 2 



^/ (2) (P) + [/0 (1) (P) + • n) - /« 



1 Q/(o) 
2 P ~dp~ 



-(2) . 



n + v • n 



/ {1) (I 



5/ (o) 
p— — v • n 

/o (p) - P 3^ /a (?) + 2 P2 ( V ' n ) 



x /i 1) (p)-P 



d/ 2 (1) (p) * 

dp 



+ v 



-(l)/ N , ^/1 (1) (P) 

2/i ; (p) +P- 



9p 



+ 3 d/3 (1) (p) 

2^ dp 



v • n 



■ 0/(°> 11 2 9 2 / (0) 
P — 5 1" ™P 



dp 20 <9p 2 



df<® 
dp + 20 P dp 2 



P- 



3 2 a 2 / (0) 



m 2 <9p 



/ ( Te ^r +/(o)(i+/( ° )) 



(121) 



(122) 



(123) 
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Notice that there is an internal hierarchy, with terms which do not depend on the baryon velocity v, 
terms proportional to v • n and then to (v • n) 2 , v and v 2 (apart from the Kompaneets terms). In 
particular notice the term proportional to #e v • n is the one corresponding to the Vishniac effect. We 
point out that we have kept all the terms up to second order in the collision term. In Refs. [531 [54] many 

(2) 

terms coming from &^ have been dropped mainly because these terms are proportional to the photon 
distribution function /W which on very small scales (those of interest for reionization) is suppressed by 
the diffusion damping. Here we want to be completely general and we have to keep them. 



6.3 The Brightness equation 
6.3.1 First order 

The Boltzmann equation for photons is obtained by combining Eq. (|83j) with Eqs. f|122|) - fjl23j) . At first 
order the left-hand side reads 



d fW ^/(o) dx i d f(o) a^(i) 

P^, EI—J- +P- 



dr] 



dp dx 1 dr] dp dr] 



(124) 



At first order it is useful to characterize the perturbations to the Bose-Einstein distribution function ([81 
in terms of a perturbation to the temperature as 



exp 



P 



n -1 



1 



Thus it turns out that 



/ 



(i) 



op 



(125) 



(126) 



where we have used the fact that (9//90|e=o — —pdf^ /dp. In terms of this variable G^ 1 ) the linear 
collision term f|122j) will now become proportional to —pdf^ /dp which contains the only explicit de- 
pendence on p, and the same happens for the left-hand side, Eq. (|124p . This is telling us that at first 
order does not depend on p but only on x l ,n l ,r], 

0(1) = 9(l)( x i )n i )T ). 

This is well known and the 

physical reason is that at linear order there is no energy transfer in Compton collisions between photons 
and electrons. Therefore, the Boltzmann equation for B^ 1 ) reads 



39(1) 3 (i) 

+ n — — -. — I — ^-^n" 



dr] 
- n e aTa 



dx l 



dx 1 



dr] 







(!) + IeWp 2 ( v • n ) - 6 (1) + v • n 



(127) 



where we made us of = —pdf^/dpQ^\ according to the decomposition of Eq. (|98p . and we have 
taken the zero-order expressions for dx 1 /dr], dropping the contribution from dn l /dr] in Eq. (|67p since it 
is already first-order. 

Notice that, since is independent of p, it is equivalent to consider the quantity 



J dpp 3 fW 
JdpptfW ' 



(128) 
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being = 40W a t this order. The physical meaning of A^ is that of a fractional energy perturbation 
(in a given direction). From Eq. (J83]) another way to write an equation for A^ 1 ) - the so-called brightness 
equation - is 



d_ 

dr) 



n e aTa 



a{, 1} + -A^ J P 2 (v • n) - A (1) + 4v • n 



4 1. + W 
077 V 



(129) 



6.3.2 Second order 

The previous results show that at linear order the photon distribution function has a Planck spectrum 
with the temperature that at any point depends on the photon direction. At second order one could 
characterize the perturbed photon distribution function in a similar way as in Eq. (]125f) 



f(x\p,n\rj) 



T{ji) 



ex P i ^73 . e ~ 1 



(130) 



where by expanding = O^ 1 ) + O^ 2 ) /2 + ... as usual one recovers the first-order expression. For example, 
in terms of 0, the perturbation of /W is given by Eq. (|126p . while at second order 

?-f^- + i(^^) («")'■ 

However, as discussed in details in Refs. I53j . now the second-order perturbation will not be 
momentum independent because the collision term in the equation for does depend explicitly on 
p (defining the combination — (pdf^/dp)~ 1 f^ does not lead to a second-order momentum indepen- 
dent equation as above). Such dependence is evident, for example, in the terms of C^ 2 \p), Eq. (I123p . 
proportional to v or v 2 , and in the Kompaneets terms. The physical reason is that at the non- linear 
level photons and electrons do exchange energy during Compton collisions. As a consequence spectral 
distorsions are generated. For example, in the isotropic limit, only the Kompaneets terms survive giving 
rise to the Sunyaev-Zeldovich distorsions. As discussed in Ref. [54], the Sunyaev-Zeldovich distorsions 
can also be obtained with the correct coefficients by replacing the average over the direction electron 
(v 2 ) with the mean squared thermal velocity {v 2 h ) = 3T e /m e in Eq. ()123p . This is due simply to the 
fact that the distinction between thermal and bulk velocity of the electrons is just for convenience. This 
fact also shows that spectral distorsions due to the bulk flow (kinetic Sunyaev-Zeldovich) has the same 
form as the thermal effect. Thus spectral distorsions can be in general described by a global Compton 
y-parameter (see Ref. [53] for a full discussion of spectral distorsions). However in the following we will 
not be interested in the frequency dependence but only in the anisotropies of the radiation distribution. 
Therefore we can integrate over the momentum p and define [541 [53] 

as in Eq. ffT28|) . 
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Integration over p of Eqs. (|83p - (|123p is straightforward using the following relations 

-47V; 



J 3 df {0) 
dpp p- 



Op 



, ,9 2 / (0) 



dp 2 



Jdpp 3 f {1) = iVAW; J dpp 3 p^— = -ANA {1) . 



(133) 



Here N = J dpp 3 is the normalization factor (it is just proportional the background energy density 
of photons At first order one recovers Eq. (I129p . At second order we find 
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where we have expanded the angular dependence of A as in Eq. ([99 
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where we recall that the superscript stands by the order of the perturbation. When going to Fourier 
space some convolutions will appear and the coefficients A^(k', rj) are related to A^(k', rj) as A^(k') = 

iV47T(2*+l)y*n(k') 

A^(k'). In Eq. (|134p we have introduced the differential optical depth 

r 1 = —n e aTa . (137) 

It is understood that on the left-hand side of Eq. f|134|) one has to pick up for the total time derivatives 
only those terms which contribute to second order. Thus we have to take 



1 d 



A( 2 ) + 4$< 2 ) 



2dr] 
1 

2 \drj dx 



+ 



d 
dq 



AW+4^ 1 ) 



(2) 



(138) 



1 ( V + n iJL ) (a( 2 ) + 4$< 2 >) + n'($W + 



:ft(A (1) + 4$ (1) ) + ($JJ } + tfJjVV - + 

34 



gAW 

dri 1 



where we used Eqs. (|69|) and ([80]) . 



6.3.3 Hierarchy equations for multipole moments 

Let us now move to Fourier space. In the following, for a given wave- vector k we will choose the coordinate 
system such that e3 = k and the polar angle of the photon momentum is with /i = cos$ = k • n. Then 
Eq. (|134p can be written as 

A< 2 >' + ikfiA^ - t'A® = 5(k, n, 77) , (139) 

where S(\&,n,r]) can be easily read off Eq. (|134p . We now expand the temperature anisotropy in the 

(2) 

multipole moments A)' in order to obtain a system of coupled differential equations. By applying the 
angular integral of Eq. (| 136j) to Eq. (| 139j) we find 

" + r'Afl + S im (140) 
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where the first term on the right-hand side of Eq. (|140p has been obtained by using the relation 
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where Gi m = (— i) \J 4tt/(2£ + l)Y^ m (n) is the angular mode for the decomposition (|135p and = 
y/P — m 2 . This relation has been discussed in Refs. |56l [57] and corresponds to the term n % dA^> /dx l 
in Eq. (|134p . In Eq. (j!40p 5; m are the multipole moments of the source term according to the decom- 
position (|136p . We do not show its complete expression here, since it is very long, and we prefer to 
make some general comments (for some specific examples about the terms Si m see Sec. ITOT) . As it should 
be already clear from the calculations involving second-order perturbations performed so far, Eq. (|140p 
has the same functional form of the Boltzmann equation at linear order (just replace the order of the 
linear perturbations (1) with that of intrinsically second-order terms (2)) with the exception of the source 
terms 5/ m which now contain both intrinsic second-order perturbations and also products of first-order 
perturbations. Therefore as expected, at second order we recover some intrinsic effects which are char- 
acteristic of the linear regime. In Eq. (I140p the relation (|14ip represents the free streaming effect: when 
the radiation undergoes free-streaming, the inhomogeneities of the photon distribution are seen by the 
obsever as angular anisotropies. At first order it is responsible for the hierarchy of Boltzmann equations 
coupling the different £ modes, and it represents a projection effect of fluctuations on a scale k onto the 

( 2^ 

angular scale £ ~ krj. The term rA\' causes an exponential suppression of anisotropies in the absence 
of the source term Si m . The source term contains additional scattering processes and gravitational ef- 
fects. The intrinsically second-order part of the source term just reproduces the expression of the first 
order case. Of course the dynamics of the second-order metric and baryon-velocity perturbations which 
appear will be different and governed by the second-order Einstein equations and continuity equations. 
The remaining terms in the source are second-order effects generated as non-linear combinations of the 
primordial (first-order) perturbations. Notice in particular that they involve the first-order anisotropies 
and as a consequence such terms contribute to generate the hierarchy of equations (apart from the 
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free-streaming effect). On large scales (above the horizon at recombination) we can say that the main 
effects are due to gravity, and they include the Sachs- Wolfe and the (late and early) Sachs- Wolfe ef- 
fect due to the redshift photons suffer when travelling through the second-order gravitational potentials. 
These, toghether with the contribution due to the second-order tensor modes, have been already studied 
in details in Ref. [23j (see also [29] ) . Another important gravitational effect is that of lensing of photons 
as they travel from the last scattering surface to us. A contribution of this type is given by the last term 
of Eq. (|138p . On the other hand examples of second-order scattering effects are the terms proportional 
to the square of the baryon- velocity fluid (v 2 ), giving rise to the quadratic Doppler effect (like those in 
the last line of Eq. (|134p ). or those coupling the photon fluctuations to the baryon velocity (second- line 
from the bottom of Eq. (|134p ). The Vishniac effect corresponds to the terms proportional to Finally 
notice that in the Boltzmann equation (I134p ) the second-order baryon velocity appears. At linear order 
the baryon velocity is irrotational, meaning that it is the gradient of a potential, and thus in Fourier 
space it is parallel to k, and following the same notation of Ref. [55] , we write 

V W(k) = -iv£\k)k. (142) 

The second-order velocity perturbation will contain a transverse (divergence-free) part whose components 
are orthogonal to k = e3, and we can write 

v( 2 )(k) = - ro f(k)e 3+ £ t^S^p, (143) 

m=±l V2 

where form an orhtonormal basis with k. The second-order perturbation is decomposed in a 
similar way, with oj±\ the corresponding components (in this case in the Poisson gauge there is no scalar 
component). Similarly for the tensor perturbation Xij we have indicated its amplitudes as x±2 i n the 
decomposition [57] 

Xij = -V oXm(ei±ze 2 )i(ei±ze 2 ) i . (144) 

m=±2 8 

In computing the source term Si m one has to take into account that in the gravitational part of the Boltz- 
mann equation and in the collsion term there are some terms, like Se v, which still can be decomposed 
in the scalar and transverse parts in Fourier space as in Eq. (|143p . For a generic quantity /(x)v one can 
indicate the corresponding scalar and vortical components with (/v) m and their explicit expression is 
easily found by projecting the Fourier modes of /(x)v along the k = e 3 and (e 2 =F iei) directions 

(/v) m (k) = / |^4 1) (ki)/(k 2 )yr m (k 1 )^. (145) 
Similarly for a term like / (x)Vg(x) one can use the notation 



|^fc l5 (k 1 )/(k 2 )F 1 * m (k 1 )^/ y . (146) 



(27T) 

6.3.4 Integral solution of the second-order Boltzmann equation 

As in linear theory, one can derive an integral solution of the Boltzmann equation (|134p in terms of the 
source term S. Following the standard procedure (see e.g. Ref. [5]) for linear perturbations, we write 
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the left-hand side as A^' + ikfiA^ - t'A® = e~ ikfir,+T d[A^e ik ^i' T }/dri in order to derive the integral 
solution 

A( 2 )(k,n,r ?0 ) = / dr ? 5(k,n,r ? )e a ^- , ' )- T , (147) 
Jo 

(2) 

where % stands by the present time. The expression of the photon moments AL: can be obtained as 
usual from Eq. (|136p . In the previous section we have already found the coefficients for the decomposition 
of source term S 



S(k, n, V )=J2J2 S im (K v)(-iY \l 2lTl Y ^ n ) ■ ( 148 ) 



m=—£ 



In Eq. (|22ip there is an additional angular dependence in the exponential. It is easy to take it into 
account by recalling that 

e ,k.x = J2(if(2e + l)j e (kx)P t (t • x) . (149) 

i 

Thus the angular integral fj 136j) is computed by using the decomposition of the source term (|148j) and 
Eq. ([MD 

AfitMo) = (-l)- m {-ir'iU + l) Tdr/e-^) 

Jo 

x E E (-if 2 s i2m2 ^uAHn-rio)) 

(150) 

where the Wigner 3 — j symbols appear because of the Gaunt integrals 
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Since the second of the Wigner 3-j symbols in Eq. (|150p is nonzero only if m = mj, our solution can be 
rewritten to recover the corresponding expression found for linear anisotropies in Refs. [57] 



21 + 



^ = dne-rW Y, ^wf '"W - r,)] , (152) 
1 Jo % 



where j^ 2 ' m ^ [k(r]Q — 77)] are the so called radial functions. Of course the main information at second order 
is included in the source term containing different effects due to the non-linearity of the perturbations. 
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In the total angular momentum method of Ref. |57j Eq. (|152p is interpreted just as the intergration over 
the radial coordinate (x = t]q — rj) of the projected source term. Another important comment is that, 
as in linear theory, the integral solution (|150p is in fact just a formal solution, since the source term 
S contains itself the second-order photon moments up to I = 2. This means that one has anyway to 
resort to the hierarchy equations for photons, Eq. (|140p . to solve for these moments. Nevertheless, as 
in linear theory, one expects to need just a few moments beyond £ = 2 in the hierarchy equations, and 
once the moments entering in the source function are computed the higher moments are obtained from 
the integral solution. Thus the integral solution should in fact be more advantageous than solving the 
system of coupled equations (|140p . 

6.4 The Boltzmann equation for baryons and cold dark matter 

In this section we will derive the Boltzmann equation for massive particles, which is the case of interest 
for baryons and dark matter. These equations are necessary to find the time evolution of number 
densities and velocities of the baryon fluid which appear in the brightness equation, thus allowing to 
close the system of equations. Let us start from the baryon component. Electrons are tightly coupled to 
protons via Coulomb interactions. This forces the relative energy density contrasts and the velocities to 
a common value, 5 e = 5 P = Sj, and v e = v p = v, so that we can identify electrons and protons collectively 
as "baryonic" matter. 

To derive the Boltzmann equation for baryons let us first focus on the collisionless equation and 
compute therefore dg/dr], where g is the distribution function for a massive species with mass m. One of 
the differences with respect to photons is just that baryons are non-relativistic for the epochs of interest. 
Thus the first step is to generalize the formulae in Section 4 up to Eq. (I80p to the case of a massive 
particle. In this case one enforces the constraint Q 2 = g^vQ^Q" = —m 2 and it also useful to use the 
particle energy E = \J q 2 + m 2 , where q is defined as in Eq. (I62p . Moreover in this case it is very 
convenient to take the distribution function as a function of the variables q l = qn l , the position x 1 and 
time rj, without using the explicit splitting into the magnitude of the momentum q (or the energy E) and 
its direction n % . Thus the total time derivative of the distribution functions reads 



We will not give the details of the calculation since we just need to replicate the same computation we 
did for the photons. For the four- momentum of the particle notice that Q l has the same form as Eq. (|65p . 
while for Q° we find 




(153) 




(154) 



In the following we give the expressions for dx 1 /drj and dq l /drj. 
a) As in Eq. (j69j) dx 1 /drj = Q l /Q° and it turns out to be 




(155) 



b) For dq % jdr\ we need the expression of Q l which is the same as that of Eq. (|65p 




(156) 
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The spatial component of the geodesic equation, up to second order, reads 
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Proceeding as in the massless case we now take the total time derivative of Eq. (|156|) and using Eq. (|157|) 
we find 
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We can now write the total time derivative of the distribution function as 
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(159) 



This equation is completely general since we have just solved for the kinematics of massive particles. As 
far as the collision terms are concerned, for the system of electrons and protons we consider the Coulomb 
scattering processes between the electrons and protons and the Compton scatterings between photons 
and electrons 



dg e , x 
^(X,q^) 



dg P 
drj 



(x,Q,f?) 



{ c ep)QQ>q> + (Ce^pp'q' 
(c e p}qq'Q' , 



(160) 
(161) 



where we have adopted the same formalism of Ref. [12] with p and p' the initial and final momenta of the 
photons, q and q' the corresponding quantities for the electrons and for protons Q and Q'. The integral 
over different momenta is indicated by 



<->, 



d 3 p f d 3 p' 



(2vr) 3 J (2vr) 3 J (2vr) 



d 3 c 



\3 • • • 



(162) 



and thus one can read c e7 as the unintegrated part of Eq. (|85p , and similarly for c ep (with the appropriate 
amplitude |M| 2 ). In Eq. (|160|) Compton scatterings between protons and photons can be safely neglected 



39 



because the amplitude of this process has a much smaller amplitude than Compton scatterings with 
electrons being weighted by the inverse squared mass of the particles. 

At this point for the photons we considered the perturbations around the zero-order Bose-Einstein 
distribution function (which are the unknown quantities). For the electrons (and protons) we can take 
the thermal distribution described by Eq. (JH7J). Moreover we will take the moments of Eqs. (|160p - (|16ip 
in order to find the energy- momentum continuity equations. 



6.4.1 Energy continuity equations 

We now integrate Eq. (|159p over cZ 3 (?/ (27r) 3 . Let us recall that in terms of the distribution function the 
number density n e and the bulk velocity v are given by 

A 



n e = I 7^3^, (163) 



and 

j 1 f d 3 q qn 



where one can set E ~ m e since we are considering non-relativistic particles. We will also make use of 
the following relations when integrating over the solid angle d£l 

dVLrj = [ dtl nVra fc = , / ^ nW = -5 ij . (165) 
J J An 3 

Finally notice that dE/dq = q/E and dg/dq = (q/E)dg/dE. 

Thus the first two integrals just brings n' e and (n e v l )j. Notice that all the terms proportional to 
the second-order vector and tensor perturbations of the metric give a vanishing contribution at second 
order since in this case we can take the zero-order distribution functions which depends only on r\ and E, 
integrate over the direction and use the fact that 5 lJ Xij = 0. The trick to solve the remaining integrals 
is an integration by parts over q l . We have an integral like (the one multiplying — 71)) 

dh ^ ' f A r . = -3n e , (166) 



(2vrp dq i J (2vr) 

after an integration by parts over q l . The remaining integrals can be solved still by integrating by parts 
over q % . The integral proportional to 3>' 1 in Eq. (|159p gives 

AS-E = - Vl n e , (167) 

where we have used the fact that dE/dq 1 = q l /E. For the integral 

d 3 q q l q k dg 
(2vr) 3 E dq 1 



(168) 



the integration by parts brings two pieces, one from the derivation of q l q k and one from the derivation 
of the energy E 

_ 4 [-fS-jL + [ l^ q q - q ---Av k n + ( ^-a^- (169) 
J {2nf 9 E + J (2tt)3 9 EE ~ 4V Ue + J {2nf 9 E^E ' UbJj 
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The last integral in Eq. (|169p can indeed be neglected. To check this one makes use of the explicit 
expression (fHTl) for the distribution function g to derive 

dg qi m e 

g—-—Vig, (170) 



- e 



and 



Thus it is easy to compute 



d 3 a 

gq 1 ^ = 5 l3 n e m e T e + n e m\v l iP . (171) 



(2tt): 



*.* f d3 1 „„2„fc_ , T ,.„.2 T e , „ T e 



,3 



s 3 SsV = — + 3¥ fc u fc ra e — + * fc Vfct/ , (172) 



which is negligible taking into account that T e /m e is of the order of the thermal velocity squared. 

With these results we are now able to compute the left-hand side of the Boltzmann equation (|160p 
integrated over d 3 q/(2ir) s . The same operation must be done for the collision terms on the right hand 
side. For example for the first of the equations in (I160p this brings to the integrals {c ep )QQ' qq > c e~f) pp' qq' ■ 
However looking at Eq. ([86]) one realizes that {c e -y) P pi qq / vanishes because the integrand is antisymmetric 
under the change q •(->• q' and p p'. In fact this is simply a consequence of the fact that the electron 
number is conserved for this process. The same argument holds for the other term {c ep )QQ> qq > . Therefore 
the right-hand side of Eq. (1160p integrated over d 3 q/(2ir) 3 vanishes and we can give the evolution equation 
for n e . Collecting the results of Eq. (|166p to (|172p we find 

7T + + ^ ~ + (*,* " 2 *^) = ■ ( 173 ) 

Similarly, for CDM particles, we find 

+ e* + ^CDM n CDM - 2* fc ) = . (174) 
6.4.2 Momentum continuity equations 

Let us now multiply Eq. (|159p by (q 1 / E) / (2ir) s and integrate over d?q. In this way we will find the 
continuity equation for the momentum of baryons. The first term just gives (n e v 1 )' . The second integral 
is of the type 

o f A a ( e n 6ij + nevi A (175) 



8x3 J (2vr) 3 ^ E E dxi \ m e 
where we have used Eq. (j!7ip and E = m e . The third term proportional to (T~L — is 

d 3 q k d 9Q i , f d 3 q q 2 q i 

An e+ / 7KZ^9-e^-b, (176) 



(2^) 3 ^ 8q k E e J (2Kf y E 2 E 
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where we have integrated by parts over q l . Notice that the last term in Eq. (|176p is negligible being the 
same integral we discussewd above in Eq. (|172p . By the same arguments that lead to neglect the term of 
Eq. (|172p it is easy to check that all the remaining integrals proportional to the gravitational potentials 
are negligible except for 



' J (27rydq k 

The integrals proportional to the second-order vector and tensor perturbations vanish as vector and tensor 
perturbations are traceless and divergence-free. The only one which survives is the term proportional to 
+Uu i in Eq. (fT59]l . 

Therefore for the integral over d 3 qq l /E of the left-hand side of the Boltzmann equation (|159p for a 
massive particle with mass m e (m p ) and distribution function (I87p we find 



d 3 q q l dg e d(n e v % ) 
dr] dr] 

T B V (T>_i-\T/ d , , 



m e J ox3 or] 



(178) 



Now, in order to derive the momentum conservation equation for baryons, we take the first moment 
of both Eq. (1160|) and (j 1611) multiplying them by q and Q respectively and integrating over the momenta. 
Since previously we integrated the left-hand side of these equations over d 3 qq l /E, we just need to multiply 
the previous integrals by m e for the electrons and for m p for the protons. Therefore if we sum the first 
moment of Eqs. (|160p and (|16ip the dominant contribution on the left-hand side will be that of the 
protons 

/ Wf Qi ^ = {Cep{ql + Ql)) QQ'^' + W ■ (179) 

Notice that the integral of the Coulomb collision term c ep (q l + Q l ) over all momenta vanishes simply 
because of momentum conservation (due to the Dirac function 5 4 (q + Q — q' — Q')). As far as the 
Compton scattering is concerned we have that, following Ref. [42], 



{Ce-yQ }pp'qq' — {Ce"/P }pp'qq' > (1^0) 

still because of the total momentum conservation. Therefore what we can compute now is the integral 
over all momenta of c e7 p\ Notice however that this is equivalent just to multiply the Compton collision 
term C(p) of Eq. (|86j) by p L and integrate over d 3 p/ (2tr 3 ) 

/d 3 v 
T^-^C(p) . (181) 

where C(p) has been already computed in Eqs. (I122p and (|123p . 

We will do the integral (|18ip in the following. First let us introduce the definition of the velocity of 
photons in terms of the distribution function 
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where p-y = p-y/ 3 is the photon pressure and p 7 the energy density. At first order we get 



l v 0»= /^AWrf 

4-7T 



(183) 



where A is the photon distribution anisotropies defined in Eq. (|132p . At second order we instead find 



4 u 



(2)i 



3 2 



2 J 4vr n 3 7 7 " 



(184) 



Therefore the terms in Eqs. (|122|) and (|123|) proportional to /^(p) and /^(p) will give rise to terms 
containing the velocity of the photons. On the other hand the terms proportional to /g (p) and /og (p) ; 
once integrated, vanish because of the integral over the momentum direction n l , J d£ln l = 0. Also the 
integrals involving P2( v • n ) = [3(v • n) 2 — l]/2 in the first line of Eq. (|122p and (|123|) vanish since 



J dttP 2 (v ■ n) n l = v k v j J dUn k rijn l = . 



(185) 



where we are using the relations (|165j) . Similarly all the terms proportional to v, (v-n) 2 and v 2 do not give 
any contribution to Eq. (I18ip and, in the second-order collision term, one can check that / d£lY2(n)n l = 0. 

Then there are terms proportional to (v • n)f(°\p), (v • n)pdf^ /dp and (v • n)pdfQ /dp for which we 
can use the rules (|133j) when integrating over p while the integration over the momentum direction is 



dQ 
4tt 



(v • n)n l = v k 



dQ 

4tt ' 



n k n l 



Finally from the second line of Eq. Q123)) we get three integrals. One is 



n 



where p 7 is the background energy density of the photons. The second comes from 



d 3 p • 
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dp 
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(186) 



(187) 



(188) 



47T 



-r^n 3 



where we have used the rules ([IMP , Eq. (fT65]) and f(&l/4n) n i nin k n l = (S^S kl + 5 ik 5 l i + S il S^ k )/15. In 
fact the third integral 



d 3 P i f w(i)/ \ 
P (v-n)/^ ; (p) 



(27T) 



(189) 



exactly cancels the previous one. Summing the various integrals we find 



dp 

(2^p 

4 v( 2 ) 
3 2 



C(p)p = n e a T p 1 



|(v«-v«) 



d^A^ 2 ) 
4vr 2 



-n 



+ ^ 1) (v (1) -v 7 1 )) + J ^AW(v.n)n + Af)v 



(190) 
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Eq. (|190p can be further simplified. Recalling that S-y' = Aq 1 ^ we use Eq. (|20ip and notice that 

J ^A« (v • njn* = t^D* + 1<M« , (191) 
where the photon quadrupole IB? is defined as 

A(2)\ 



(192) 



Thus, our final expression for the integrated collision term (|181|) reads 

d 3 p 



(2tt) 3 
4 
3 



C(p)p l = n e a T p 1 



3 l 7 J 3 I 2 2 



+ l(5W+AS 1 J)(t;W i -t;^) + t;«n? 



(193) 



We are now able to give the momentum continuity equation for baryons by combining m p dg p /dr] from 
Eq. (fTT8l) with the collision term (fj"8TD 



d{p b v % 
drj 



+ 4(n - V'fav* + $*e* + V + e*+* f p b ^ 



+ e 



4 /V 2 )* ^ 2)i ' 



3^2 2 

+ | + A« + 3«) (««< - «W«) + ««Di' 



(194) 



where p b is the baryon energy density and, as we previously explained, we took into account that to a 
good approximation the electrons do not contribute to the mass of baryons. In the following we will 
expand explicitly at first and second-order Eq. (|194D . 



6.4.3 Second-order momentum continuity equation for baryons 

At first order we find 



3 Oh V 



L>Pi( v m*- v ( 



i)< 



(195) 



drj 3' p b 

At second order there are various simplifications. In particular notice that the term on the right-hand 
side of Eq. (|194p which is proportional to 6 b vanishes when matched to expansion of the left-hand side 
by virtue of the first-order equation (|195j) . Thus, at the end we find a very simple equation 
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with r' = —n e (JTa. 



6.4.4 Second-order momentum continuity equation for CDM 

Since CDM particles are collisionless, at first order we find 



ch] 

At second order we find 



^ M +^cdm + ^ (1) ' 4 = 0. (197) 
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7 CMB anisotropies at second-order at all scales: analyt- 
ical approach 

7.1 Towards a second-order CMB radiation transfer function 

As pointed out in Sec. [6] various non-primordial sources of non-Gaussianity for the CMB anisotropies can 
arise from the non-linear evolution of the cosmological perturbations, including the Sunyaev-Z'eldovich, 
ISW (Rees-Sciama) effects, and the gravitational lensing, with possible correlations between these con- 
tributions. Here we will focus on another relevant source of non-Gaussianity: the non-linear effects 
operating at the recombination epoch. The dynamics at recombination is quite involved because all the 
non-linearities in the evolution of the baryon-photon fluid at recombination and the ones coming from 
general relativity should be accounted for. 

The following section can be considered as an application of the Boltzmann equations found previ- 
ously. Despite they depend on a high numbers of terms, because of the apperance within second-order 
perturbation theory of products of fist-oder perturbations, it is remarkable that an analytical study of 
the complicated recombination dynamics is possible. This allows to account for those effects that at the 
last scattering surface produce a non-Gaussian contribution to the CMB anisotropies that add to the 
primordial one. Such a contribution is so relevant because it represents a major part of the second-order 
radiation transfer function which must be determined in order to have a complete control of both the 
primordial and non-primordial part of NG in the CMB anisotropies and to gain from the theoretical side 
the same level of precision that could be reached experimentally in the near future [5"! HI [58]. 

In order to achieve this goal, we consider the Boltzmann equations derived in Sec. [6] at second-order 
describing the evolution of the photon, baryon and CDM fluids, and we manipulate them further under the 
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assumption of tight coupling beteween photons and baryons. This leads to the generalization at second- 
order of the equations for the photon energy density and velocity perturbations which govern the acoustic 
oscillations of the photon-baryon fluid for modes that are inside the horizon at recombination. The 
evolution is that of a damped harmonic oscillator, with a source term which is given by the gravitational 
potentials generated by the different species. An interesting result is that, unlike the linear case, at 
second-order the quadrupole moment of the photons is not suppressed in the tight coupling limit and it 
must be taken into account. 

The analytical solutions for the acoustic oscillations of the photon-baryon fluid at second-order are 
derived adopting some simplifications which are also standard for an analytical treatment of the linear 
CMB anisotropics, and which nonetheless allow to catch most of the physics at recombination. One 
of these simplifications is to study separately two limiting regimes: intermediate scales which enter the 
horizon in between the equality epoch (r] e q) and the recombination epoch (j] r ), with r/" 1 <C k <C T]~ q , 
and shortwave perturbations, with k 3> , which enter the horizon before the equality epoch. Here our 
main concern is to provide a simple estimate of the quantitative behaviour of the non-linear evolution 
taking place at recombination, offering at the same time all the tools for a more accurate computation. 
We find that the second-order CMB anisotropies generated on the last scattering surface do not reduce 
only to the energy density and velocity perturbations of the photons evaluated at recombination, but a 
number of second-order corrections at last scattering arise from the Boltzmann equations of Ref. [38] in 
the form of products of first-oder perturbations. 

We will see that the dynamics at recombination is indeed dominated on small scales by the non- 
linear evolution of the second-order gravitational potentials feeded by the Cold Dark Matter density 
perturbations. The gravitational potentials determine the energy density fluctuations of the photons at 
recombination and their effects show up in the CMB anisotropies as 

^r = \*®+*®, (199) 

(2) 

which is the usual term due to the intrinsic fractional temperature fluctuation Aq /4 on the last scattering 

(2) 

surface (Ag is the monople of the photon distribution function) and the gravitational redshift due 
to the gravitational potential. However the analysis of the remaining contributions that come in the 
form of products of first-oder perturbations is equally important. The reason is that one of the central 
quantities we are interested in is the contamination to the primordial non-Gaussianity that is produced 
by the secondary effects. In that respect the reasonable question to ask is which kind of primordial 
non-Gaussianity a given secondary effect can contaminate most, wether it is of the so called "local type" 
or of the "equilateral type" , for example. Therefore it might well be the case that, even if some secondary 
effects appear to be the dominant ones, they might give a high contamination to a given type of primordial 
non-Gaussianity, but a low contamination to a different kind of primordial non-Gaussianity, for which 
"subdominant " terms, on the other hand, represent the strongest contaminant. This section and the 
following two provide a clear example: we anticipate here that while the term in Eq. (|199p will mainly 
mimic an equilateral type of primordial non-Gaussianity, second-order effects that come as products of 
first-order times first-order perturbations actually mainly contaminate local primordial non-Gaussianity. 

Notice that the case k 3> rj~g has been treated in two steps. First we just assume a radiation dom- 
inated universe, and then we provide a much more complete analysis by solving the evolution of the 
perturbations from the equality epoch onwards taking into account that the dark matter perturbations 
around the equality epoch tend to dominate the second-order gravitational potentials. As a byproduct, 
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this last step provides the Meszaros effect at second-order. In deriving the analytical solutions we have ac- 
curately accounted for the initial conditions set on superhorizon scales by the primordial non-Gaussianity. 
In fact the primordial contribution is always transferred linearly, while the real new contribution to the 
radiation transfer function is given by all the additional terms provided in the source functions of the 
equations. Let us stress here that the analysis of the CMB bispectrum performed so far, as for example 
in Ref. [59], adopt just the linear radiation transfer function (unless the bispectrum originated by specific 
secondary effects, such as Rees-Sciama or Sunyaev-Zel'dovich effects, are considered, see on this Sec. [8]). 
The formalism to provide in a systematic way the second-order CMB radiation transfer function will be 
reviewd in Sec. [TO]. 



7.2 The Boltzmann equations in the tightly coupled limit 

We now derive the moments of the Boltzmann equations for photons in the limit when the photons 
are tightly coupled to the the baryons (the electron-proton system) due to Compton scattering. This 
leads to the governing equations for the acoustic oscillations of the photon-baryon fluid. The well-known 
computation at linear order is briefly reviewed in IC.ll under some standard simplifying assumptions. 
Here we will focus on the derivation of the equations at second-order in the perturbations, pointing out 
some interesting differences with respect to the linear case. In particular note that , while we already 
know that the L.H.S. of the equations at second-order will have the same form as for the linear case, 
the source term on the R.H.S. of the moments of the Boltzmann equations will also consist of first-order 
squared terms. 



7.2.1 Energy continuity equation 

Let us start by integrating Eq. (|134p over dCl n /4-K to get the evolution equation for the second-order 
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photon energy density perturbations Aq 
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where we have used the explicit definition for the second-order velocity of the photons [38] 
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We can now make use of the tight coupling expansion to simplify Eq. (|200p . In the L.H.S. we use 
d iV ^ ]i = dM 1 ^ = 3*W - <5p = 31» W - 3AS74 obtained in the tightly coupled limit from Eqs. ([C?fjl 
and (ICT9D . On the other hand in the R.H.S. of Eq. (f200D 
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by using Eq. (|C.4p and the evolution equation for the photon veleocity (|C.5p . Here we introduce the 
baryon-photon ratio 



4p 7 



We thus arrive at the following equation 



(203) 
(204) 



where the source term is given by 
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Notice that the integral in the second line of Eq. (I200p has been computed by expanding the linear 
anisotropies as = £<(2M-1)A? } P^v-n) and expressing explicitly the dependence on as Pi(v-n) = 
^th, P 2 (v • n) = (3(v • n) 2 - l)/2 and so on. This allows to perform the derivative dA^ /dri 1 . It turns 
out that the term proportional to the first-order quadrupole vanishes, and higher order terms can be 
neglected because of the tight coupling approximation. Therefore the only term that is non negligible is 
the dipole term and using A^ = 4i>/3 in the tight coupling limit, one obatins that the integral is equal 
to (16/3)($( 1} + ^M)' i v l . 



7.2.2 Velocity continuity equation 



We now derive the second moment of the Boltzmann equation (|134|) and then we take its tight coupling 
limit. The integration of Eq. (11341) over d£l n n l /47r yields the continuity equation for the photon velocity 
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The difference between the second-order baryon and photon velocities (v^ 1 — u-y ) appearing in Eq. (|206|) 
is obtained from the baryon continuity equation which can be written as (see Ref. [38]) 

P 



,(2)< 



18 



- 2(A[ ) 1 ) +$( 1 ))(^ 1 ) l -^ 1 ) i ). (207) 

We want now to reduce Eq. (|206p in the tightly coupled limit. We first insert the expression (|207p in 
Eq. (|206p . Notice that the last three terms in Eq. (|207p will cancel out. On the other hand in the tight 
coupling limit expansion one can set v^ 1 = and v^ 1 = in the remaining terms on the R.H.S. 
of Eq. (12071 Thus Eq. ([206]) becomes 

(v®* + 2jy + n-^—{v^ + 2J) + 1^ + = 
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1 + R 1 + R 

where in the tightly coupled limit we are neglecting the first-order quadrupole and (higher-order moments) 
of the photon distribution since it is suppressed by 1/r with respect to the other terms. Moreover the 
integral in the fourth line of Eq. (|206p has been computed following the same steps used for the analogous 
integral in Eq. (|200p . In this case it turns out that the contribution from the first-order dipole term 
vanishes, while the one from thequadrupole is non zero, but this and higher-order terms can be neglected 
because of the tight coupling approximation, so that the integral is in fact negligible. Next for the term 
like t'5^\v^ 1 — v^p*) we employ the relation previously derived in Eq. (|202p with 6fP = 3Aqq^/4 and 
we use the first order tight coupling equations (|C.ip and (|C5[) in order to further simplify Eq. (|208|) . We 
finally obtain 

v ? 1 ' + n T^ v ? 1 + 7^75 + * (2),< = S v . ( 209 ) 
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We have spent some time in giving the details of the computation for the photon Boltzmann equations 
at second-order in the perturbations. As a summary of the results obtained so far we refer the reader 
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to Eqs. (|204[) and (|209p as our master equations which we will solve in the next sections. In particular 
Eq. (|209p is the second-order counterpart of Eq. (IC.5|) for the photon velocity in the tight coupling regime. 
Notice that there are two important differences with respect to the linear case. One is that, in Eq. (|209p . 
there will be a contribution not only from scalar perturbations but also from vector modes which, at 
second-order, are inevitably generated as non-linear combinations of first-order scalar perturbations. In 
particular we have included the vector metric perturbations u/ in the source term. Second, and most 
important, we have also kept in the source term the second-order quadrupole of the photon distribution 
niy 2 ^' 7 . At linear order we can neglect it together with higher order moments of the photons since they 
turn out to be suppressed with respect to the first two moments in the tight coupling limit by increasing 
powers of 1/r. However in the next section we will show that at second order this does not hold anymore, 
as the photon quadrupole is no longer suppressed. 

Finally following the same steps leading to Eq. (IC.6j) at linear order we can derive a similar equation 

(2) 

for the second-order photon energy density perturbation Aq which now will be characterized by the 
source terms 5a and S v 

(A$" - WW) + H-^ (A£> - 4*( 2 )') - c 2 V 2 (AS - 4*( 2 )) 
™ \ R 4 




where we have introduced the photon-baryon fluid sound speed c s = 1/^/3(1 + R). 

7.2.3 Second-order quadrupole moment of the photons in the tight coupling limit 

Let us now consider the quadrupole moment of the photon distribution defined in Eq. (|192p and show 
that at second-order it cannot be neglected in the tightly coupled limit, unlike for the linear case. We 
first integrate the R.H.S. of Eq. ()134j) over dVt^r^n 3 — 5*- ? /3)/47r and then we set it to be vanishing in 
the limit of tight coupling. 

The integration involves various pieces to compute. For clarity we will consider each of them sepa- 
rately. The term A 00 does not contribute. For the third term we can write, from Eq. f|135[) 
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so that the integral just brings IL^^/IO, since the only contribution in Eq. (|212p comes from A( 2 )/10 
with all the other terms vanishing. The following nontrivial integral is 
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where the baryon velocity appearing in i"2( v • n) is first order and we make use of the relations f|165|) 
together with 



^ n i n j n k n l = — (5 ij 5 kl + 5 ik 5 lj + 5 il 5 jk ) . (214) 
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The integrals of 5^ Aq 1 ^ , 6e (v • n) and v^ 2 ) • n vanish and 



uA^ y ^^2(v • n) (nW - l#A = ^vA? (W - V' 



— ( cV - -<fV 
15 V 3 



(215) 



where in the last step we take A^ = 4^/3 in the tight coupling limit. Similarly the integral of 14(v • n) 2 
brings 

Uv k v e J ^n k n e (nW - ^6 ij \ = ^ (W - ^S ij v 2 ^j . (216) 

The integral of 2(v • n)A^ 1 ) can be performed by expanding the linear anisotropies as A^ 1 ) = ^2i(2£ + 
lJA^P^v • n). We thus find 
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where we have used Eq. (|165p and 0£>2 indicates all the integrals coming from the multipoles £ > 2 in 
the expansion (for £ = and £ = 2 they vanish.) In fact we have dropped the 0{£ > 2) since they are 
proportional to firts-order photon moments £ > 2 which turn out to be suprressed in the tight coupling 
limit. Finally the term proportional to (v • n)A 2 (1 — -P2( v • n )/5 gives a vanishing contribution. 

Collecting all the various pieces we find that the third moment of the R.H.S. of Eq. fj 134j) is given by 
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(218) 

Therefore in the limit of tight coupling, when the interaction rate is very high, the second-order quadrupole 
moment turns out to be 

IL™ ~ | (vV - , (219) 

by setting Eq. ()218p to be vanishing (the term multiplying 6^ goes to zero in the tight coupling limit 
since it just comes from the first-order collision term). At linear order one would simply get the term 

Qr'lI^^/lO implying that, in the limit of a high scattering rate r', 11^ goes to zero. However at 
second-order the quadrupole is not suppressed in the tight coupling limit becasue it turns out to be 
sourced by the linear velocity squared. 
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7.3 Second-order CMB anisotropies generated at recombination 

The previous equations allow us to follow the evolution of the monopole and dipole of CMB photons at 
recombination. As at linear order, they will appear in the expression for the CMB anisotropies today 
A( 2 )(k, n, t/o ) together with various integrated effects. Our focus now will be to obtain an expression for 
the second-order CMB anisotropies today A^ 2 ^(k, n, r/o) from which we can extract all those contributions 
generated specifically at recombination due to the non-linear dynamics of the photon-baryon fluid. This 

(2) 

expression will not only relate the moments A)' today to the second-order monopole and dipole at 
recombination as it happens at linear order, but one has to properly account also for additional first- 
order squared contributions. Let us see how to achieve this goal in some details. 

As we have seen previously, it is possible to write down an integral solution of the photon Boltzmann 
equation (1134p in Fourier space. We wrote 
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in order to derive a solution of the form 
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Here fi = cosi? = k • n is the polar angle of the photon momentum in a coordinate system such that 
e3 = k. At second-order the source term has been computed in Ref. [38] and can be read off Eq. (|134p 
and Eq. ([rgB]) to be 
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(222) 



The key point here is to isolate all those terms that multiply the differential optical depth r'. The reason 
is that in this case in the integral (|22ip one recognizes the visibility function g(rj) = — e~ T r' which is 
sharply peaked at the time of recombination and whose integral over time is normalized to unity. Thus 
for these terms the integral just reduces to the remaining integrand (apart from the visibility function) 
evaluated at recombination. The standard example that one encounters also at linear order is given by 
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the first term appearing in the source S, Eq. (|222p . that is — r Aqo . The contribution of this term to the 
integral (12210 just reduces to 

A( 2 )(k, n, no) = r d7]e ik ^-vo) e -T(_ T ^ A m ^ e ikn{ v *-r, ) A^(^) , 

JO 

(223) 

where 77* is the epoch of recombination and, in the multipole decomposition (|136p . Eq. (|223p brings the 
standard result 

Aglfo) oc Aw MMHv* - m)) , (224) 

having used the Legendre expansion e = ^2g{i) (2i + l)j^(A;x)P^(k • x). In Eq. (|224p the monopole 
at recombination is found by solving the Boltzmann equations Eqs. (|204p - (l209p derived in tight coupling 
limit. 

Looking at Eq. (j222H we recognize immediately some terms which multiply explicitly t' (the first one 
discussed in the example above and the last two lines of Eq. (|222p ). However it is easy to realize from 

the standard procedure adopted at the linear-order that such terms are not the only ones. This is clear 

(2) 

by focusing, as an example, on the term — 4n*<I> \ in the source S which appears in the same form also 
at linear order. In Fourier space one can replace the angle /i with a time derivative and thus this term 
gives rise to 

- 4ik r drie^-Vo)^ = _ 4 P° d $(2) e ~r d f jkrfri-vo)) 
Jo Jo W] V / 

= =A£°dr l e ik ^-^ ) e- T (^^'-T'^), (225) 

where, in the last step, we have integrated by parts. In Eq. (|225p the time derivative of the gravitational 
potential contributes to the Integrated Sachs- Wolfe effect, but also also a t' results implying that we 
have also to evaluate $^ 2 ^ at recombination. Thus, in the following, we look for those terms in the 
source (|222|) which give rise to a r 1 factor in the same way as for — 4n l< 5 \ . In particular let us consider 
the combination in Eq. (I222p 

C = 8AW(*W-n i $J ) )-2n i ($ (1) +* (1) )(AW+4$W) ii 

= SAW^W' - Sn'CAW^Wj.i + A^rt&f - tot&WXi , (226) 

where for simplicity we are setting ~ We already recognize terms of the form n l di(-). Moreover 
we can use the Boltzmann equation (|129p to replace ntA p in Eq. (|226p . This brings 
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In fact we will not be interested for our purposes in the first three terms of Eq. (|227p . since they will not 
contribute to the anisotropies generated at recombination. 
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Therefore, as a result of Eqs. (|222p . (|225p and (|227p . we can rewrite the source term (|222p as 
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In Eq. (|228p .S* contains the contribution to the second-order CMB anisotropics created on the last 
scattering surface at recombination, while S' includes all those effects which are integrated in time from 
the last scattering surface up to now, including the second-order Integrated Sachs- Wolfe effect and the 
second-order lensing effect. Since the main concern of this Section is the CMB anisotropies generated at 
last scattering, from now on we will focus only on the contribution from the last scattering surface S*. 
In particular notice that, following the same steps that lead to Eq. (I223p . the first two terms of Eq. (|229p 
give rise to the CMB temperature anisotropies as written in Eq. (|199p . 



7.4 Tightly coupled solutions for the second-order perturbations 



In this section we will solve the tightlty coupled limit of the Boltzmann equations (j204p and ()206p at 
second-order in perturbation theory. We will proceed as for the linear case, focusing on the two limiting 
cases of perturbation modes entering the horizon respectively much before and much after the time of 
equality. The solution of Eq. (|21ip can be written as 
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where for simplicity we use the shorthand notation sin^, 77'] = sin [A; (r s (77) — r s (rj'))], and r s (rj) is the 
sound horizon 

T,{rj) = fdn'csirf). (232) 
Jo 

The source terms are given in Eq. ()205p and (|210p . Notice that we can write S' A + j^Sa = (5a (1 + 
R))'/l + R so that we can perform an integration by parts in Eq. (|231h leading to 

[1 + R(ji)] 1/4 (A$ - 4^ (2) ) = Acos[kr s (r])\ + B sm[kr s (7/)] 
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In order to give an analytical solution that catches most of the physics underlying Eq. (|23ip and which 
remains at the same time very simple to treat, we will make some simplifications which are also usually 
adopted in linear theory, see e.g. Ref. |87j . In particular we will treat R = R* as a constant evaluated at 
the time of recombination. In this way, with the presence of R in the varying cosines and sines, we can 
keep track of a damping of the photon velocity amplitude with respect to the case R = which prevents 
the acoustic peaks in the power-spectrum to disappear. Treating R as a constant is justified by the fact 
that for modes within the horizon the time scale of the ocillations is much shorter than the time scale on 
which R varies. If R is a constant the sound speed is just a constant 

(234) 



^3(1 + R*) 

and the sound horizon is simply r s {rf) = c s rj. 

There is also another good reason not to forget about R, even tough R < 1, which is particularly 
relevant when estimating the second-order contributions to the CMB anisotropics in (|199p . As we will 
see in Sec. [9l that combination on small scales comes actually as proportional to R, so that neglecting R 
would lead to miss one of the most significant CMB contribution at second order. 

Another simplification consists in solving the evolutions of the perturbations in two well distinguished 
limiting regimes. One regime is for those perturbations which enter the Hubble radius when matter is 
the dominant component, that is at times much bigger than the equality epoch, with k <C k eq ~ r\~ q , 
where k eq is the wavenumber of the Hubble radius at the equality epoch. The other regime is for those 
perturbations with much smaller wavelenghts which enter the Hubble radius when the universe is still 
radiation dominated, that is perturbations with wavenumbers k S> k eq ~ f]eq- ^ n ^ ac ^ we are interested 
in perturbation modes which are within the horizon by the time of recombination 77*. Therefore we will 
further suppose that 77* 3> ij eq in order to study such modes in the first regime. Even though 77* S> r] eq is 
not the real case, it allows to obtain some analytical expressions. 
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7.4.1 Setting the initial conditions: primordial non-Gaussianity 

The integration constants A and B are fixed according to the initial conditons for the second-order 
cosmological perturbations. These refer to the values of the perturbations on superhorizon scales deep 
in the radiation dominated period. We will consider the case of initial adiabatic perturbations, for which 
there exist some useful conserved quantities on large scales which as such carry directly the information 
about the initial conditons. 

As explained in Sec. H] such a conserved quantity is given by the curvature perturbation Q and its 
conserved value allows to set the initial conditions for the metric and matter perturbations accounting 
for the primordial contributions. At linear order during the radiation-dominated epoch and on large 
scales = — 2\pW /3. On the other hand, after some calculations, one can easily compute A<^ 2 ) for a 
radiation dominated epoch 

AC^ = 7 -(^)\ (235) 

where in Eq. (|38[) one uses that on large scales S^p^/p-y = —2^^ and the energy continuity equation 
<5 (1) 'p 7 + mS^Pi - 4^«V 7 = 0. Therefore we find 

C (2) = _tf(2) + ^L + ^(l)2 (0)j (236) 

where we are evaluating the quantities in the large scale limit for rj — > 0. Using the parametrization (|39p 
at the initial times the quantity Aqq — 4^( 2 ) is given by 

Ag - 4^ 2 > = 2(9a NL - 7)* (1)2 (0) . (237) 

Since for adiabatic perturbations such a quantity is conserved on superhorizon scales, it follows that the 
constant B = and A = 2(9a NL - 7)^( 1 ) 2 (0). 

Eqs. (I231D and (|233l) are analytical expressions describing the acoustic oscillations of the photon- 
baryon fluid induced at second-order for perturbation modes within the horizon at recombination. In the 
following we will adopt similar simplifications already used for the linear case in order to provide some 
analytical solutions. In particular, if in Eq. (|233|) we treat R as a constant we can write, using the initial 
conditions determined above, 

(Ag - 4tt< 2 >) = 2(9a NL - 7)*« 2 (0) cos[kr s (r))} 

1*7) 

- 4kc s J drf ($ (2) (v') + ^ (2) (??') + R<Z> {2) ) sin*fa,»/] 

+ / di]'S A (v') cos fe [??,?/] - ^Sa(O) sm[kr s (rj)] 
Jo k 

~ \t~ r^Svtf) rinkhr/]. (238) 
3 kc s Jo 

7.5 Perturbation modes with k k eq 

In order to study the contribution to the second-order CMB anisotropies coming from perturbation 
modes that enter the horizon during the radiation dominated epoch, we will assume that the second-order 



56 



gravitational potentials are the ones of a pure radiation dominated universe throughout the evolution. 
Though not strictly correct, this approximation will give us the basic picture of the acoustic oscillations 
for the baryon-photon fluid occuring for these modes. Also for the second-order case, in Section [7.71 we 
will provide the appropriate corrections accounting for the transition from radiation to matter domination 
which is indeed (almost) achieved by the recombination epoch. Before moving into the details a note of 
caution is in order here. At second order in the perturbations all the relevant quantities are expressed as 
convolutions of linear perturbations, bringing to a mode-mode mixing. In some cases in our treatment for 
a given regime under analysis (k 3> k eq or k <C k eq ) we use for the first-order perturbations the solutions 
corresponding to that particular regime, while the mode-mode mixing would require to consider in the 
convolutions (where one is integrating over all the wavenumbers) a more general expression for the 
first-order perturbations (which analytically does not exist anyway). For the computation of the CMB 
bispectrum this would be equivalent to consider just some specific scales, i.e. all the three scales involved 
in the bispectrum should correspond approximately to wavenumbers k 3> k eq or k <C k eq , and not a 
combination of the two regimes (a step towards the evaluation of the three-point correlation function 
has been taken on Ref. [37] where it was computed in the in so-called squeezed triangle limit, when one 
mode has a wavelength much larger than the other two and is outside the horizon). 

In Appendix C we show how to solve the Boltzmann equations at the linear level in the various 
regimes. In fact for k S> k eq the acoustic oscillations can be solved in an alternative way (see Ref. [39]) 
for this procedure at linear order). One can start directly from Eq. (|204p where we can neglect the 
gravitational potential term ^>( 2 ) . The reason is that, as it happens at linear order, the second-order 
gravitational potentials decay at late times as f]~ 2 , while the second-order velocity iffl* oscillates in time. 
Let us now see that in some details. 

The evolution equation for the gravitational potential is given by Eq. (|B.13|) and is characterized 
by the source term Sy, Eq. (|B.14|) . In particular the source term contains the second-order quadrupole 

moment of the photons liffl^ . We saw in Section 17.2.31 that at second-order the quadrupole moment 
is not suppressed in the tight coupling limit, being fed by the non-linear combination of the first-order 
velocities, Eq. (1219j) . For the pertubation modes we are considering here the velocity at late times is 
oscillating being given in Appendix C in Fourier space. Since the linear gravitational potential decays in 
time and for a radiation dominated period T~L = 1/r], it is easy to check that the dominant contribution 
at late times to the source term iS 7 simply reduces to 

S, - ^ 2 ^nf * i ee F(kl ^ 2 ' k) C (0)<> (0) srn(kic sV ) sinfec^) , (239) 



where 



C =-AA (240) 



and the sound speed is c s = 1/ \/3(l + R). Before proceeding further let us explain the notation that 
we are using. The equivalence symbol will be used to indicate that we are evaluating the expression in 
Fourier space. At second-order in perturbation theory most of the Fourier transforms reduce to some 
convolutions. We will not indicate these convolutions explicitly but just through their kernel. For example 
in Eq. (1239P by i ? (ki,k2,k) we actually indicate the convolution operator 

F = ^_ y^ 3 A; 1 d 3 M (1) (ki+k2-k)F(k 1 ,k2,k). (241) 
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In the specific case of Eq. (|239|) the kernel is given by 

F(k 1 ,k 2 ,k) = (k - k f 2 k - k2) -lk 1 -k 2 . (242) 

The choice of these conventions is due not only for simplicity and to keep our expressions shorter, but 
also because at the end we will be interested to the bispectrum of the CMB anisotropies generated 
at recombination, and the relevant expressions entering in the bispectrum are just the kernels of the 
convolution integrals. 

Having determined the leading contribution to the source term at late times, we can now solve the 
evolution equation (|B.13|) . Since the source term scales like r/~ 2 it is useful to introduce the rescaled 
variable x = ri 2 ^ 2 \ Eq. (IB.13j) then reads 



X" +[k 2 c 2 s --)x = il 2 S~ i . (243) 



2 
rf 

For perturbation modes which are subhorizon with krj 3> 1 the solution of the homogeneous equation is 
given by 

Xhom. = Acos(kc s rj) + B sm(c s ki]) , (244) 
from which we can build the general solution 

X = X hom . + X + [ drj'^S^') - x- jf drf^&yirf) , (245) 

where W = —kc s is the Wronskian, and x+ = cos(kc s r]), X- = sin(c s A;?7). Using Eq. (|239p the integrals 
involve products of sines and cosines which can be performed giving 

X = Xhonu-^ff (0)0) 

k [2kik 2 cos(kic s i])cos(k2C s r]) - 2kik 2 cos(kc s r]) + (k\ + k 2 - k 2 ) sm(kic s T]) sm(k 2 c s r])] 
X t-4 i 1.4 i 1.4 _ ipp _ 91,22,2 _ ?k 2 k 2 ' ' ' 

Thus the gravitational potential at late times is given by 
= _ 3 *«(0)^-^<>(0)<>(0) 



[2/eiA:2Cos(A;ic s 77)cos(A:2C s 77) — 2kik2Cos(kc s r]) + (k 2 + k\ — k 2 ) sin(kic s r)) sin(k2C s r])] 
X i-4 _i_ J.4 , p oi,2i,2 91,21,2 91,2 1,2 ' 

(247) 

where we have set the integration constant B = and j4 = — 3^( 2 \0)/(kc s ) 2 in order to match the 
homogenoeus solution at late times. Here ^( 2 )(0) is the intial condition for the gravitational potential 
taken on large scales deep in the radiation dominated era which will be determined in Section 17.5.21 

Eq. (|247p shows the result that we anticipated: also at second order the gravitational potential varies 
in time oscillating with an amplitude that decays as r\~ 2 . Let us then take the divergence of the {i — 0) 



58 



Einstein equation (|3ip expanded at second-order 



1 u 



+ ^00 w 7 



(248) 



which, using the first-order (i — 0) Einstein equation and ~ reduces to 



1^(2)' , ^$(2) _ ^(l)'^(l) 

2 2 



2 V 1 + ^00^7 



(249) 



Since during a radiation dominated period is given by Eq. (IB. 11}) and at late times it decays 
oscillating, it is easy to see that (\J/( 1 )'9 i ^/( 1 )) will be oscillating and decaying as r]~ 4 and thus can be 
neglected with respect to ^v 2 )', which oscillates with an amplitude decaying as i]~ 2 . Also T-L^ 2 ^ turns out 
to be subdominant. Recall that <1?( 2 ) = — and is dominated by the second-order quadrupole 
of the photons in Eq. ()B.14p . so that $( 2 ) scales like \E r( - 2 ) but there is the additional damping factor of 
the Hubble rate % = 1/rj. Thus the dominant terms give 



(250) 



Eq. (|250p allows to proceed further in a similar way as for the linear case by using the results found so far, 
Eqs. (|247p and (|250p . in the energy continuity equation (|204p . In Eq. (|204p the first- and second-order 

gravitational potentials can be neglected with respect to the remaining terms given by A ^ and which 
oscillate in time. Thus, replacing the divergence of the second-order velocity by the expression (|250p . 
Eq. (|204p becomes 



.(2)/ 

*00 



L)2\' 



4 00 



which, using the first-order equation (jC.ip . further simplifies to 



A 



(2)/ 
00 



m 2 



: V 2 ^(2)' 



(251) 



(252) 



where we have kept only the dominant terms at late times. 

The gravitational potential is given in Eq. (|247p . so the integration of Eq. (|252p gives 



A 



6¥ 2 \o)cos(kc s rj) + 2^< ) (0)^(0)^ 



(2) 



[2kik2Cos(kic s rj)cos(k2C s r]) — 2kik2Cos(kc s r]) + (kf + k\ — k 2 ) sin(A; 1 c s r/) sin(/c 2 c s ?7)] 

X 1.4 _i_ 1.4 _i_ 1.4 1,2 u2 _ OjU2jU2 _ 9Z~2l2 ' 

(253) 

Needless to say, modes for k 3> ko, where k~j^ indicates the usual damping length, are supposed to be 
multiplied by an exponential e 



-{k/k D f 



(see, e.g. @2]). 
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7.5.1 Vector perturbations 



So far we have discussed only scalar perturbations. However at second-order in perturbation theory an 
unavoidable prediction is that also vector (and tensor) perturbation modes are produced dynamically 
as non-linear combination of first-order scalar perturbations. In particular notice that the second-order 
velocity appearing in Eq. (]229p . giving rise to a second-order Doppler effect at last scattering, will contain 
a scalar and a vector (divergence free) part. Eq. (|250p provides the scalar component of the second-order 
velocity. We now derive an expression for the velocity that includes also the vector contribution. 

The (second-order) vector metric perturbation uo l when radiation dominates can be obtained from 
Eq. (jETTo) 



,(2)j 



+ A 



00 ^7 



(254) 



where we have dropped the gravitational potentials *$>^ ~ $W which are subdominant at late times. 
On the other hand from the velocity continuity equation (|209p we get 

3 



(2)i/ 1 A (2),< 
7 + 4 n 



^00 



iKo' 2 )' 



+ 



-vWdjvWi - 2oj 1 - -u k 



fk^ ki 



(255) 



neglecting the term proportional to R and the decaying gravitational potentials. Using the tight coupling 
equations at first-order and integrating over time one finds 



-2cj 1 



l oo 



We can thus plug Eq. (|256j) into Eq. f|254[) to find that at late times (for kr] 3> 1) 



VV 



-3W 



(256) 



(257) 



We will come later to the explicit expression for the term on the R.H.S. of Eq. (|257p . Here it is enough 
to notice that the second-order quadrupole oscillate in time and thus w ! will decay in time as Ti 2 = 1/r] 2 . 
This shows that oj 1 in Eq. (|256p can be in fact neglected with respect to the other terms giving 



,(2)i 



(258) 



It can be useful to compute the combination on the R.H.S. of Eq. (|257p (<5*- — d l dj /V 2 )d}Jl^ k ^ ■ The 
second-order quadrupole moment of the photons in the tightly coupled limit is given by Eq. (|219p . and 



(2)fcj 



d k (v k v j ) - 2v k d j v k 



v j d k v k - v k d j v k 



(259) 



where in the last step we have used that the linear velocity is the gradient of a scalar perturbation. We 
thus find 



"V 2 " 



8 d l 
3V2 



v l d k v k — v k d l v k 

(d k v k ) 2 + vid 3 d k v k - djV k &v k - v k V 2 v k 



(260) 
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Notice that if we split the quadrupole moment into a scalar, vector (divergence-free) and tensor (divergence- 
free and traceless) parts as 

= njpi _ V^nf + nf^' + nf^ + n® k * , (26i) 

o 

then it turns out that 

(*i - ^r) ^ n ? )fcJ = V X 2) * > ( 262 ) 

where II^ 2 ^ is the vector part of the quadrupole moment. Therefore one can rewrite Eq. (|257|) as 

J = -m 2 J drj'U^ . (263) 

7.5.2 Initial conditions for the second-order gravitational potentials 

In order to complete the study of the CMB anisotropics at second-order for modes k 3> k eq we have 
to specify the initial conditions ^f^(0) appearing in Eq. (|253p . These are set on super-horizon scales 
deep in the standard radiation dominated epoch (for rj — > 0) by exploiting the conservation in time of 
the curvature perturbation £. On superhorizon scales (^ 2 ) is given by Eq. (|236p during the radiation 
dominated epoch and, using the (0 — 0)-Einstein equation in the large scale limit Aq 2 ^ = — 2<E>( 2 ) +4<3>( 1 ) 2 , 
we find 

C (2) = -|*(2)( ) - I ($(2)( ) - *( 2 )(0)) + ^ (1)2 (0) . (264) 

The conserved value of is parametrized by = 2cinlC^ 2 ) where, as explained in Section [7.41 the 
parameter onl specifies the level of primordial non-Gaussianity depending on the particular scenario for 
the generation of the cosmological perturbations. On the other hand , during radiation-domination 

$(2) (o) _ ^(2) ( 0) = _q(2) (q) ( (2g5) 

where 

q(2) (0) = -2v- 2 a fc d>( 1 )(o)5^( 1 )(o) + (d^(o)d^(o)) + 9 -n 2 ^nf) , 

(266) 

where we are evaluating Eq. (|B.14p in the limit krj <C 1. The contribution from the second-order 
quadrupole moment in this limit reads 

5« 2 ^<»' = („V - S -3^<'(0)<'(0)sin( fe o„) S in( fac „, 

~> ^*k'( )*k>). P 67 ) 
where F and C are defined in Eqs. (|24ip and (|240p . Therefore we find that in Fourier space 

F 
'¥ 

and from Eq. (I264p we read off the intial condition as (convolution products are understood) 

*( 2 )(0) = f-3(a NL - 1) + ll F(kl ;> k) 1 <> (0)*g>(0) . (269) 



q(2) (0) = 33 ^(0)^(0), (268) 



k 2 
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7.6 Perturbation modes with k <C k 



eq 



Let us consider the photon perturbations which enter the horizon between the equality epoch and the 
recombination epoch, with wavelenghts r]~ l < k < rj^g. In fact, in order to find some analytical solutions, 
we will assume that by the time of recombination the universe is matter dominated r\ eq <C r/* . In this case 
the gravitational potentials are sourced by the dark matter component and their evolution is given in 
Sec IB. 11 At linear order the gravitational potentials remain constant in time, while at second-order they 
are given by Eq. (|B.6l) . In turn the gravitational potentials act as an external force on the CMB photons 
as in the equation (I21ip describing the CMB energy density evolution in the tightly coupled regime. 

For the regime of interest it proves convenient to use the solution of Eq. (|21ip found in f|238[) . The 
source functions 5a and Sy are given by Eqs. (1205 \ and (|210p . respectively. In particular 5a at early 
times - 5a (0) appearing in Eq. ()238j) vanishes. For a matter dominated period 



S A (R = 0) 



3 



(270) 



where we have used the linear evolution equations (|C1| and (|C.5p with $W = SSfV-', and 
S l v (R = 0) = ^ i d j v^ + ]d i A^ 2 -2d i ^ 2 -^d i A^ 



a/ - -a-n {2fe ' . 

4 J 



(271) 



As at linear order we are evaluating these expressions in the limit R = 3/9b/4/9 7 — > 0, while retaining 
a non-vanishing and constant value for R in the expression for the photon-baryon fluid sound speed 
entering in the sines and cosines. In fact this approximation gives the dominant contributions to the 
source terms 5a and S v made by first-order squared terms. Using the linear solutions ()C.15P and ()C.17p 
for the energy density and velocity of photons, the source functions in Fourier space read 



S A (R = 0) 



108 



I ^ J k 2 c s cos(kic s 7]) sm(k 2 c s T]) + ~^k 2 c s sm(k 2 c s rj) 



32 / 9 \ - ki , . oz / y 



32/9 



ki 



k 2 c s sm(kic s rj) 



^(0)^(0), 



(272) 



and 



s l v (R = o) 



g \ 2 jj 

TK c l T~ k 2 sin(fcic s ?7) sin(fc 2 c s 77) 
10 / k\ 



| A;* (-cos{kic s r]) 



9 



18 



6 (u \ 18 
-cos(k 2 c s rj) - — 

5 5 



18 



2ik l ( — - i — k\ -cos(k 2 c s ri) ) - 2oj' l 

10/ 10 2 V5 v " 5 
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+ 



9 \ 2 ok 



3 V 10 



c 2 ^ ■ kiyisin(fci c s rf) sm(k 2 c s i]) 
k 2 k\ 



< } (0)^(0). 



(273) 



In Sy we have used the expression (|219|) for the second-order quadrupole moment niy J of the photons 
in the tight coupling limit, with the velocity Notice that, for the modes crossing the horizon 

at rj > rj eq , we have expressed the gravitational potential during the matter dominated period in terms 
of the initial value on superhorizon scales deep in the radiation dominated epoch as VE^ 1 ) = 9^ 1 )(0)/10. 

As for the second-order gravitational potentials we have to compute the combination <1?( 2 ) + \p( 2 ) 
appearing in Eq. (I238p . The gravitational potential \I^ 2 ) is given by Eq. (|B.6j) . while $^ 2 ^ is given by 



<J>(2) = ^(2) _ q(2) ; 

according to the relation (j33|) . where for a matter-dominated period 



We thus find 



q( 2 ) = 5V- 4 d i 6» i (dM 1 )d^W) - -v~ 2 (d k ^d k ^) 

3 



$(2) + ^(2) = _5v- 4 a i ^'(a i * (1) 3j* (1) ) + 5v- 2 (^* (1) s** (1) ; 

3 



(274) 



(275) 



(276) 



which in Fourier space reads 

$(2) _|_ \I/( 2 ) 



2tt£>(0) + 



9 



io ) *£(o)*£(°) 



jG(k 1 ,ii 2 ,-k)rf i -^F(k 1 ,-k 2 ,-k) 



where the kernels of the convolutions are given by Eq. (|242p and 

10 (k-k 1 )(k-k 2 ) 



G(ki,k 2 ,k) = ki-k 2 



k 2 



(277) 



(278) 



In Eq. (|276p (0) is the initial condition for the gravitational potential fixed at some time rfi > r\ eq . 
For the regime of interest it corresponds to the value of the gravitational potential on superhorizon scales 
during the matter-dominated epoch. 

Notice a property that will be useful later on. By looking at Eqs. (I275j) and the explicit solution 
for ^( 2 ) (|B~6l) which grows as rf, it is easy to realize that on very small scales, for kr] S> 1, the two 
gravitational potentials are equal with 



$( 2 ),f( 2 ) = ^(0) + iG(k 1) k 2 ,k) )) 2 . 



(279) 
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We are now able to compute the integrals entering in the solution (|238p . The one involving the 
second-order gravitional potentials is straightforward to compute 



- Akc s jT drf (V 2) + 1- (2) + iM> {2) ) sm[kc s (r] - rf)] 

- 8^)(0)(l-cos(fcc s 77))-^ 



1 



G (-2 + {kc s rf) 2 + 2cos(fcc s 7?) 



7k 2 c 2 s 

4R*^(0) (1 - cos(kc s 7])) - 4R 
1 



< } (0)*£(0) 



5F 

--p- (1 - cos(kc s rj)) 



^- 2 G (-2 + (kc sV ) 2 + 2cos(/cc s r ? ))l (±Y <>(0)*«(0) 



(280) 



where for brevity F and G stand for the convolutions (|242p and ()278p . Notice also the we have isolated 
the terms proportional to R as it will be useful later on. 

For the two remaining integrals, in the following we will show only the terms that in the final expression 
for Aqq 1 and the second-order velocity utf give the dominant contributions for kr] 3> 1, even though we 
have perfomed a fully computation. The integral over the source function 5a yields a sum of oscillating 
functions (cosines) which turn out to be subdominant, so we do not display the full result. For the last 
integral we find ((1 «->■ 2) stands by an exchange of indices) 



4 iki 
3 kcq 



drj' Sy sm[kc s (r) - rf)} 



27 2k • k 2 + k 1 
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k 2 c 2 



+ (!<-> 2) 



^£(0)^(0) 



(281) 



where the terms that have been dropped vary in time as sines and cosines. We have written the contri- 
bution in Eq. (|28ip because, upon integration over time, it will give a non-negligible contribution to the 

velocity v^ z . 



From the general solution f|238|) and the expression ()B.6P for the second-order gravitational potential 
1 I'( 2 ) we thus obtain 



A 



(2) 
00 



-4(l + i?) tf£)(0) 

2(9a NL - 7)^(0)^(0) + (8 + 4i2)*£)(0)l cos^r?) 



(282) 



9 V,T,(D, 



^G(k 1 ,k 4> k)i/ , (l + fl) ( *£(0)<>). 



We warn the reader that in writing Eq. (|282p we have kept all those terms that contain the primordial 
non-Gaussianity parametrized by onl> and the terms which dominate at late times for kr] S> 1. 



7.6.1 Initial conditions for the second-order gravitational potentials 

The initial condition ^fm\o) for the modes that cross the horizon after the equality epoch is fixed by 
the value of the gravitational potential on superhorizon scales during the matter dominated epoch. To 
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compute this value we use the conservation on superhorizon scales of the curvature perturbation £( 2 ) 
defined in Eq. (|37|) . For a matter-dominated period the curvature perturbation on large-scales turns out 
ot be 

C ( 2 ) = -^(0) + i^ + f^) 2 (0), (283) 

6 Pm y 

where we used the energy continuity equation S^p' m + 3H5^ p m — 3p m ^^' = and the (0 — 0) Einstein 
equation <$W p m j ' p m = —2^^ in the superhorizon limit. 

From the (0 — 0) Einstein equation on large scales 5^ p m j p m = — 2<£>( 2 ) + 4<I>( 1 ) 2 bringing 



C (2) = _jj*£) (0) _ 2 ^ )(0) _ ^ (0) ) + ^W2 (0) 



(284) 



The conserved value of £( 2 ) is parametrized as in Eq. p]), C (2) = 2a NL C (1)2 = (50a NL /9)^ (1)2 , with 
£W = — 5\l/W/3 on large scales after the equality epoch. At second-order the two gravitational potentials 
in a matter dominated epoch differ according to Eq. (|275p and using Eq. (|284j) we find 



= -^(a Nh - 1) (^(O)) 2 + ^ [2V" 4 ^(a^«(0) 

x s,-*W(o)) - ^v- 2 (a fc ^ (1) (o)a fe * (1) (o)) 



(285) 



we have expressed the gravitational potential during the matter dominated period in terms of the 
initial value on superhorizon scales after the equality epoch as = 9^'^ 1 ^(0)/10. In Fourier space 
Eq. (|285p becomes 



where .F is the kernel defined in Eq. (|242p . 



27 /9\ 2 F(k 1; k 2 ,k) 
_( ONL -l) + 2(-l p 



< } (0)^(0), 
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We can use the explicit expression for *$(0) in Eq. 



still keeping only the terms that contain 



the primordial non-Gaussianity parametrized by onL; and the terms which dominate at late times for 
kr\ ~3> 1 to find 
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{l + R)(a NL -l) 
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- -r(9a N L - ^cos^c^) + — i?(a NL - l)cos(kc s r]) 



|(1) G(k 1 ,k 2 ,k) ?? 2 ( 



(0) 



(1) 



(287) 



7.6.2 Second-order photon velocity perturbation 

The second-order velocity of the photons can be obtained from Eq. (|209j) 



,(2)< 
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Notice that for simplicity in writing Eq. (|288p we have dropped off the dependence on R. In fact 
the main conclusions of this Subsection remains unchanged. The second-order gravitional potential in 
matter-dominated universe can be obtained from Eqs. (|274p - (|275p and Eq. (]B.6P as 



$(2) = #(2)(o) -If d k yMd k vW _ ^v-idi&id^d^W)) rf 



(289) 



In Fourier space this becomes 

$(2) = *£)(()) + 



^G(k 1 ,k 2 ,k)7 ? 2 --AF(k 1 ,k 2 ,k) 



(290) 



where the kernels of the convolutions are given by Eqs. (|242p and ()278p . The integral over $( 2 ) in Eq. (1288]) 
is then easily computed 
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10 



0)^(0) 
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where as usual the equivalence symbol means that we are evaluating a given expression in Fourier space. 
For the integral over the source function S v we use its expression in Fourier space, Eq. (|273p . and the 
dominant terms for krj S> 1 are 



dr/S 



v 



9 



.■81 



2ik l — + ik l — +i- — — (fcf - fcf )k • kifc 
V 10/ 25 3 V 10/ A; 4 2 1; 



\ 2 1 



(292) 



Notice that, in order to compute this integral, we must know the second-order vector metric petturbation 
aA This is easily obtained for a matter-dominated universe from Eq. (|B.7p . Using Eqs. (|B.4p and (|B.2p 
one finds 



- (- ) v-a, 

3V10' J 



dV* (1) (0)d% (1) (0) - cVV 2 ^ (1) (0) 



x 5^ (1) (0) 



(293) 



giving rise to the third term in Eq. ()292p . 

Finally for the integral over Aqq some caution is nedeed. Since in the final expression for the 

(2) 

dominant terms at late times turn out to be proportional 77, one has to use an expression for A$ that 
keep track of all those contributions that, upon integration, scale like 77. Thus we must use the expression 
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written in Eq. (|282p . plus Eq. (|28ip . and some terms of Eq. (|280p that have been previously neglected in 
Eq. (1282]) . Then we find for krj > 1 



+ 4 



1 f'^'^A(2) - ik " 



dr/ffA™ = ^ -W&H0)r, + (2(9a NL - 7)^(0)^(0) 



+ 8*W(0) 



s'm(kc s rj) 
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_9_\ 2 G_ 54 fc 2 + k-(k!+k 2 
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Using Eqs. (p9TD . ([292]) and (gMt we get 
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/9\ 

To obtain Eq. ()295|) we have also used the explicit expression ([286]) for ^ ; (0) and we have kept the 
terms depending on onl parametrizing the primordial non-Gaussianity and the terms that dominate at 
late times for ki] S> 1. 



7.7 Perturbation modes with k ^> k eq : Cold Dark Matter perturba- 
tions at second order and improved analytical solutions for CMB 
anisotropies 

In Sec. 17.51 we have computed the perturbations of the CBM photons at last scattering for the modes 
that cross the horizon at rj < rj eq under the approximation that the universe is radiation-dominated. 
However around the equality epoch, through recombination, the dark matter component will start to 
dominate. In this section we will account for its contribution to the gravitational potential and for the 
resulting perturbations of the photons from the equality epoch onwards. This leads to a more realistic 
and accurate abalytical solutions for the acoustic oscillations of the photon-baryon fluid for the modes of 
interest. 

The starting point is to consider the density perturbation in the dark matter component for subhorizon 
modes during the radiation dominated epoch. Its value at the equality epoch will fix the magnitude of 
the gravitational potential at rj eq and hence the initial conditions for the subsequent evolution of the 
photons fluctuations during the matter dominated period. At linear order the procedure is standard (see, 
e.g |42|). and we will use a similar one at second-order in the perturbations. Thus this section serves 
also as a guide through the evolution of the CDM density perturbations at second-order accounting for 
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those modes that enter the horizon during the radiation dominated epoch. This allows to determine the 
second-order transfer function for the density perturbations with a generalization at second-orde of the 
Meszaros effect. 



7.7.1 Subhorizon evolution of CDM perturbations for r\ < r) eq 

From the energy and velocity continuity equations for CDM it is possible to isolate an evolution equation 
for the density perturbation 5 d = Sp^/pd, where the subscript d stands for cold dark matter. In Ref. [3S] 
we have obtained the Boltzmann equations up to second-order for CDM. The number density of CDM 
evolves according to [38] 

dn± + e *+*^M + m _ _ n d 

+e*+*$ jk v k d n d = 0. (296) 
At linear order n d = n d + 5^n d and one recovers the usual energy continuity equation 

# } ' +«g i -3*W.* = 0, (297) 
with 5^ = SWpd/pd = 5 w n d /n d . The CDM velocity at the same order of perturbation obeys [38] 



it 



(i^ + ^M* = (298) 



Perturbing tt-cdm up to second-order we find 

*f+^-3* (2) ' = -2(*« + *«)t,^-2t,^-2t,f)*dg 

+ e^ 1 ^ + (wf - 2$f)v<p k . (299) 

The R.H.S. of this equation can be further manipulated by using the linear equation (|297p to replace 
v d i yielding 



jt(2)' , (2)t 
S d + V d A 



+ 2^%f\ (300) 



where we use <3?( 1 ) = Vl^ 1 ). In Ref. [38] the evolution equation for the second-order CDM velocity 
perturbation has been already obtained 

v^ + Hvf* + + 2HJ + ^ = 2^'vf 1 -2v [ d 1)] d. J v d 1)i 

- 4$W$W.i. (301) 

At linear order we can take the divergence of Eq. (|298p and, using Eq. (|297p to replace the velocity 
perturbation, we obtain a differental equation for the CDM density contrast 

a ( 3 ^(i)' -8W)]' = -aV 2 $W, (302) 
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which can be rewritten as 

+ K5 { j Y = S« , (303) 

where 

5(1) = 3^(1)" + 37^(1)' + v 2 $« . (304) 

When the radiation is dominating the gravitational potential is mainly due to the perturbations in the 
photons, and a(rj) oc r\. For subhorizon scales Eq. (I303j) can be solved following the procedure introduced 
in Ref. [?]. Using the Green method the general solution to Eq. (|303j) (in Fourier space) is given by 

5 { }\k,ri) =d + C 2 Hv) ~ r dr/S^i^T/ihxikrf) - ln(krj)) , (305) 

Jo 

where the first two terms correspond to the solution of the homogeneous equation. At early times the 
density contrast is constant with 

4 1} (0) = |^(0) = -1*^(0), (306) 
having used the adiabaticity condition, and thus we fix the integration constant as 

C 1 = -3$£ ) (0)/2, (307) 

and C2 = 0. The gravitational potential during the radiation-dominated epoch starts to decay as a given 
mode enters the horizon. Therefore the source term behaves in a similar manner and this implies 
that the integrals over rj reach asymptotically a constant value. Once the mode has crossed the horizon 
we can thus write the solution as 

= A^^{0)ln[B^k V ], (308) 
where the constants and B^ 1 ' are defined as 

roo 

AW$W(0)= / drfS^\r/)rf, (309) 



and 

O POO 

A(i)$(°)ln(.B(i)) = _2$(i)(o)_ / drfS w (Tf)rf]n(hf). (310) 
2 Jo 

The upper limit of the integrals can be taken to infinity because the main contribution comes form when 
krj ~ 1 and once the mode has entered the horizon the result will change by a very small quantity. 
Performing the integrals in Eq. (j309]) and ([310]) one finds that A^ = -9.6 and ~ 0.44. 

Before moving to the second-order case, a useful quantity to compute is the CDM velocity in a 
radiation dominated epoch. From Eq. (|298j) it is given by 

iff* - 1 = -m^ (1 \o) kCsV ; 3 ff Csr,) , (3ii) 

a Jo k 6 c 6 s rf 

where the last equality holds in Fourier space and we have used Eq. (IB.llj) (and the fact that a(rj) oc r/ 
when radiation dominates). 
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Combining Eq. (|300p and (|301[) we get the analogue of Eq. f|303|) at second-order in perturbation 
theory 

4 2) " + ^4 2) ' = s(2) > ( 3l2 ) 

where the source function is 

5 (2) = 3^(2)" + 3^(2)' +V 2 $ (2)_ 2^(*«4 1}i )+V 2 4 1)2 

+ 2V 2 $( 1 ) 2 + 1 [a (^W'tfW - 6 (*« 2 )' + (4 1)2 )' 

" + 2*ff«r)] '■ (313) 

In fact we write Eq. (|312|) in a more convenient way as 

S f _ 3^(2)" _ s > + n{6 f _ 3^(2)' _ Sl) = S2 ; (314) 

where for simplicity we have introduced the two functions 

Sl = 44 1} '*« - 6 (^ 2 )' + (4 1)2 )' - + 2*«„?>* , (315) 

and 

S2 = VV) _ 28^' v^) + V 2 4 1)2 + 2V 2 ^ 2 . (316) 

In this way we get an equation of the same form as (|3U3p in the variable [5^ — — drj'si(rj')] with 
source S2 on the R.H.S.. Its solution in Fourier space therefore is just as Eq. ()305p 

5f -&BV> - T drfaxtf) = Ci + C 2 ln(r ? ) (317) 
J o 

rv 

— / dr]'s2(r]') rf [ln(fcr/) — ln(fcr/)] . 

As we will see, Eq. (|317|) provides the generalization of the Meszaros effect at second-order in perturbation 
theory. 

7.7.2 Initial conditions 

In the next two sections we will compute explicitly the expression (|317|) for the second-order CDM density 
contrast on subhorizon scales during the radiation dominate era. First let us fix the constants C\ and C2 
by appealing to the initial conditions. At 77 — > the L.H.S. of Eq. (|317|) is constant, as one can check by 
using the results of Sec. 17.5.21 and the condition of adiabaticity at second-order (see, e.g., Ref. [?,[5]) which 
relates the CDM density contrast at early times on superhorizon scales to the energy density fluctuations 
of photons by 

4 2) (0) = |a8>(0) - i (^(O)) 2 = f A$(0) - \ (^(0)) 2 , (318) 
where in the last step we have used Eq. (|306p . Therefore we can fix C2 = and 

d = 4 2) (°) -3^ (2) (0). (319) 
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Eq. (|237p gives Aqq^(O) — A^^(0) in terms of the primordial non-Gaussianity parametrized by onl 5 
and the expression for ty( 2 \o) have been already computed in Eq. (|269l) . Thus we find (in Fourier space) 



A 



(2) 
00 



2(3a NL - 1) + 44 



F(ki,k 2 ,k) 

k 2 



<>)<>), 



and from Eqs. ()318j) we derive the initial density contrast for CDM at second-order 



<# } (o) 



3 QQ F(k 1; k 2 ,k) 3 
-(3a NL - 1) +33 p 



0)0) 



Eq. (|32ip togheter with Eq. (|269[) allows to compute the constant Ci as 



d = 4 2) (0) -3^(0) 



27 / ,N 9 

y(«NL - 1) + ^ 



OK>) 



(320) 



(321) 



(322) 



7.7.3 Meszaros effect at second order 



We now compute the integrals over the functions s± and s 2 appearing in Eq. (|317p . Let us first focus on 
the integral Jq drf s\{rf). 

Notice that, using the linear equations (|297p and (|298p for the CDM density and velocity perturba- 
tions, the function s\{if) can be written in a more convenient way as 



and then 



(323) 



(324) 



+ / -2^ 1) ^+2(^^ K ), fc -6^ 1 ^ 



In Eq. (|324|) all the quantities are known being first-order perturbations: the linear gravitional potential 
for a radiation dominated era is given in Eq. (jB.lip . the CDM velocity perturbation corresponds to 
Eq. (|311|) and the CDM density contrast is given by Eq. (|308|) . Thus the integral in Eq. (|324p reads (in 
Fourier space) 



drf 



+ (9(k • ki) - 27k( 



2 -.k 1 c s if - sin{kic s rf) 



kfCgT]' 2 



sm(k s c s 7]') - k2C s rj'cos(k2C s r]') 



k^CgT]' 3 



<V (0)^(0) 



(325) 



Let us recall that we are interested in the evolution of the CDM second-order density contrast on sub- 
horizon scales during the radiation dominated epoch. Therefore once we compute the integrals we are 
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interested in the limit of their expression for late times (krj S> 1)- For the first contribution to Eq. f|325|) 
we find that at late times it is well approximated by the expression 



JO k^C s Tj 



2.2 



1.2 



ln(feic s 7?)] 2 + ln(B^k 2 v) Hhc s rj) 



(326) 



We have computed also the remaining integral in Eq. (|325p . but it turns out to be negligible compared 
to Eq. (|326p . The reason is that the integrand oscillates with an amplitude decaying in time as ?]~ 3 , 
and thus it leads just to a constant (the argument is the same we used at linear order to compute the 
integrals in Eq. (|305f) ) . Thus we can write 



i) ki ■ k 2 



2.2 (~ [\n{k x c s ri)] 2 + \n{B^k 2 r,) ln(fcic s Tj) 



(327) 



< } (0)*g(0). 



We now compute the integrals over the function s 2 (?7) given in Eq. ()316|) . Since at late times (p^ 2 ~ 
1/r/ 4 and (fW'tij ),i ~ l/rp the main contribution to the integral will come from the two remaining 
terms, <3?( 2 ) and v^ 2 , whose amplitudes scale at late times as l/?? 2 



s 2 ~ V 2 $ {2) + V 



2 (1)2 



(328) 



Two are the integrals that we have to compute 

rn 



o 



and the one multiplying ln(fcr/) 



df] S2(t))t] ki(kr] ) 



drf 82(rf)rf 



(329) 



(330) 



Let us first consider the contributions from V 2 i>!j . The second integral is easily computed using the 
expression (|31 1 1) for the linear CDM velocity. We find that at late times the dominant term is (for 
krj > 1) 



9 , 2 k r k 2 

h 2 h 2 



-k 



]n(krj) 



< } (0)*{^(0). 



(331) 



The first integral can be computed by making the following approximation: we split the integral between 
< krj < 1 and krj > 1 and for < krj < 1 we use the asymptotic expression 



,(!)* 



1 



(332) 
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while for krj > 1 we use the limit 



,(!)* 



3 k l 1 T (ii , , 

*k (°) 



c 2 k krj 



(krj > 1) . 



(333) 



The the integral for < kr] < 1 just gives a constant, while the integral for krj > 1 brings the dominant 
contribution at late times being proportional to [ln(A:r/)] 2 so that we can write (krj S> 1) 

jT d/Vt^V = ^| KM] 2 *k>) *k>) • (334) 

As far as the contribution to the integrals (|329|) and (|330p due to V 2 ^ 2 ** is concerned we have just to 
keep track of the initial condition provided by the primordial non-Gaussianity. We have verified that all 
the other terms give a negligible contribution. This is easy to understand: the integrand function on 
large scale is a constant while at late times it oscillates with decreasing amplitudes as r]~ 2 , and thus the 
integrals will tend asymptotically to a constant. We find that 



dr{V 2 & 2) rf ~ -9$ {2) (0) 



and 



dr/V 2 $ (2 V ln(/b/) 



-9 + 97-9- 



ln3 



$ (2) (0), 



(335) 
(336) 



where 7 = 0.577... is the Euler constant, and <3?( 2 )(0) is given by Eq. (I265D . 
Therefore, from Eqs. (l3Mj) . (133TT) . and (13351)- (1336"]) we find that for kr] > 1 



ki • k? 



rrj g 
dr]'s 2 {v'W[H kr l') ~ m(fcr/)] = -^-jfc 
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[ln(A;77)] 2 
In(jfe»7)*W(0) 



(337) 



Let us collect the results of Eqs. (|322|) . (13271) and (l337|) into Eq. (j3T7|) . We find that for kr] > 1 



<5i 2) (fcr ? > 1) 



3(a NL - l)Ai ln(BiA;7/) + ^ 2 m^^r?) ln^A^) 



+ 



-A 



ki ■ k 2 



1.2 



x ln(fcic s 7 ? )J + (1 -H- 2) 
x < } (0)< } (0). 



2.2 - — []n(k lCs rj)Y + ln^fer/) 



+ 2? fc2 W^ 2 



(338) 



Notice that in Eq. (|317p we have neglected \I'( 2 ), which decays on subhorizon scales during the radiation 
dominated epoch (see Eq. (|247j) . and we have used Eqs. (j306H and (j308[) . Eq. ()338p represents the second- 
order Meszaros effect: the CDM density contrast on small scales (inside the horizon) slowly grows starting 
from the initial conditions that, at second-order, are set by the primordial non-Gaussianity parameter 
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onl- As one could have guessed the primordial non-Gaussianity is just transferred linearly. The other 
terms scale in time as a logarithm squared. We stress that the computation of these terms allows one to 
derive the full transfer function for the matter perturbations at second order accounting for the dominant 
second-order corrections. In the next section we will use (j338H to fix the initial conditions for the evolution 
on subhorizon scales of the photons density fluctuations Aq after the equality epoch. 

7.7.4 Second-order CMB anisotropies for modes crossing the horizon during the 
radiaton epoch 

(2) 

In this section we derive the energy density perturbations Aq of the photons during the matter dominated 
epoch, for the modes that cross the horizon before equality. In Sec. l7.6l we have already solved the problem 
assuming matter domination for modes crossing the horizon after equality. Thus it is sufficient to take 
the result (|282p and replace the initial conditions 

A$ = -4(1 + J R)^)(0) (339) 



+ 



A + (8 + 4fl)*^)(0)l cos(kc s r]) +Bsm(kc s r)) 



9 



2 



- =G(k ll k a ,k)^(l + /J)^j <>)<>). 
where we have restored the generic integration constants A and B, \pW is the linear gravitational potential 

(2) 

(which is constant for the matter era) and V&m (0) represents the initial condition for the second-order 
gravitational potential fixed at some time rji > rj eq . Eq. (|338p allows to fix the proper initial conditions 
for the gravitational potentials on subhorizon scales (accounting for the fact that around the equality 
epoch they are mainly determined by the CDM density perturbations). At linear order this is achieved 
by solving the equation for 5^p which is obtained from Eq. (I302|) and the (0 — 0)-Einstein equation which 
reads (see Eq. ([50])) 

3W* «' + mW> - v 2 \p« = -\n 2 {^bf + £ia$) . (340) 



+H#> = (341) 



On small scales one neglects the time derivatives of the gravitational potential in Eqs. (|302p and (|340p 
to obtain 

where we have also dropped the contribution to the gravitational potential from the radiation component. 
The solution of this equation is matched to the value that 5^ has during the radiation dominated epoch 
on subhorizon scales, Eq. (|308p . and one finds that for rj S> rj eq on subhorizon scales the gravitational 
potential remains constant with 



We skip the details of the derivation of Eq. (|342p since it is a standard computation that the reader can 
find, for example, in Refs. [32] • Since around r\ eq the dark matter begins to dominate, an approximation 
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to the result (|342p can be simply achieved by requiring that during matter domination the gravitational 
potential remains constant to a value determined by the density contrast (|308|) at the equality epoch 



V W Wleq — 2 n °d IVeq I 



from Eq. f|34Q|) on small scales, leading to 



K (V > Veq) 



(kr), 



eq I 



(1) _ HgjfaTeg) (i) 

,2^ he, " ( . 13fc77e?) 2 W k W> 



(343) 



(344) 



where we used a(r?) oc rf during matter domination and Eq. (|308p with = —9.6 and B\ = 0.44. 

At second-order we follow a similar approximation. The general solution for the evolution of the the 
second-order gravitational potential for rj > rj eq is given by Eq. (|B.6p . We have to determine the 
initial conditions for those modes that cross the horizon during the radiation epoch. The (0 — 0)-Einstein 
equation reads 

% n * ( e± s f + ^a^) = sn^y + m 2 ^ - v 2 ^ 



- m 2 - im*v>*w - 3 (*«') + d&MtfvW 

- 4$(i)v 2 $( 1 ) . (345) 
We fix the initial conditions with the matching at equality (neglecting the radiation component) 

(346) 



V 2 * (2) - d^Wd^M +4#Wv 2 *«| r?eg ~ 



3 a/2x( 2 )| 

2 d Weq I 



(2) 

where for small scales we neglected the time derivatives in Eq. (|345p . Using Eq. (|338p to evaluate 5^ L 
and Eq. (|342|) to evaluate ^f^\r] eg ) we find in Fourier space 



<:<-)) 



3(a NL -l)MS^- 



ki • k 2 
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ln(0.15fci7? eq ) 



(0.l3k Veq ) 2 



(Q.27^^) 2 



(0.13^) 



ki • k 2 [ln(/c7? eg )]_ 3 ki • k 2 ^ 2 



fc 2 /c 2 (A;?/ 



2 c 2 fc 2 



1.2 [ln(feic s r/ eg )]' 



2 (0.l3k Veq ) 2 



(347) 



In Eq. (p39~P the initial condition *^(0) is given by Eq. flMZD and *W is given by Eq. ([532]) . The 
integration constants can be fixed by comparing at ~ ?7 e g the oscillating part of Eq. (1339P to the 
solution Aq obtained for modes crossing the horizon before equality and for r] < r] eq , Eq. (|253p . Thus 
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-4(1 + R)^ 2 \r] eq ) + Acos{kc s r]) 



2 -(l + ^Gi^M^W^SiVe^M , 



(348) 



where 



A 




x 



cos(k 2 c s rj eq ) - 2kik 2 cos(kc s r) eq ) + {k\ + k\ - k 2 ) sva(kic s rj eq ) 



x 




(349) 



and tf( 2 >(0) is given in Eq. (|269p . 

8 Secondary effects and contamination to primordial NG 

Given that a detection of a sizable primordial non-Gaussianity (and its shape) would represent a real 
breakthrough into the understanding of the dynamics of the universe during its very first stages, it is 
crucial that all sources of contamination to the primordial signal are well understood and kept under 
control. In fact any non-linearities can make initially Gaussian perturbations non-Gaussian. Such non- 
primordial effects can thus complicate the extraction of the primordial non-Gaussianity: we have to be 
sure we are not ascribing a primordial origin to a signal that is extracted from the CMB (or LSS) data 
using estimators of non-Gaussianity when that signal has a different origin. Moreover, as stressed in 
Section. 17. 1[ we must always specify of which primordial non-Gaussianity we study the contamination 
from the non-primordial sources (e.g. primordial non-Gaussianity of "local" , "equilateral" or "folded" 
shape). 

Broadly speaking, non-primordial sources of non-Gaussianity can be classified into four categories: 
instrumental systematic effects; residual foregrounds and unresolved point sources; some well known 
secondary CMB anisotropics, such as Sunyaev-Zel'dovich (SZ) effect, gravitational lensing, Rees-Sciama 
effect, and finally previously unknown effects coming from non-linearities in the Boltzmann equations, 
which are related to the non-linear nature of General Relativity and to the non-linear dynamics of 
the photon-baryon system. This Review focuses on the last category, but in this Section we offer also a 
summary of some recent results of the other types of possible contaminations. Before doing that, however, 
it is important to be more precise about how secondary effects may impact the extraction of a primordial 
non-Gaussian signal. Mainly they act in two ways. They might "mimic" a three-point correlation function 
similar in shape to the primordial one. This would produce a bias or a contamination to the estimator 
of primordial non-Gaussianity. On the other hand secondary effects might increase the variance of the 
estimator without contributing to the signal-to-noise ratio. In other words, in this case they degrade the 
signal-to-noise ratio. Let us define in a quantitative way how to characterize these effects. 
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8.1 Signal-to- noise ratio, shapes, and contamination to primordial 
non- Gaussianity 

The angular bispectrum, B^i^, the harmonic transform of the angular three-point function, of the CMB 
anisotropies is often used to measure non-Gaussianity (see for example [5]). Prom the usual spherical 
harmonic coefficients of the temperature anisotropies 

a tm = J rf 2 nF; m (n)^ , (350) 

the angle-averaged bispectrum is given by (for more details see, e.g., Ref. |35j ) 

B h£2i 3 =Yl ( ^ ^ ) Kmi%mj a 4m 3 ) , (351) 

all rn \ / 

where the brackets indicate ensemble average. 

We shall quantify the degree to which the primordial and the secondary bispectra are correlated, as 
well as their expected signal-to-noise ratio following a method that as now become standard [35 1 1331165] . 
Namely, in the limit of weak non-Gaussianity, one can introduce the Fisher matrix for the amplitudes of 
the bispectra, F^, given by 

Fij = £ g^M , (352) 
where the variance of the bispectrum is 



rr 

2<e 1 <e 2 <£ 3 u t\t2h 



G \hh = ( B hhi 3 ) ~ ( B htih) ~ C£ 1 Ci 2 Ce 3 Ai 1 i 2 £ 3 , (353) 

and A^^g takes values 1, 2, and 6 when all ts are different, two of them are equal and all are the same, 
respectively. The power spectrum, Ce, is the sum of the theoretical CMB and the detector noise. 
The signal-to-noise ratio is given by 

^ ' = , (354) 

FH 1 



N 

and we define the cross-correlation coefficient between different shapes i and j, fij, as 

r v = / Ji i ' ( 355 ) 

a 33 

which does not depend on the amplitudes of the different bispectra, but on the shapes, thus measuring 
how similar are two bispectra i and j. One can also define a degradation parameter di = FaF^ 1 for a 
degraded (S/N) due to r^. 

If a secondary bispectra has some correlation with the primordial one, one expects that the signal-to- 
noise ratio will be degraded: if one does not account for any secondary bispectra (S/N)q = ^/-Fprim.prim 
while, marginalizing over a (single) secondary bispectrum, (S/N) pv i m gets modified from its zero-order 
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value to (S/N) pr i m = (S/N)q(1 — r 2 ^) 1 / 2 . This means that the minimum detectable value of the primordial 
non-linearity parameter /nl which is obtained by imposing (S/N) pT i m = 1 (or the 1-a uncertainty on 

/nl given by 5/ NL = ^/(^)p r L,priml/ N L=i) g ets shifted by a quantity A/ NL /(/ NL ) = (1 - r^)" 1 / 2 - 1. 

How much does a given secondary bispectra contaminate the extraction of the primordial bispectrum? 
If, for example, the predicted shape of the secondary bispectra is sufficiently different from that of 
the primordial bispectrum, then one would hope that the contamination would be minimal. We can 
quantify the contamination of the primordial bispectrum as follows: we fit the primordial bispectrum 
template to the second-order bispectrum, and find the best-fitting /™£ ( "con" stands for contamination) 
by minimizing the x 2 given by 

f*TT uVrim _ R 2nd \ 2 



X = }^ 2 ' ( 356 ) 

2<£ 1 <£ 2 <£ s titzh 

with respect to /nl- Here, Bf^^ is the primordial bispectrum with amplitude /nl = 1 [35]. We obtain 



/] 



con 



i r>2nd r>prim 



iy 2<e 1 <e 2 <h 

B p - Hm \ 2 



N = £ ^ ■ (357) 

2<£ 1 <i 2 <h a hhh 

This is the value of /nl one would find, if one did not know that the primordial bispectrum did not exist 
but there was only the secondary bispectrum. In other words the effective non-linearity parameter 
gives the value that a bispectrum estimator designed for constraining primordial non-Gaussianity would 
measure due to the presence of the secondary bispectrum B 2 ^^. As one can see turns out to be 
proportional to the mixed entry of the Fisher matrix, since this governs the correlation between different 
types of bispectra on the basis of their shape dependence. 



8.2 Some "well-known" secondary sources of non-Gaussianity 

Instrumental systematic effects are one of the principal issues in the search for primordial non-Gaussianity. 
Since they depend on the specific instrument, here we do not intend to go into any detail, but we refer 
the reader to Ref. [60] as an interesting example of this kind of effects for the currently flying mission 
of the satellite Planck. Residual foregrounds and unresolved point sources represent other important 
possible contaminants and particular attention have been devoted to them in the analysis of the WMAP 
data [591 ESI HE]- For an experiment like WMAP secondary effects such as Sunyaev-Zel'dovich effect, 
gravitational lensing, Integrated-Sachs- Wolfe (or Rees-Sciama) effect should be not so relevant for £ < 500 
with a bias of just 1.5 to 2 |70j for local primordial non-Gaussianity (see also|35j and [36]). El However 
for higher resolution and more sensitive experiments like Planck, such effects must be taken into account 

5 In order to understand the relative importance of the bias from a given secondary effect let us recall that, 
including the information from polarization, the minimum detectable (local) for an experiment like WAMP is 
of the order of 10 while for Planck it is of the order of 3; the minimum detectable value for an equilateral primordial 
non-Gaussianity is J^l — 30 for Planck, see, e.g., Refs. [35 ] l80 l l4l [58 ] , 
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with care. Here we report the results of some recent analyses about the contamination to primordial 
non-Gaussianity from some secondary effects. Such a brief summary is by no means exhaustive. An 
interesting case is given by the cross-correlation of the ISW and the lensing effects, see Ref. [33] . In the 
detailed analysis of [7T] it is shown that the ISW-lensing bispectrum produces a bias of /ni?* — 9 to 
local primordial non-Gaussianity (see also [70]), while the bias to equilateral primordial non-Gaussianity 
is negligible (see also [361 170]). The impact is mainly on the local type of primordial non-Gaussianity 
because the ISW-lensing correlation correlates the large-scale gravitational potential fluctuations sourcing 
the ISW effect with the small scale lensing effects of the CMB, thus producing a bispectrum which peaks 
in the squeezed configurations as the local shape. As far as the increase of the variance of the estimator of 
primordial non-Gaussianity induced by CMB lensing is concerned it turns out to be of the order of 20% 
for an experiment like Planck (see Refs. [711 180] ). Notice that, interestingly enough, gravitational lensing 
can also have another peculiar effect on the extraction of primordial non-Gaussianity: it can modify its 
shape by smoothing its acoustic features. However, according to the analysis of Ref. [71] (and contrary 
to the results of [61]) such a "distorsion" of the shape should be small, of the order of 10% for I < 2000. 

A similar analysis has been performed looking also at scales where non-linear effects can become 
relevant. In this case a bispectrum from the correlation of the Rees-Sciama and lensing effects has been 
studied in Ref. [63], and it has been found that in this case the bias to a local primordial non-Gaussianity 
corresponds to f^ t — 10, agreeing with the results of [71] where the two analyses can be compared. 

Another secondary source of contamination that has been studied quite recently in Ref. [62] deals with 
bispectra generated by correlations of number density and lensing magnification of radio and SZ point 
sources with the ISW effect. Also in this case the large-scale modulation of small-scale number density 
due to fluctuations which source the ISW generates a bispectrum which peaks on squeezed configurations 
and it has been found that it corresponds to a contamination of local primordial non-Gaussianity of 
jcont ^ y .5 ; which should be not so relevant for an experiment like Planck. The novelty of this effect is 
that it is different from the usual bispectrum of point sources when treated with a Poisson distribution 
(see, e.g., [35]). 

Other secondary effects include, for example, the SZ-SZ-SZ bispectrum [6l] and non-Gaussianities 
from the kinetic SZ and Ostriker-Vishniac effect [661 165] , However the SZ-SZ-SZ bispectrum can be 
considered as an extra Poisson contribution to the unresolved point sources for I < 1500 |70| . Moreover 
recently an analysis of the bispectrum generated from inhomogeneous reionization has been performed 
in [67] . where it has been shown that it can give rise to a very small contamination to the local primordial 
non-Gaussianity of f$£^ ~ —0.1. 

Finally let us also recall some recent studies of a specific effect arising at second-order in perturbation 
theory, when looking at the fluctuations in the Boltzmann equations for the photon-baryon fluid. In 
Refs. [72} [68l 173"! 173] the perturbations to the phase of recombination between electrons and baryons has 
been considered ("inhomogeneous recombination") . It turns out that the electron density perturbation 
can be of the order of 5 times bigger than the one in the density of the baryons. This gives rise to the 
prospect of a non-Gaussianity which corresponds to | /nl 1 * I — 5- In fact the detailed study of this effect 
shows that the contamination is much smaller. The bispectrum from the perturbed recombination peaks 
in the squeezed configuration with a contamination to the primordial non-Gaussianity of the local type 

of i/£°ri^o.7. 

Let us conclude this section by mentioning that the non-linearities emerging from secondary effects 
can be very interesting by themeselves, since they carry a lot of information about the evolution of the 
universe and the growth of structures after recombination till very low redshift. Detailed examples and 
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reviews about this aspect can be found, e.g., in Refs. [75" 1 176 ] 177]. 



9 How to analytically estimate the signal-to-noise ratio 
from the NG at recombination and its contamination to 
the primordial NG 

As we have seen in Sec. [7J the dynamics at recombination is quite involved because all the non-linearities 
in the evolution of the baryon-photon fluid at recombination and the ones coming from general relativity 
should be accounted for. Such a contribution is so relevant because it represents a major part of the 
second-order radiation transfer function which must be determined in order to have a complete control 
of both the primordial and non-primordial part of NG in the CMB anisotropics. The NG generated 
at the surface of last scattering comprises various effects. In this Section, we devote our attention to 
one particular relevant contribution, the one coming from the non-linear evolution of the second-order 
gravitational potential which grows in time on small scales. The analysis of this contribution offers the 
opportunity to understand how some secondary effects can contaminate mostly a given type of primordial 
NG, while leaving almost untouched other different forms of primordial NG. In this case this effect is 
a causal one, developing on small scales, so we expect that the NG it generates will be mainly of the 
equilateral type, rather than of local type. Therefore, a reasonable question is to which extent the NG 
from recombination alters the possibile detection of the primordial NG of the equilateral type. The goal 
of this Section is therefore to illustrate how to estimate in a semi-analytical way the contribution to NG 
from recombination. 



9.1 Signal-to-Noise ratio for the primordial equilateral bispectrum 

In this subsection we wish to recover the estimate for the signal-to- noise ratio (S/N) given in Ref. [36] for 
the primordial bispectra of "equilateral" type [75] by adopting a simple model. In other words, we test 
the goodness of the semi-analytical model we will be using in the next Section to estimate the bispectrum 
from the recombination era. 

Our starting point is the primordial equilateral bispectrum |78] 



($(k 1 )$(k 2 )$(k 3 )) = (2vr) 3 ^ 3 )(k 1 +k 2 + k 3 )B cquil (A :i ,A ;2 ,A : 3), (358) 

where 

e equil RA 2 



B equ il(k 1 ,k2,k S ) - f NL ■ 6A ■ ( 33 33 33 ,21,21.2 

\ i\' y i<- 2 \ 3 2 3 1 2 3 

+ 7TT2T3 + (5 perm.)) , (359) 

K/<2 J 

and the permutations act only on the last term in parentheses. The parameter J^l"' quantifies the level 
of NG while A = 17 AQ x 10~ 9 is the amplitude of the primordial gravitional potential power spectrum 
computed at first-order 
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(k 2 )) = (27r) 3 5< 3 >(ki + k 2 )P(h) , (360) 

with P(k) = A/k 3 . Since the signal-to- noise ratio (S/N) will be some function of the maximum multipole 
a given experiment can reach, ^ max 3> 1, we can use the flat-sky approximation [791 [80] and write for 
the bispectrum 

(a(4)a(f 2 )a(4)) = (2tt) 2 5^ (£i 2 3)B(£ 1 ,e 2 ,e 3 ) , (361) 
where £ 123 = h + h + 4 with [H] 

^2 
(2vr) 5 " 

x A T (£ 1 ,k z 1 )A T (e 2 ,k z 2 )A T (e 3 ,kl), (362) 
where k' means k evaluated such that P = £/(rjo — rj r ) and 



B cqui i(£i,£ 2 ,£ 3 ) = { \ ^ [dkfdkldklS^Xki^B^iki,^,^ 



A T (£,k z ) = I" j s^ik^+iy^-^,^^ , (363) 



o (vo ~ v) 

is the radiation transfer function defined by the CMB source function S(k,rj). In this notation, 770 and 
r) r represent the present-day and the recombination conformal time, respectively and k z and fc" are the 
momentum components perpendicular and parallel respectively to the plane orthogonal to the line-of- 
sight. 

The (S/N) ratio in the flat-sky formalism is }79|[80] 

I) 2 - / w^>&3) 6 ^;;^ 3) • (364 » 

where / s k y stands for the portion of the observed sky. In order to compute the bispectrum B equ n(£i,£ 2 ,£ 3 ) 
and the power spectrum C(£) we adopt the following model 

/dk z 
—e ikZ ^ - Vr ^(k')A T (£,k z ) (365) 
2ir 

where we mimic the effects of the transfer function on small scales as 

A T (£, k z ) = a( m - 1]r )-2 e -i/Wt*y\-i/m*\/^ t (366) 

i.e. a simple exponential and a normalization coefficent a to be determined to match the amplitude of the 
angular power spectrum at the characteristic scale £ ~ £* = k*(rjo — 77,. )H It is important to make clear 
what are the reasons underlying the choice of such a model. When computing the (S/N), Eq. (|365j) with 
£* = k* (770 — f]r) — 750 and a ~ 3 is able to account for the combined effects of "radiation driving" , which 
occours at £ > £ eq — 160 and boosts the angular power spectrum with respect to the Sachs- Wolfe plateau, 
and the effects of Silk damping which tend to suppress the CMB anisotropics for scales £ > £jj ~ 1300. 



6 We could equally choose a transfer function as A T (£, k z ) = a (770 — T) r ) 2 e 1 / 2 ^/^*) 6(k* — \k z \), the relevant 
approximation being that the integral over k z is cut at the scale fc* . 
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The combination of these effects produces a decrease in the angular power spectrum from a scale £* ~ 
7500 The power spectrum in the flat-sky approximation is given by (a(li)a(l2)) = (2ir) 2 5^ (lvi)C(li) 
with 

C (£) = {7]0 {2 J ) r)2 J dk z \A T (£, k z )\ 2 P(k) . (367) 

The exponential of the transfer function for Eq. (|365|) allows to cut off the integral for k — fc* and one 
finds (see also Ref. [80] ) 

where the last equality holds for I 3> £*■ To compute the bispectrum we proceed in a similar way. One 
first uses the Dirac deltas, S^(kf 23 ) and <S (2) (^i23)- Then it proves to be useful the change of variable k\ = 
xi£\/(rjQ — rj r ), k\ = x 2 £2/(vo — Vr)- In this way the transfer functions become A T (^j, &f) oc e - 1 / 2 (\ x i\ e i/ 1 *) 
which allows to cut the integrals over Xi (i = 1,2) at £*/£i. Now, as a good approximation to see the 
effects of the transfer functions, we can take I* and thus the integral over Xi can be easily computed 
by just evaluating the integrand in x% = times \{t* j j t-i) . With this approximation the integral in 
kf is easily obtained and we get for the bispectrum 

(2iry 



x 



— — — _L ^ 

/3/-3 /3*3 /3/3 p2p2p2 + /,/2/3 + P erm -J ) ) 
t 1 t 2 t 1 t 3 C 2 C 3 ^1*2^3 £ 1>2 £ 3 / 



(369) 



where 



£^ = £{+t 2 + 2£ l - £ 2 . (370) 

The coefficient f\ ~ 1/1.4 = 0.7 is a fudge factor that improves the matching between our approximation 
for the bispectrum and numerical results that have been consistenly checked. Notice that, according to 
our approximation, the equilateral structure of Eq. (|359p is preserved in i spaceJl £\ = £2, in correspon- 
dence of the maximum following expression In computing the signal-to-noise ratio, consistency with our 
approximation (|369j) requires that we integrate over £\ , £2 starting from a minimum £ m \ n > £* up to £ max 
and paying attention to the fact that even £3 in Eq. (|370p must be larger than £ m in- The scaling with 
4nax with respect to the case of a local type bispectrum turns out to be much milder. While for the 
local type (S/N) 2 oc ^ax [HO], for the equilateral bispectrum (|358|) we fincH (S/N) 2 oc £ max and, setting 



7 The choice of the exponent 1.2 derives from the study of the diffusion damping envelope in Ref. [81] . 

8 The expression (|369j) can be also written as 
B cqu n(£i,£ 2 ,£ 3 ) = (2^)- 2 48/ 1 /° q L uil a 3 A 2 ^ e -(^- 2 +4- 2 +4- 2 )/2^- 2 (1 + cos 9)(h + i 2 - 4)/^|4 6 being the angle 
between £\ and £ 2 . 

9 These scalings can be easily understood by analyzing the expressions (S/N) 2 for the local and equilateral 
primordial NG. In the local case, (S/N) 2 is proportional to / d 2 £ 1 d 2 £ 2 d 2 h8^\li 2 ?,)(£l+ll+£l) 2 / (^^f ED]; since 
the squeezed configuration, e.g. l\ <C £2, £3, is dominating the local bispecrum, the integral becomes proportional 
to / d£ x d£ 2 (£ 2 /£\) oc In the equilateral case, however, (S/N) 2 receives contributions from the configuration 

which is peaked at £\ ~ £ 2 ~ £3 and therefore it can be written as J d 2 £\d 2 £ 2 8(£\ — £ 2 )/£\ oc £ max . 



N2 



4 = 750 and £ min ~ 1200, 

(371) 

By choosing / s k y = 0.8 and £ max = 2000 we find a minimum detectable 

/£i uil ^ 66, (372) 

obtained imposing (S/N) cqu n = 1. Both the estimate of the minimum value of /nl"' an d the scaling 
(S/N) 2 oc £ max are in remarkable agreement with the result obtained in Ref. [36] where the full transfer 
function is used and a value of J^l" 1 = 67 is obtained Notice that our estimate is independent from 
the coefficient a and the exponetial e -1 / 2 ^^*-* 1 ' 2 introduced below Eq. (|365p to mimic the full transfer 
function. This is because there is an equal number of transfer functions in the numerator and denominator 
of the expression (|364p for the signal-to-noise ratio and their effect tend to cancel despite they are not 
simple multiplicative factors (see discussion in Ref. [80] ). 



9.2 Non-Gaussianity from recombination 

Comforted by the goodness of our model, in this Section we wish to estimate the level of NG generated at 
the recombination era. One can check that on small scales the second-order anisotropies are dominated 
by the second-order gravitational potential $( 2 ) which grows as ry 2 , as first pointed out in Ref. [82]. Let 
us first briefly show how this result can be obtained using the calculations of Sec. Consider the tight 
coupling regime, and small scales, well inside the diffusion Silk-damping scale Xd- Then look at Eq. ([229D . 
One realizes that, apart from the combination Aqq + 4<3?( 2 \ all the other terms scale at most as 7] or 
they are suppressed either because in the tight coupling regime or because of the damping diffusion on 
small scales. On the other hand one finds that the main contribution to the bispectrum generated at 
recombination comes from 

e& = iAg+*« (373) 
which is the usual term appearing in the CMB anisotropies due to the intrinsic photon energy density 

(2) 

fluctuations Aq and the gravitational redshift due to the potential. In fact such a term on large scales 
reduces to the Sachs- Wolfe effect, while on small scales at recombination, using Eq. (I279P and (1282j) . we 
find 



= -^G(k 1 ,k 2 ,k)T(A: 1 )<I>( 1 )(k 1 )T(A: 2 )cl»( 1 )(k 2 ), 



4 oo 



(374) 



10 We thank M. Liguori for discussions about the minimum value of detectable by Planck and for its scaling 

With £max. 
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where we have evaluated the expression at the recombination time r/ r and R = 3/9b/4/3 7 is the baryon-to- 
photon energy density ratio. In writing Eq. (13740 we have dropped off the contribution from primordial 
non-Gaussianity (which, anyway, always propagates linerarly) since here we are focusing on the relevant 
secondary effects that can constitute a source of contamintation to it (see Eq. (ll99p and the discussion 
below). Notice also that this expression has been obtained assuming all the momenta much larger than 
k e q. 

Eq. (|374p has been confirmed numerically in Ref. [82]. There is also another way to obtain the 
analytical expression in the first line of Eq. (|374p , as discussed in details in Ref. [82] , which exploits some 
well-known results at linear order extending them at second-order. At second-order in the perturbations 
the solution for the acoustic oscillations will have a form very similar to the linear solution written in 
Eq. (IC.18p . except, as usual, of some source terms S made by products of fistr-order perturbations. 



A 



(2) 
00 



A 



(2) 
00 



(0) + (1 + R)& 2 ^ cos(kr s ) - (1 + fl)$( 2 ) + S . (375) 

Now on sufficiently small scales the products of first-order terms are indeed suppressed, due either to 
Silk damping, or because gravitational potentials are suppressed as (k/k equ ) 2 . Recall that only acoustic 
oscillations are damped by Silk damping diffusion, so that the cosine in Eq. (|375p is multiplied by e~^ k / k °^ 
k S> kn, where k~^ = indicates damping length, but the baryon induced shift (1 + i?)^ 1 ) is left 
untouched. Therefore on small scales the combination of the damping effects and the growth of the second- 
order gravitational potential as rj 2 single out the dominant effect in Eq. (l375p to be — —R& 2 ^. 

9.2.1 Bispectrum from recombination due to the non-linear growth of the gravita- 
tional potential 

Let us start to go into some details of the form of the bispectrum induced by the contribution (I374p . In 
Eq. (I374p the kernel G is given by 

1, i,\ 1, 1, 10 (k-k!)(k-k 2 ) 
G(ki,k 2 ,k) = ki -k 2 - y p . (376) 

From the form of the kernel we see that the NG at recombination is dominated by an equilateral con- 
figuration, as expected from the fact that its origin is gravitational. Here and in the following we are 
implicitly assuming that a convolution is acting on the kernel as 

-^3 J dk 1 dk 2 j( 3 )(k 1 +k 2 + k 3 )G'(k 1 ,k 2 ,k)T(A ;i )cI»( 1 )(k 1 )r(fc 2 )ci>W(k 2 ). 

The reader should remember that, at first-order in perturbation theory, the combination + R&^ 
is exponentially suppressed by the Silk damping, but still greater than the term R^ 1 ^ (which does not 
suffer the damping) for the maximum multipole of interest, £ max ~ 2000. This is meanly due to the fact 
that the first-order gravitational potential rapidly decays on small scales. On the contrary, at second- 
order in perturbation theory, the gravitational potential grows like the scale factor on small scales and it 
turns out that the i?<E>( 2 ) dominates on small scales (see Ref. |82|). 

The gravitational potential at linear order can be expressed as usual in terms of the transfer function 
T(k) 

T(k) ^ 12 Mfc/Sfceq], (377) 
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where the last step is an approximation valid on scales smaller the the equivalence scale, k 3> k cq . In the 
following we will account for the logarithmic growth just with a coeffcient To(k) = 121n[fc/8fc e q] ~ 11 for 
the scales of interest. 

In the flat-sky approximation one arrives at an expression similar to (|362h . where now one of the 
linear transfer functions must replaced by a transfer function at second-order. Specifically one finds 

0?o - Vr) 2 



Brec(£ 



(2*) s 



dkfdk z 2 dk z 3 5^(k z 123 ) G(UM,U 3 )T(k{)T(k> 2 )P(k{)P(k> 2 ) 



x A T (£i,A;f)A T (4,A ; |)A T ( 2 )(4,A ; |) + cycl 
By using our model (|365h and 



1C 



A T ( 2 \e,k z ) 



R 



14 (770 - Vr) 2 

for the second-order radiation transfer function, we find 



(378) 
(379) 



B iec {£i,e 2 J 3 ) = - 

x f dkfdk^dk^ (fcf 23 



R (VP - Vr) 4 

14 (2tt) 2 '" c 



G(ki,k' 2 ,k^) 5 5 e 



-l^dfc^l/fc,) 1 



-l^dfca,!/*.) 1 



h'°h' 



+ cyclic 



(380) 



At this point we proceed further by employing the same approximation described after Eq. (|368p . We 
use the Dirac delta to replace the variable k 3z , and the exponential allow us to evaluate the integral for 
ku = k 2z = 0, for scales t{ S> £*. This leads to 



B Tec (hA2,£ 3 ) 



4 /2 R A2 2rr2 



(2vr) 2 14' 



10 (4-4)(4-^ 2 ) 



+ cyclic . 



(381) 



Again here ji is a coefficient to better calibrate our approximations with numerical results that we have 
performed in order to test the validity of our approach. Not surprisingly, it turns out that jfe — j\ — 1/1.4. 



9.3 Contamination to primordial non-Gaussianity from recombina- 
tion: Fisher matrices 

Our goal now is to quantify the level of NG coming from the recombination era and to estimate the 
level of degradation it causes on the possible measurement of the equilateral primordial bispectrum. The 
reader should keep in mind that, given the form of the kernel function (|376p . the NG from recombination 
is expected to be of the equilateral type. Following Sec. 18.11 a rigorous procedure is to define the Fisher 
matrix (in flat-sky approximation) as 
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where i (or j)= (rec, equil), and to define the signal-to-noise ratio for a component i, (S/N)i = 1/y F~ , 
and the degradation parameter di = FuF^ 1 due to the correlation bewteen the different components 
rij = F~- 1 / J F^ l F~^ . The first entry -F C quil,cquil of the Fisher matrix corresponds to the (S/N) 2 ratio 
computed in Eq. (|37ip which does not account for any kind of cross-correlation. Due to the equilateral 
form of the NG generated at recombination we expect that the minimum value detectable for f^ 11 will 
be higher that the one reported in Eq. (|372p . For the mixed entry we find 



rcc, equil 



/sky 1 



7T (2tt) 



d 2 £ 1 d 2 £ 2 d 2 £ 3 8^{£ 123 ) 



iW4,4,4)£ equi i(4,^,4) 



6C(4)C(4)C(4) 



-3/i/ : 



/ sky 4R 48 T 2 



vr 3 14 2 5 6 a 
d£ 1 d£ 2 (l + £i-£ 2 /ht 



(k CCL r] r ) 2 £ 2 £* j4/nl 



l^s//.) 1 " 2 



1 



(4 + ^ 2 - 4) 



'1 • t2 



io (4^i)(4-4) 



(383) 



where £3 is given by Eq. (|370p . The factor 3 in front of this expression comes from cyclic permutations. 
The integral can be performed numerically and, integrating from a minimum £ m { n ~ 1200 up to £ m3iX = 
2000, and by taking R ~ 0.3 when evaluated at recombination, a ~ 3, Tq ~ 11, (k eq r] r ) 2 ~ 26, £ eq = 150, 



4 = 750, we find F, 



roc, equil 



9.4 x 10 . Finally for the entry -F re c,rec we get 



/sky 1 /" ,2« j2« A 1n s(2)(p \ B 2 cc (£i , £ 2 , £3) 



f: 



2 /sky /^Wljf 
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x Ui • 



1 • ^3 
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(384) 



0.014. We are now able to compute the entries of inverse of the Fisher 
matrix, F^ . In the following we report our results for the signal-to- noise ratios and the degradation 



and W6 find £1 Valll6 -^rec,rec 
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12.6 x 10- 3 /N q L U , 
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(385) 
(386) 
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F" 1 



r rec , equil = :~-0.53, (387) 

/ 1?— 1 z^— 1 



'^equil,equil^ lrec i rec 

dree — Frec,recF rec rec — 1.4 , (388) 
dequil = F equ ii iCqu iiF cquil cquil ~ 1.4 . (389) 

As a confirmation of our expectations, we find that the NG of the type given by Eq. (|376|) has a quite 
high correlation with an equilateral primordial bispectrum. This translates into a degradation in the 
mimimum detectable value for J^l" 1 with respect to the value given in (|372p . In fact from the signal-to- 
noise ratio (|385p we find a minimum detectable value of 

/Nl Uil -79, (390) 

imposing that (S/N) cqn n = 1 with an an increase of Af^jf = O(10). This corresponds to an increase of 
the l-o- uncertainty on /^ uil of 20% (see Sec. KTbN 

Following Sec. 18. 1( in order to measure the contamination to the primordial bispectra one can define 
that effective non-linearity parameter which minimizes the x 2 defined as 

2 



X 1 = I d z £ 1 d z £ 2 d 2 £ 3 5 (2) (£ 12 3) - 



6C(h)C(£ 2 )C(£ 3 ) 



(391) 



to find 

fcon _ -^rec.equil 



^con _ ^ rcc, cqmi 

/NL — p f equil_ 1 ' [6\)Z) 

" cquil, cquil Jnl — 

and an analogous expression to compute the contamination to the local primordial bispectrum. In this 
case we find an effective non-linearity parameter f^{f — 5@ 

Similarly we can compute the Fisher matrix accounting for the NG generated at recombination and 
the primordial NG of the local type 

<$(ki)$(k 2 )$(k3)> = (27r) 3 5( 3 )(k 1 +k 2 + k 3 ) J Bi oc (A: 1 ,A :2 ,A ; 3), (393) 

where 

B loc (h, h, h) = • 2A 2 • (-I3 + -I3 + * ) . (394) 

\ t\j-t Fury t\j-i t\jn r\jij rvn J 



11 Due to a non- vanishing correlation, ry, (S/N) gets modified from its zero-order value to (S/N) = (S/N)q(1 — 
rfj) 1 / 2 , so that the minimum detectable value of gets shifted by a quantity A/^"" 1 / (f^ a )o = (1— Ty) — 1 / 2 — 1. 

12 Notice that sometimes in the literature one can find also an effective non-linearity parameter /^f; defined in 
such a way that the equilateral (or the local) bispectrum has the same Fisher matrix errors as the recombination 
bispectrum 



/•off. V^" 
/NL 



jequil 1 

Jnl — 1 



quil, cquil 1 

However this is not the proper quantity to define the contamination level to primordial NG; here we are just 
comparing signal-to-noise ratios, while the contamination /^™ defined in Eq. (|400[) contains a somewhat richer 
information, since we are asking what is the value of equilateral (local) /nl which best mimics the bispectrum from 
recombination. 
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The bispectrum and the signal-to-noise ratio as defined in Eq. (|364p have already been computed in 
the flat-sky approximation in Ref. [80]. The result is that (S/N)f oc = 47r~ 2 /sky(^*/^mm)(/N°£) 2 -4^m a x 5 
corresponding to a minimum detectable value of = 0(7) for £ mSLK = 2000 (when other possible sources 
of NG are ignored). We can compute the off-diagonal entry of the Fisher matrix in a similar way to what 
we have described in this section, and we get F reCjloc ~ 8 x 10- 3 /nl- Finally the entry F reCjrec ~ 0.014 
has already been computed above. From inverting the Fisher matrix, we get the following signal-to-noise 
ratios and the degradation parameters 

S \ _ . 2 j^loc 



loc ^ F loc,loc 



N J ro „ l F -i 



= ^14xHTX L c , (395) 
0.1, (396) 



rcc.crcc 



p- 1 

r rcCjloc = _ i reC ''° C _ 1 ^-0-44, (397) 

V ^loc,loc-^ rec > rec 

dree — -^rec,rec-^rec,rec — ' (398) 
^loc = ^equil.equil-P'equi^equil — • (399) 

In particular, from Eq. (|395p we see that now the minimum detectable value of /^l remains basically 
unchanged in the presence of the recombination signal. Similarly the effective /^if rea ds 

f eon = ^oc 3> (40Q) 

^loc,loc J NL 

which is much smaller than the effective non-linearity parameter (|400|) for the equilateral case. We 
have also checked the cross-correlation between the primordial local and equilateral bispectra finding 
a value of noc,equil — 0.23, which is in agreement with the value reported in Ref. [36] • This reflects 
the fact that the primordial local and equilateral signals are not fully uncorrelated. The reason is due 
to the fact that the equilateral and local bispectrum (|394p and (|359p approach the same shape in the 
equilateral configuration. This is also the reason why the cross-correlation between the primordial local 
and recombination bispectra is not so small. 

10 How to perform a numerical analysis of the CMB bis- 
pectrum produced by second-order perturbations 

In the previous section, through a specific example, we have learned some basic principles to determine 
a reasonable and quite fast analytic estimate for the contamination to the primordial NG. In this section 
we provide all the tools necessary to a full numerical implementation of the second-order Boltzmann 
equations which allow to obtain a systematic numerical evaluation of the CMB angular bispectrum 
produced by second-order cosmological perturbations. This section is mainly based on the results of 
Ref. [85. 

In particular we will apply this formalism to numerically evaluate the contamination to primordial NG 
from the second-order fluctuations of the Boltzmann equations that come as products of the first-order 
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perturbations, and ignore the intrinsically second-order terms, or the effects of the perturbed recombina- 
tion [721 EE G3]. El Here therefore we come to the discussion of the second example mentioned in Sec. 17.11 
As anticipated there we will see that, unlike the intrisically second-order contribution considered in the 
previous section, the products of the first-order perturbations give a CMB bispectrum that peaks in the 
squeezed configuration. Therefore we shall study the contamination of the primordial NG of the local 
type. 



10.1 CMB Bispectrum from second-order perturbations 
10.1.1 Definitions 

Again, we expand the temperature fluctuation into the linear (first-order) part and the second-order part 

as 

AT(n) AT«(n) AT< 2 )(n) /A . 

—fr 1 = + jr U - + ---- ( 401 ) 

The spherical harmonic coefficients of temperature anisotropy, 
&tm = T j (i 2 fiy^ l (n)AT(n), are therefore expanded as 

a /m = a£>+a£ + .... (402) 

Recall that to expand the Boltzmann equation up to the second order in perturbations, we had to expand 
the distribution function, 



f(x,p,h,i]) = 2 



(403) 



up to the second order in perturbations: 6 = 6^ + /2 + . . . , and accordingly / = /(°) + /W + 
/( 2 )/2 + . . .. We then had to compute the fractional perturbation in photon's energy density at the i-th 
order in perturbations, AW, by multiplying /W by p, and integrating over p 2 dp: 

A (0 = f d PP 3 f il) m) 

At the linear order, we recovered the usual relation between the linear fractional temperature fluctuation, 
O^ 1 ) = AT^ 1 ) /T, and the linear fractional energy density perturbation, A^ 1 ) = 5p^ /p 7 , i.e., A^ 1 ) = 40W. 
At the second order we have 

A (2) = 4@ (2) + 16 [ (1)]2 ) ( 405 ) 

which is related to the second-order temperature fluctuation as 

^ (2) = U^-(A^))- 3 -([e w ] 2 -([e {i m 



I (@m _ (e( 2 )) + [e«] 2 - ([Q w ] 2 )) , (406) 



13 At the moment of writing this review the calculations that also include the intrinsically second-order terms 
and the perturbed recombination are in progress and will be presented in Ref. [85) . 
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where we have subtracted the average of the temperature fluctuation so that the average of AT^/T 
vanishes. 

We compute af^ from AT^/T using 

^ - ^nY^n)^ (407) 
(2) 3 V- V- / -i\mr>— mm' m" („iX) / „(1) „iX) 



f'm' i"m" 



where we define 



/ , / j \ 1 ) y ee'£" \ a l'm' a i"m" \ a l'm' a l"m"l /' 
I'm 1 t"m" 

a S = g / d 2 ^; m (n)(A( 2 )(n) - <A( 2 )(n)», (408) 
TO"* = / d 2 nFWfi)yWn)>Wn) 



'(2^ + 1) (2^2 + 1) (2^3 + 1) I h h h\( h 



1-2 



V 47r I J \ mi ni2 m% J 

Here the matrix is the Wigner 3j symbol. 

The CMB angular- averaged bispectrum, B^^^, is related to the ensemble average of a£ imi a£ 2m , 2 a^ 3m3 

as 

B txhh - ( m \ m 2 m 3 ) ^i m i a ^m2 a ^m 3 )- ( 409 ) 

all m \ / 

This definition guarantees rotational invariance for the bispectrum, and the Wigner 3j symbol ensures 
that the bispectrum must satisfy triangle conditions: £i — £j\ < £k < h + £j for all permutations of indices, 
and selection rules: mi + mi + 1713 = 0. 
The ensemble average is given by 

/ \ / (!) C 1 ) ( 2 ) \ 1 7- 

(a fimi a £2m2 a f3m3 ) = (a^o^a^) + cyclic = 

(!) „(!) ^( 2 ) \ 3 t ■\\m 3n - rn3m 'z m 3 



\ a i 1 m 1 a l 2 m 2 a l i m i l r, 2^ V X J ^ 3 £'£" 



4 m 3 « 



x ( 4i L 2 a ^L 3 ) - 4il 2 > <4 3 L 3 4>m>> > ) + c y clic > 



(410) 



where cyclic means that we have to sum the cyclic permutations of Eq. (I410j) for indices (1,2,3) — > 
(3,1,2) -> (2,3,1). 

As we assume that o^'s are Gaussian random variables, the four-point function of in Eq. ()410p 
is given by the sum of products of all possible pairs. Each pair gives the angular power spectrum, Cf. 

($141) = (-l) m C e 5 u ,5- mm ,. (411) 
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We obtain 

(a {1) a (1) a {1) a (1) ) - (a {1) a (1) Va (1) a (1) ) 

\ a hmi a t 2 m 2 a e' 3 m' 3 a l%m' 3 ' ' \ a e imi a £ 2 m 2 1 \ a £' 3 m' 3 a e%m' 3 ' ' 

= (-ir i+Tn3 QiCi 2 [^5_ mi ^^<y_ ma ^ + (1 O 2)]. (412) 

Substituting the right hand side of equation (|412p for the second term of equation (|410p . and using 
4 + 4 + 4 = even, we obtain the angular averaged bispectrum, 

Bhi 2 i 3 = Btxi 2 h - ^h^i^C^C^ + cyclic), (413) 

where we have defined the quantities, 



'(24 + 1) (24 + 1) (24 + 1) ( h 



■2 



4tt V 



(414) 



and 




= E ( * Z X ) (4L a ZAs) + ^ lic - ( 415 ) 




10.1.2 The second-order CMB radiation transfer function 

In this section we show how it is possible to define in a rigorous way the radiation transfer function 
for CMB anisotropics at second-order in the perturbations. It is a generalization of the well known 
quantity used at linear-order and as such it allows us to develop a systematic numerical analysis of the 
angular bispectrum produced by second-order perturbations in the very same way as at linear-order 
various numerical codes are nowadays available for the computation of the CMB power spectrum, such 
as CMB FAST or CAMB. 

The starting point is the Boltzmann equation that governs the evolution of A^ 1 ) (k, fi, rj) and A^ 2 -* (k, n, rj), 
where fi = k ■ n and n is the direction of propagation of photons. Note that for the linear perturbation 
there is azimuthal symmetry such that A^ 1 ) depends only on the angle between k and n; however, for 
the second-order perturbation there is no such symmetry. We write again the Boltzmann equations in 
Fourier space 

AW' + i^AW-T'AW = S^{k,n, v ), (416) 
A( 2 )' + ikfiA^ - t'A® = S (2) (k> fi) 7?)) (417) 

where the primes denote derivatives with respect to the conformal time d/drj, and are the source 
functions at the first and the second orders, respectively, and r' is the differential optical depth which is 
defined by using the mean electron number density, n e , the Thomson scattering cross-section, or, and 
the scale factor, a, as 

t' = -n e a T a. (418) 
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We expand again the angular dependence of A" as 



and that of the source terms as 



'21 + 1 

Air 



'21+1 

47T 



d 2 hY; m (h)S^(k,h,rj), 



(419) 



(420) 



where i = 1,2. 

The source functions relate the observed a&n's to the primordial curvature perturbations in comoving 
gauge, C(k)- The relations contain the linear radiation transfer function, gi{k), and the second-order 
radiation transfer function, Ff™ (k), and are given by 



a 



(i) 



d 6 k 



a 



fan 



(2vr)3- 
(fk 
(2i 



,(k)Y; m (k)c(k), 

d 3 k' 



x Y: Ff^'iy, k", k)y; m ,(k)c(k')c(k 

I'm' 



(421) 



(422) 



The linear transfer function is given by 

g t (k) = J^d^^\k, V )+S^{k, V )^ + Sg\k,r,) 



1' 

2 du 2 + 2 



where u = k(r]Q — rj) and sj^ is the standard linear source function 

sg\k, v ) = 4*( i ) , (m)-t , a; ) i) (m), 

S[ l j(k, v ) = 4fe$«(A;,r ? )-4T / 4 1) (A ; ,r ? ), 
sg\k,r,) = yA«(M), 



(423) 



(424) 
(425) 

(426) 



where ^^(A;, r/) and ^(^(/c, 77) are the metric perturbations at the linear order in the longitudinal gauge 
and A^\k,r]), Aj (k,r)), and A^ 1 (k,rj) are the coefficients of the expansion in Legendre polynomials of 
AW(/c,/j,7?). These coefficients A.f\k,r)) are related to A^ (Eq. ([415]) ) via 



47T 



r y/ m (k)A«(k)(2£ + i) 



2^ + 



(427) 



The first-order velocity perturbation, Vq (k,rj), is the irrotational part of the baryon velocity defined by 
v(k) = — ivo(k)\i. 
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The new piece, the second-order radiation transfer function, is the line-of-sight integral of the second- 
order source terms in the Boltzmann equation: 



F/;-'(k',k",k) = ^£(-in-;) A -^ir' M i 



4tt 



A// 

'"no 



2A + 1 



x jf ° d^sg (k' ; k", k, r?)^ [fcfa - %)] . 

(428) 

Here, we have introduced a new function, 5L:(k , k , k, 77), which is denned by the following equation: 

sj2(M) = /^dk'V 3 )(k' + k''-k)5S(k / ,k'',k, ?? )c(kOC(k / '). 

(429) 

/n't 

Basically, S\ m '(k', k", k, 77) is the second-order source function divided by 
C(k')C(k"). The explicit expression of the second-order source function can be be read from Eq. (|134p . 
Using equation (I42ip and (I422j) . we calculate the first term in Eq. (I413p . Bi ± i 2 i 3 , as follows: 

\ ' £ 3 M 3 i 

x y; 2m2 (k 2 )y/ 3M3 (k 3 )^ 1 (A ;i )^ 2 (fc 2 )P c (A :i )P c (A ;2 ) 

x {F£^(k l5 k 2 , k 3 ) + F^(k 2 , kx, k 3 )} , 

(430) 

where, as in the previous section we will sometimes use the notation ki 2 3 = ki + k 2 + k3. P^(k) is the 
power spectrum of £ given by the usual definition: 

<C(ki)> =0, (C(ki)C(k 2 )> = (27r) 3 5( 3 )(k 1 + k 2 )P c (A: 1 ). (431) 

In order to perform the integral over angles, k, we expand the three-dimensional 5-function using 
Rayleigh's formula, 



^)(k 1 + k 2 + k 3 ) = 8 £ ^-^^^y^cti)^^) 

allZ'm' 

x ^ m ^(k 3 ) y drr 2 j^(rfa)je 2 (rfa)j^(rfa) , 

(432) 

and also expand the angular dependence of <S^(ki, k 2 , k 3 , rf) by introducing the transformed source 
function, S^jg (fa , k 2 , k 3 , rf) , as 

4-7T / 47T 



5g ) w ,(k 1) k a) ka,t / ) = £ EH^ Al+A2 

Ai,£ti A2,^2 
:MiA t 2At3 



V 2Ai + 1 V 2A 2 + 1 

x s^(fa, fa, fc 3 ,r ? )y AlAtl (k 1 )y A2M2 (k 2 ) . 



(433) 
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This result shows that <S\ 3 ^ 3 

(ki,k 2 ,k 3 ,?7) = 

^3/^3(^1' ^2 j &3) r /)> an( i thus i 7 /^™ (ki,k2,A;3) follows (see 

Eq. (gZBD )- 

Now we can perform the angular integration of Eq. (|430p to obtain 

all m all £'m' all X/i 



4ir 



V (2Ai + 1)(2A 2 + 1)(2A 3 + 1) V mi m 2 m 3 J 

xy e' 1 e , 2 £' 3 y £' 1 e 1 M y e 2 e 2 x 2 y e 3 i 3 x 3 

3 

i=i 

xi*+4 I d^{S^(k 1 M,h,v) + l^2}j^[k 3 (r,- Vo )] 

+cyclic, (434) 



where for brevity we indicate -P\iA 2 a 3 = \/47r/[(2Ai + 1)(2A 2 + 1)(2A3 + 1)], and we have used the fol- 
lowing relation of the Wigner 9j symbol, 



1)444+4 yw* 1 t2 t3 )^ m ^;rA miM1 (435) 

z — ' . \ ?7li 777,2 m.3 / '1*2*3 Vi-M 



all mm 

xy £'£ 2 x 2 y tu 3 x 3 ~ ( -~ L) h'^hxt'^h^Mh^Xi 



C £ 2 ^3 

x { *-\ ^2 ^3 
^ Ai A2 A3 



Ai A2 A3 

Ml M2 M3 



with R = £± + £2 +£3 + ^"1 + ^2 +^3 + A1 + A2 + A3. The Wigner 9j symbols have the permutation symmetry: 

f 4 *2 4 1 

(-1) B 4 ^2 ^3 = 
Ai A2 A3 ^ 



«2 


4 


4 1 




f 4 


£3 


<2 








hi 




l 3 




A 2 


Ai 


A3 J 






A3 


A 2 




< 2 








^2 


«3 








hi 




A 2 


A3 


Ai 


A 2 


A3 J 




I £[ 


' 2 


f 



(436) 

and the coefficients I? p p , Ie 2 e' \ 2 , and If g* \ , ensure ^ + £ 2 + ^3 = even, £1 + £[ + Ai = even, 

12 3 12 3 

£2 + £2 + A2 = even, and £3 + £ 3 + A3 = even, respectively, which gives R = even. Hence the Wigner 9j 
coefficients are invariant under the permutations. 

Finally, we obtain the angular averaged bispectrum, 



4 / 47T 



j 3 -e'+n 



n \tPx\ (2Ai + 1)(2A 2 + 1)(2A 3 + 1) 
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* Ai A2 A3 > 



i=l 



x J dr'e- T ^j e3 {r'k 3 )S XlX2X3 (k 1 ,k 2 ,k 3 ,r / ) + perm., (437) 

where r' = (rjo — rf) and we have used the relation of the spherical Bessel function, ji(—x) = (— 1) jt(x). 
We also define the "angular-averaged source function," 

Sx 1 x 2 x 3 (h,k 2 ,h,r) = W^ 1 i 2 ^^KffiZfa'bM 



„. Al+A / 2A 1 + 1 / 2A 2 + 1 ^/ Ai A 2 A 3 
V 4vr V ^ a% V m ^ 2 ^ 3 



(438) 



where we have used the inverse relation of Eq. (|433p . Note that cyclic terms in Eq. (|434p have become 
permutations because of invariance of the Wigner 9j coefficients under the permutations. 

The final analytic formula (|437p we have obtained is a general formula which can be applied to any 
second-order perturbations. The information about the specific second-order terms is contained in the 
angular-averaged source term, 5aiA 2 a 3 (see Eqs. (|438p and (|433p for the definition). 

For terms that are products of the first-order perturbations, there is indeed a further simplifi- 
cation. One can show that 5a 1 a 2 a 3 (ki, &2> k 3 , rj) does not depend on k 3 , i.e., S\ 1 \ 2 \ 3 (ki,k2,k 3 ,7]) = 
<Sx 1 x 2 x i (ki, k<ii rf). This property enables us to integrate Eq. (|437h over k 3 . We obtain 



* IPx V ( 2A i + !)( 2A 2 + l)(2As + 1) 
f h h h \ 

^h^Ejh^xJiii'^hz^xA t'l ^2 ^3 \ 

* Ax A2 A3 ' 

x y dre~ r JJ J dkikiP^(ki)jir.(rki)g£ i (ki)Sx 1 X2X 3 {ki,k2,r) + perm. 



i=l 



(439) 



where r = (r/o — r/), i? = l\ + £2 + ^3 + ^"1 + ^"2 + ^3 + Ai + A 2 + A3, and we have used 



dk 3 k 2 3 j e , (rk 3 ) Je , {r'k 3 ) = -^«5(r - r'). (440) 



2r 
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Finally, by adding the remaining term in the full bispectrum, Eq. (|413p . we obtain 



d _ 2 y> / 4vr J 3 -i' 3 +R 

* £px V ( 2A i + !)( 2A 2 + 1)(2A 3 + I)' 



1 «2 £3 



Ai A2 A3 > 

dre~ T ^\ I dkik?P{(ki)j e >.(rki)g^(ki)S Xl \ 2 \ 3 (h,k2,r) 
i=i J 
3 

-^fi&fe^i^ +perm. (441) 

The remaining task is to calculate the angular-averaged source term, 
<5aiA 2 a 3 (&i 5 &3> t) according to Eq. (|438p . Up to now our scope has been to offer the reader all the 
general tools that allow a systematic treatment of the CMB angular bispectrum from second-order pertur- 
bations. In the following we will give just some examples of how such tools can be implemented in order 
to obtain numerical results about signal-to-noise ratios and contamination to primordial non-Gaussianity. 



10.1.3 Source Term: an example 

The explicit expression for Se'Qs., rf) can be obtained from Eq. (|134p . Here we do not want to report the 

(2) 

complete expression 5^(k, rj) of the source function at second-order, rather, for the goal of this Review, 
we think it is more instructive to show explicitly how the calculation of 5a 1 a 2 a 3 proceeds focusing on just 
one simple example. 

Therefore consider, for example, the term of the second-order source term from Eq. fjl34H given by 

2 r'^W x A« . (442) 

(2) 

First we compute the multipole coefficients of the source term S^(k) from this contribution, as 
defined in Eq. (|420p . They are given by the convolution 

d 3 k 



Now for A^(k) we use 



Afi(k) =zY^ T y; m (k)A«(k)(2£ + l), (444) 



so that 



x 2r',5 e (k 2 )A( 1) (k 1 )y; m (k 1 ). (445) 
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We now compute the corresponding "angular-averaged source function" coefficients S\ 1 \ 2 \. A (ki,k2,r) 
denned by Eq. (|438H . From Eq. (|445p you read the kernel defined in Eq. (|429|) 



S^M^) = i\l^- l (2£ + l)2r'5 e (k 2 )Aj 1) (k 1 )y/ m (k 1 ) . (446) 



From here 



3 



S i2) (kik,r) - ^. ^1 + 1 / 2^ + 1 ^ A t A 2 A 



d 2 h 2 Y^ (k)Y^ 2 (k 2 )5^ 3 (k l5 k 2 , r) 



V 4vr V 4vr V 2A 3 + f 



nV 2 (k2)^ 3 ( fc 2)2r'<5 e (A ;i )Ag(A ;2 ) 



I /i 2 — A*2 



2r'«5 e (A ;i )A«(fc 2 )^(-iy 



/'2 



(447) 



where we have used Yg* = (—1) m Yi_ m , and the orthonormality of the spherical harmonics. Using the 
property of the Wigner 3-j symbols 

D-ir(!J £ _^)=i;^ t = (-i)'-^tt. (« 8) 

we find 

^A 2 A 3 (^^2,r) = 2r' (2A 2 + 1)( 3 / 2 ) ^JA^fe, ^o^Aa • (449) 

This term therefore corresponds to coefficients 5oa 2 a 2 {ki, & 2 , r). 

In general the perturbation variables of the source term can be split into two parts (see, for example, 
Eq. (ll34p ). A part containing perturbations that are intrinsically second-order (these perturbations have 
superscripts (2), and uj m and Xm are also intrinsically second-order). Solving for these terms requires 
solving the full second-order Boltzmann equations coupled with the Einstein equations. 

Another part contains terms that are products of two linear variables, as the example that we have 
just considered. Evaluation of these terms is much easier than that of the intrinsically second-order 
terms, as the first-order variables have already been calculated using the standard linearized Boltzmann 
code such as CMBFAST. 
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10.2 Second-order bispectrum from products of the first-order terms 



10.2.1 A worked example 

We shall now focus only on the products of the first-order perturbations and we choose to analyze just 
some of these contributions to offer the reader an example of the analysis one can perform. We warn 
the reader that the full analysis is under progress and that the full results will be presented in Ref. [85] . 
including the intrinsically second-order perturbations which are equally important, and the contribution 
from perturbing the recombination history [721 1731 ^74]- 

For the products of the first-order perturbations, from now on we will consider the following non-zero 
four cases for the source terms, S\ 1 \ x \ 3 , which have been analyzed in Ref. [84] (for notational simplicity 
we shall omit the superscripts (1)): lj 



Sooo = 4*T / « (A;i)A 1 (fc 2 ) + [2r'(5 e + $)(fci) + S^^Ao^), 

Sno = J-{-5 r 'v (k 1 )v (k 2 ) + 2k 1 (y + $)(k 1 ) £ (2£ + 1)A,(£; 2 )}, 

e=odd 

Sioi = 2iV3{r , v (k 1 )(45 e + ^ + 2A -A 2 )(k 2 ) 
+Ak 1 <^(k 1 )(A - *)(fea) + fciA (fci)(* + 

Sin = 2J^{7T , v (k 1 )v (k 2 )-k 1 (y + <Z>)(k 1 ) J2 (2€ + l)A,(fc2)}. 

£=odd 



(450) 



In particular it is recognizable the contribution to 5ooo from the example just discussed. From these 
results we find that S\ 1 \ 2 \ 3 does not depend on i.e., S\ 1 x 2 \ s = S\ 1 \ 2 \ 3 (ki,k 2 ,r). Note also that 
S n(ki,k 2 ,r) = S 10 i(k 2 ,ki,r). 

10.2.2 Bispectrum from products of the first-order terms 

For the four non- vanishing combinations of Ai , X 2 , and A3 in Eq. (|450p , we rewrite the expression for the 
bispectrum, Eq. ([441 p . as 

R \ " R (Ai,A 2 ,A 3 ) u ci 

Ai A2A3 

o(0,0,0) „(1,1,0) „ B (1,0,1) , R (l,l,2) , B Cl M k 1 x 



where we have used B^f^ = B^f^ , and defined 



B £iM 3 = -MhbtaChCit + cyclic, (452) 



14 The following terms are in fact the only ones that have been fully analyzed numerically; they correspond to 
the source term S^JQs., 77) as given in Ref. [84], Eq. (3.2). However notice that that expression does not include all 
the products of first-order perturbations; the complete expression for the source term will be given in Ref. [85] . 
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and 




ik x ;: x ^ = -T.J 47r . ..... . 



(2Ai + 1)(2A 2 + 1)(2A 3 + 1) 



1)1 pi Ol 



^ Ai A2 A3 ) 
x J dre~ T f[ I ' dk i k 2 i P Ci (ki)j Vi {rki)g h (ki)Sx 1 \ 2 \ 3 {ki,k 2 ,r) 

i=l 

+ perm. (453) 

To proceed further, it turns out to be useful to simplify the expression by introducing the following 
notation for the integral over k that appears many times: 

[x]$(r) = ~J dkk 2+n P ( (k)j e ,(rk)g e (k)x(k,r). (454) 

This function corresponds to the existing functions in the literature in the appropriate limits. For 
example, for x(k, r) = tt/2, this function is the same as /3«!J (r) introduced in [86]. In fact, we find that an 
order-of-magnitude estimate of [a:]^/(r) is given by [a;]^?(r) ~ 2f3^)(r)/ir x x(k = £'/r,r) for a smooth 
function of x(k,r). As p^, (r) is a sharply peaked function at the decoupling epoch, r = r*, we find 
that [x]^)(r) is also sharply peaked at r = r*. With these tools in hand, one can calculate B^f^ , 

^h'lth > ^h%t£a ' an< ^ ^1^3 ^ e ^° no ^ ^°^ ner here the reader with the details of this computation 
whose details can be found in Ref. [S3]. Rather here we prefer to go straight to the results of the analysis 
of these bispectra which can be particularly instructive as far as their shape and their contamination to 
primordial non-Gaussianity are concerned. 



10.3 Shape and signal-to-noise of the second-order bispectrum from 
products of the first-order terms 

One of the motivations for calculating the second-order bispectrum is to see how much the second-order 
effects in gravity and the photon-baryon fluid contaminate the extraction of the primordial bispectrum. 
If, for example, the predicted shape of the second-order bispectrum is sufficiently different from that of the 
primordial bispectrum, then one would hope that the contamination would be minimal. To investigate 
this, we shall compare the numerical results of the second-order bispectrum with the so-called "local" 
model of the primordial bispectrum. 

We extract the first-order perturbations from the CMBFAST code. We use the following cosmological 
parameters: f?A — 0.72, Q m = 0.23, Qf, = 0.046, h = 0.70, and assume a power law spectrum, oc k n ~ 4 , 
with n = 1. We determine the decoupling time, r]*, from the peak of the visibility function. In this model 
we have crjo = 14.9 Gpc and cr/* = 288 Mpc. While the most of the signal is generated in the region 
of the decoupling epoch, in the low-£ regime we must also take into account the late time contribution 
due to the late integrated Sachs- Wolfe effect; thus, we integrate over the line-of-sight, r, in the following 
regions: c(t]q — 57/*) < r < c(rjo — 0.77/*) for i > 100, and < r < 0(770 — 0.77/*) for I < 100. The step size 
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is Ar = O.l??* around the decoupling epoch, and we use the same time steps used by CMBFAST after the 
decoupling epoch. 

The local primordial bispectrum is given by |35] 

POO 

B he2i3 = 2I hh£3 / r 2 drb L h (r)bf 2 (r)bg L (r) + cyclic , 
J 

(455) 

where 

bf(r) = - / k 2 dkPz(k)g$!(k)j e (kr), 
vr Jo 

bf L (r) = ~ / k 2 dkh L g$!(k)j e (kr). (456) 

Note that our linear transfer function, ge(k), is related to that of [35], gj^f(k), by gi{k) = §<?yf (k). 

Figure [T] shows the shape of the bispectrum generated by the products of the first-order terms f 
selected in Eq. (|450p . and compares it to the primordial local bispectrum, for £3 = 200. Both shapes 
(second-order and primordial) have the largest signals in the squeezed triangles, t\ <C £2 ~ £3- This is 
an expected result: both the local primordial bispectrum and the second-order bispectrum that we have 
computed here arises from the products of the first-order terms, also products in position space. However, 
these two shapes are slightly different when £1/^3 is not so small (^1/^3 = 0(0.1)): the ways in which the 
radiation transfer function (which gives the acoustic oscillations) enters into the bispectrum are different 
for the products of the first-order terms and the primordial bispectrum. The primordial bispectrum 
contains ji{kr*)gi{k), whereas the second-order bispectrum contains ji(kr*)gi(k)x(k,r*) where x = Aq, 
vq, etc., also has the oscillations. Therefore, the second-order bispectrum has more interferences between 
multiple radiation transfer functions. Moreover, the second-order effects contain derivatives that the local 
primordial effects do not have, which also makes the details of the two shapes different. 

Notice, in particular, that most of these gradients in the source term, Eq. (|134p . are contracted with 
the direction vector, n. There is only one term that has a scalar product of two wave- vectors, ki - ks, 
which vanishes in the squeezed limit. The resulting bispectrum, Eq. (|45ip . resembles that of a local 
form, except for the extra powers of k coming from the derivatives. These extra powers of k will affect 
the scale-dependence of the bispectrum, i.e., the second-order bispectrum is no longer scale- invariant. 
Nevertheless, the largest signal of the bispectrum still comes from the squeezed configurations, as the 
number of extra powers of k from the derivatives in the source term is not large enough to change the 
fact that we have the largest contribution when one of ki, &2, and k^ is very small. In other words, 
schematically the bispectrum looks like B(ki,k2,kz) ~ {h" 11 k™ 1 ) / {k^k^) + cyclic, where m\ and m2 
are the extra powers of k from the derivatives. Therefore, the largest contribution is in the squeezed 
configurations as long as rrii < 3. Figure [2] shows the same for £3 = 1000. The results are similar to those 
for £3 = 200, but the acoustic oscillations are more clearly visible. 

We can quantify the degree to which the second-order and the primordial bispectra are correlated, as 
well as the expected signal-to-noise ratio of the second-order bispectrum, following the general definitions 
summarized in Sec. El Notice that we are ignoring the noise contribution. In other words, we shall only 
consider ideal cosmic-variance limited experiments with full sky coverage, which however, at least for the 
multipole maximum multipole of l ma x ~ 2000 we will consider, is a good reference for an experiment like 
Planck (see, e.g., [35]). 
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Figure 1: Shape dependence of the second-order bispectrum from products of the first- 
order terms (top) and that of the local primordial bispectrum (bottom). We show 

W4LALA,){QTXh m3 )' X l^f x 1022 as a Unction of h/£ 3 and £ 2 /£ 3 where £ 3 = 200. 
Both shapes have the largest signals in the squeezed triangles, l\ <C £ 2 ~ £3. 




Figure 2: Same as Fig. [T] for £3 = 1000. The acoustic oscillations are clearly seen. 
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Figure 3: Signal-to-noise ratios for the local primordial bispectrum for /nl = 1 (dashed), and the 
second-order bispectrum from the products of the first-order terms (solid), for an ideal full-sky 
and cosmic- variance-limited (noiseless) experiment. 



The signal-to- noise ratio is given by Eq. (1354p . In Fig. [3] we show the cumulative signal-to- noise 
ratio, summed up to ci maximum multipolc of l m ax) 

of the primordial bispectrum, assuming /nl = 1 and 
) 1 ^ 2 , as well as that of the second- 



{F 2 nd,2nd) 1/2 ■ In both cases S/N 



ignoring the second-order bispectrum, i.e., (S/N) pr i m = (Fprim^rim) 1 
order bispectrum, ignoring the primordial bispectrum, i.e., (S/N) 2n d 
increases roughly as S/N of imax (or oc yj Npi x where Npi x is the number of independent pixels in the 
map). A larger contribution to the second-order bispectrum at £ < 50 comes from the terms involving the 
Integrated Sachs- Wolfe effect. The signal-to- noise ratio of the second-order bispectrum reaches ~ 0.4 at 
4nax = 2000; thus, this signal is undetectable. While our calculation includes the temperature anisotropy 
only, including polarization would increase the signal-to-noise by a factor of two at most, which would 
not be enough to push the signal-to-noise above unity. 

While the total signal-to-noise does not exceed unity, it may still be instructive to show which terms 



of B 



and are more important than the others. To do this we show the following quantity: 



S_ 



E 



2<e 1 <£ 2 <i 3 



(J 



1/2 



(457) 



where a, b = 1, 2, 3, 4, and correspond to (0,0,0), (1, 1,0), (1,0, 1), (1, 1,2), and Ci, respectively. 



The results are shown in Fig. [H We find that (S/N)2 n d is dominated by 



whereas it is dominated by Bf^^ for I > 100 (see the top panel). Among B^^' 
inant term is (1,0,1) (the bispectrum from the second-order dipole created by the first-order dipole 
and monopole). The second most dominant is (0,0,0) (from the second-order monopole created by the 



(Ai,A2,A3) 



for £ < 100, 



(Ai,A2,A3) 



the most dom- 
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Figure 4: Absolute values of the contributions to the signal-to-noise ratio from each component, 
(S/N) ab , as defined by Eq. ( 14571) . 
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Figure 5: The cross-correlation coefficient between the second-order bispectrum from the products 
of the first-order terms and the local primordial bispectrum. 

first-order monopole) for I < 400 and (1,1,0) (from the second-order monopole created by the first-order 
dipole) for I > 400. The cross terms (middle and bottom panels) are sub-dominant compared to the auto 
terms (top panel) at all multipoles. 

How similar are the second-order and the primordial bispectra? In Fig.[5]we show the cross-correlation 
coefficient, between the local bispectrum and second-order bispectrum from the products of the first-order 
terms given in Eq. (|450p . The cross-correlation coefficient (see Eq. (I355p . reaches ~ 0.5 for l max = 200, 
and the shapes for £3 = 200 are shown in Fig. [TJ After £ max = 200 the correlation weakens, and reaches 
~ 0.35 at ^ max = 1000, and the shapes for £3 = 1000 are shown in Fig. [2 These results show that 
the second-order bispectrum from the products of the first-order perturbations and the local primordial 
bispectrum are fairly similar, with a sizable correlation coefficient. 

How large is the contamination of the primordial bispectrum? The level of contamination is measured 
by the effective non-linearity parameter given by Eq. (135711 . 

In Fig. [6] we show as a function of the maximum multipoles, £ ma x- We find that /j^i? reaches 
the maximum value, ~ 0.9, when the correlation coefficient reaches the maximum at £ max ~ 200, but 
then decreases to ~ 0.5 at £ max ~ 2000. Therefore, we conclude that the contamination of the primor- 
dial bispectrum due to the second-order bispectrum from the terms in Eq. (|450p is negligible for CMB 
experiments. 

Finally, one can also calculate the \-a uncertainty of /nl 5 <5/nl 5 with the second-order bispectrum 
marginalized over. This is given by <5/nl = J (-^)p r im,priml/NL=i- Fig. [7] shows that an increase in the 
uncertainty of /nl due to marginalization is totally negligible. 
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Figure 6: Contamination of the local primordial bispectrum as measured by (Eq (I357P ). 



11 Conclusions 

In these review we have addressed a basic question in cosmology: how a primordial NG propagates into 
an observable like the CMB anisotropy. Answering this question is fundamental as it will help us in 
getting some knowledge about the way the primordial cosmological perturbation was generated at the 
very early stages of the evolution of the universe. In the first sections we have shown how to set the initial 
conditions at second-order for the (gauge-invariant) CMB anisotropy when some source of primordial NG 
is present. This was more or less straightforward because on large angular scales is basically gravity which 
dictates the non-linear dynamics. On small angular scales the computation of the second-order effects in 
the CMB anisotropy is far more difficult, as so many sources of non-linearities appear. In this review, 
we have focussed ourselves on the study of the second-order effects appearing at the recombination era 
when the CMB anisotropy is left imprinted. We have shown how to derive the equations which allow 
to evaluate CMB anisotropics, by computing the Boltzmann equations describing the evolution of the 
baryon-photon fluid up to second order. This allows to follow the time evolution of CMB anisotropics (up 
to second order) on all scales, from the early epoch, when the cosmological perturbations were generated, 
to the present time, through the recombination era. Through some analytical and simplified example 
we have also shown how estimate the contamination of the recombination secondary effects onto the 
detection of primordial NG. More refined numerical results confirm our estimates and are also reported 
here. It goes without saying that this line of research should be pursued until the level of accuracy in the 
theoretical prediction is reached and is comparable to the one provided by the current or future satellite 
experiments. 
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Figure 7: Projected uncertainty of /nl with (dashed) and without (solid) the second-order bis- 
pectrum marginalized over. 
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A Energy-momentum tensor 

A.l Energy- momentum tensor for photons 

The energy-momentum tensor for photons is defined as 



where g is the determinant of the metric ([T]) and / is the distribution function. We thus obtain 




(A.l) 




(A.2) 



(A.3) 



(A.4) 



106 



where p 7 is the background energy density of photons and 
is the quadrupole moment of the photons. 

A. 2 Energy-momentum tensor for massive particles 



The energy-momentum tensor for massive particles of mass m, number density n and degrees of freedom 
9d 

9d_ f &Q_ Q^Q V 
-g) (2vr) 3 Q° 



T mu~ r— \ s 3 n0 9m , (A.6) 



where g m is the distribution function. We obtain 



-^m ~~ Pm — Pm 



l + Stt + btt) , (A.7) 



2 



± m — 



p m vl = -e* + *p m + + 5(1)^(1)^ (A .8) 

T !n j = Pm (f s ^ + v m v m = p m ^ + , (A.9) 

where p m is the background energy density of massive particles and we have included the equilibrium 
temperature T m . 

B Solutions of Einstein's equations in various eras 
B.l Matter-dominated era 

During the phase in which the CDM is dominating the energy density of the Universe, a ~ rj 2 and we may 
use Eq. (|35|) to obtain an equation for the gravitational potential at first order in perturbation theory 
(for which $W = *W) 



+m<5> {1) ' = 0, (B.l) 

inc 

CDM velocity can be read off from Eq. (|31|! 



which has two solutions = constant and = T-L/a? . At the same order of perturbation theory, the 



= - A^ 1 ) . (B.2) 
3 rl 



The matter density contrast 5^ satisfies the first-order continuity equation 



Going to Fourier space, this implies that 
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4 I, =4 1, (0) + ^ 1 », (B.4) 

where 8fp (0) is the initial condition in the matter-dominated period. 

At second-order, using Eqs. (|35p and (|33|) and the fact that the first-order gravitational potential is 
constant, we find and equation for the gravitational potential at second-order 

*(2)" + 3 ^(2)' = Sm} (B .5) 
Sm = -d k ^d k ^ + N = -d k ^d k ^ + y ^ , 
whose solution is 



(B.6) 



with W{rj) = Wo/a s (do = 1) the Wronskian obtained from the corresponding homogeneous equation. 
In Eq. (|B.6|) (0) represents the initial condition (taken conventionally at rj — >■ 0) deep in the matter- 
dominated phase. 

From Eq. (|36p . we may compute the vector perturbation in the metric 



^(d^V^-^W 1 )*) , (B.7) 



V 2 

where we have made use of the fact that the vector part of the CDM velocity satisfies the relation 



(3 



B.2 Radiation-dominated era 

We consider a universe dominated by photons and massless neutrinos. The energy-momentum tensor 
for massless neutrinos has the same form as that for photons. During the phase in which radiation is 
dominating the energy density of the Universe, a ~ rj and we may combine Eqs. (|30p and (|35p to obtain 
an equation for the gravitational potential at first order in perturbation theory 



*(D" + 4H*«'--V 2 *W = uqQ-Y + -V 2 Q (1) 
3 3 



9 7/2^' n (l)< 
2 rt v 2 J ' 



(Bi 
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where the total anisotropic stress tensor is 

n Ti = f^i + f 111 ^- (B-9) 

We may safely neglect the quadrupole and solve Eq. ()B.8P setting u± = ^±7]. Then Eq. (jESJ), in 
Fourier space, becomes 

u" + -u'+(^--^)u = 0. (B.10) 

This equation has as independent solutions = ji(kr)/y/3), the spherical Bessel function of order 1, 
and U- = ni(kr]/^/3), the spherical Neumann function of order 1. The latter blows up as r\ gets small 
and we discard it on the basis of initial conditions. The final solution is therefore 

T m T /n, , sin (kr]/V3) - (/cr ? /V3)cos(A;r ? /v / 3) 

where ^(^(O) represents the initial condition deep in the radiation era. 

At the same order in perturbation theory, the radiation velocity can be read off from Eq. (|3ip 



«P = -^l^— 2 - ■ (B.12) 
At second order, combining Eqs. (f30|) . ([35]) . we find 



^r(2)" + 4 ^^(2)' _ I V 2 ^ (2) = 5 7 , (B.13) 
3 



^ ' 3 

- -(d fc $ (1) d fe $ (1) + d fc ^ (1) d fc ^ (1) - d fc $ (1) d fc ^ (1) ) 
3 

+ W q(2)' + -V 2 Q( 2 ) + + ^( 1 ))V 2 Q( 1 ) , 

3 3 



3 a * 



v 2 

+ ^ 2 ^^-9H 2 ^(^)n^,), (B.14) 



whose solution is 
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(2) / " j j ®? (v)$ { - ] (V')-$ W (V)*? (rf) 



* {2) = nl. + / d V 

Jo 



W(rf) 



(B.15) 



where W{rj) = (a(0)/a) 4 is the Wronskian, and ^\.J m is the solution of the homogeneous equation. 
The equation of motion for the vector metric perturbations reads 



- - V V + 4% V 



2$(i)>'$(i)+H 2 ^^ 



PT 



_ 2-H 2 ( + ^v^ + 2^A { £\^ + 2^A$%« j 

VPt 7 />t r PT 

+ 2 ($( 1 ) - *(i))^Lt,(i)i + 2($« - xpW)^^ 



PT 



PT 



(B.16) 



where />r is the total background energy density. The Einstein equations for a universe filled by CDM 
and a relativisitc component can be foun in Ref. 



C Linear solution of the Boltzmann equations 

In this section we will solve the Boltzmann equations at first order in perturbation theory. The interested 
reader will find the extension of these formulae to second order in Ref. [39J . The first two moments of the 
photon Boltzmann equation are obtained by integrating Eq. (|129p over d£l n /4:ir and dQ a n l /4-7T respectively 
and they lead to the density and velocity continuity equations 



(C.l) 
(C.2) 



v^> + + + = -r' («W 

where IP- 7 is the photon quadrupole moment, defined in Eq. (|192p . 

Let us recall here that fffl = = J d£lA^ /Air and that the photon velocity is given by Eq. (|183p . 

The two equations above are complemented by the momentum continuity equation for baryons, which 
can be conveniently written as 



' T 



(C.3) 



where we have introduced the baryon-photon ratio R = 3pb/(Ap y ). 

Eq. (lC.3j) is in a form ready for a consistent expansion in the small quantity r" 1 which can be 



performed in the tight-coupling limit. By first taking v^ 1 
relation in the left-hand side of Eq. (|C3|) one obtains 



.(!)< 



at zero order and then using this 



R 



(C.4) 
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Such an expression for the difference of velocities can be used in Eq. (|C.2|) to give the evolution equation 
for the photon velocity in the limit of tight coupling 

+ n-^—v^ + - A ^- + = . (C.5) 
7 1 + R 7 41 + R 

Notice that in Eq. ([05]) we are neglecting the quadrupole of the photon distribution TlWv (and all the 
higher moments) since it is well known that at linear order such moment (s) are suppressed in the tight- 
coupling limit by (successive powers of) 1/r with respect to the first two moments, the photon energy 
density and velocity. Eqs. (IC.lj) and (lC.5j) are the master equations which govern the photon-baryon fluid 
acoustic oscillations before the epoch of recombination when photons and baryons are tightly coupled by 
Compton scattering. 

In fact one can combine these two equations to get a single second-order differential equation for the 
photon energy density perturbations Aqq' . Deriving Eq. (IC.ip with respect to conformal time and using 
Eq. ()C.5P to replace diV^ 1 yields 

R 



(A«''-4*«'')+^ T f^(A«'-4*«' 
-c 2 s V 2 (A&> - 4*W) = V + J^M , (c.6) 



where c s = l/y / 3(l + R) is the speed of sound of the photon-baryon fluid. Indeed, in order to solve 
Eq. (|C,6ft one needs to know the evolution of the gravitational potentials. We will come back later to the 
discussion of the solution of Eq. (|C,6I) . 

A useful relation is obtained by considering the continuity equation for the baryon density perturba- 
tion. By perturbing at first order Eq. (|173p we obtain 

S b V + v \i ~ 3 ^ (1) ' = • ( C - 7 ) 
Subtracting Eq. ([GTj) form Eq. ([C3]) brings 

^'-ffi' + t^-v (1)i ),* = 0, (C.8) 
which implies that at lowest order in the tight-coupling approximation 

Affi = , (C9) 

for adiabatic perturbations. 



C.l Linear solutions in the limit of tight coupling 

In this section we briefly recall how to obtain at linear order the solutions of the Boltzmann equa- 
tions (|C.6|) . These correspond to the acoustic oscillations of the photon-baryon fluid for modes which are 
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within the horizon at the time of recombination. It is well known that, in the var iable (A$ - 4tfW) 
the solution can be written as (HTJ [88] 

[1 + i?(r/)] 1 / 4 (A£ ) - 4fW) = Acos[kr s {r])] + B sin[fcr s (r/)] 




where the sound horizon is given by r s {r]) = Jq drf c s (rf) , with R = 3p&/(4p 7 ). The constants A and B 
in Eq. (IC.lOj) are fixed by the choice of initial conditions. 

In order to give an analytical we will use some simplifications following Ref. [891 I42j . First, for 
simplicity, we are going to neglect the ratio R wherever it appears, except in the arguments of the varying 
cosines and sines, where we will treat R = R* as a constant evaluated at the time of recombination. In this 
way we keep track of a damping of the photon velocity amplitude with respect to the case R = which 
prevents the acoustic peaks in the power-spectrum to disappear. Treating R as a constant is justified by 
the fact that for modes within the horizon the time scale of the ocillations is much shorter than the time 
scale on which R varies. If R is a constant the sound speed is just a constant c s = 1/^/3(1 + R*), and 
the sound horizon is simply r s (rf) = c s T]. 

Second, we are going to solve for the evolutions of the perturbations in two well distinguished limiting 
regimes. One regime is for those perturbations which enter the Hubble radius when matter is the dominant 
component, that is at times much bigger than the equality epoch, with k <C k eq ~ rj~ q , where k eq is the 
wavenumber of the Hubble radius at the equality epoch. The other regime is for those perturbations with 
much smaller wavelenghts which enter the Hubble radius when the universe is still radiation dominated, 
that is perturbations with wavenumbers k 3> k eq ~ f]eq- ^ n f ac t we are interested in perturbation modes 
which are within the horizon by the time of recombination rj*. Therefore we will further suppose that 
r/* S> Tjeq in order to study such modes in the first regime. Even though 77* 3> ij eq is not the real case, it 
allows to obtain some analytical expressions. 

Before solving for these two regimes let us fix our initial conditions, which are taken on large scales 
deep in the radiation dominated era (for 77 — >- 0). During this epoch, for adiabatic perturbations, the 
gravitational potentials remain constant on large scales (we are neglecting anisotropic stresses so that 
$(!) ~ xfC 1 )) and from the (0 — 0)-component of Einstein equations 




On the other hand, from the energy continuity equation ([C.lj) on large scales 

A$ - = const. ; (C.12) 



from Eq. (|C.llh the constant on the right-hand side of Eq. (|C,12j) is fixed to be — 6$^^(0); thus we find 
B = and A = -6*W(0). 




where ujq = kc s . 
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C.2 Perturbation modes with k <C k eq 

This regime corresponds to perturbation modes which enter the Hubble radius when the universe is 
matter dominated at times rj 3> rj eq . During matter domination the gravitational potential remains 
constant (both on super- horizon and sub- horizon scales), as one can see for example from Eq. (|B.ip . and 
its value is fixed to ^^(k, rj) = ^^^(O), where $^^(0) corresponds to the gravitational potential on 
large scales during the radiation dominated epoch. Since we are interested in the photon anisotropies 
around the time of recombination, when matter is dominating, we can perform the integral appearing in 
Eq. (IC.lOp by taking the gravitational potential equal to its value during matter domination so that it is 
easily computed 

2 / driWtf) sinMr? - rj')] = K -L (1 - cos(^)) . (C.14) 

Jo 1U COQ 



Thus Eq. ()C.13P gives 

A$ - 4^« = £$(i)( ) C os(u r]) - ^$W(0) . (C.15) 

The baryon-photon fluid velocity can then be obtained as diV^ = -3(A^ ) - 4$rW)'/4 from Eq. (^Hj) . 
In Fourier space 

9 



ihvM = -$«(0)sin(a;o7 ? )a;o, (C.16) 
where, going to Fourier space, e^t^ — > i/^u^ i (k) and 



v?" = -ij y^ (1) (0) sin( Wo r?)c s , (C.17) 



,(!)* 



since the linear velocity is irrotational. 

Notice that, under the approximations that R = const, and = VfW = const., it is easy to find a 
more accurate solution to Eq. ()C.6| ) which better accounts for the presence of the baryons (giving rise to 
the so called baryon-drag effect, as clearly explained in Refs. [8T 1 |88 | IHT]) 



A 



(i) 
oo 



A 



(i) 
no 



■(0) + (1 + R)$W] cos(kr s ) - (1 + . (C.18) 



C.3 Perturbation modes with k ^> k eq 

This regime corresponds to perturbation modes which enter the Hubble radius when the universe is still 
radiation dominated at times r\ <C rj eq . In this case an approximate analytical solution for the evolution 
of the perturbations can be obtained by considering the gravitational potential for a pure radiation 
dominated epoch, given by Eq. (|B. 1 1|) . For the integral in Eq. (|C.13p we thus find 

V ^(^smicooirj-rj')} = - J- COS (u, or? ) , (C.19) 

where we have kept only the dominant contribution oscillating in time, while neglecting terms which 
decay in time. The solution (|C.13j) becomes 

A ( q - 4tfW = 6$ (1) (0) cos(w ?7) , (C.20) 
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and the velocity is given by 

= -i^|* (1) (0)sin(a;o»7)c a , (C.21) 

Notice that the solutions ()C.20|) - (|C.21|) are actually correct only when radiation dominates. Indeed, 
between the epoch of equality and recombination, matter starts to dominate. Full account of such a 
period is given e.g. in Section 7.3 of Ref. [12], while its consequences for the CMB anisotropy evolution 
can be found e.g. in Ref. |90j . We refer to [39], where an alternative way to solve for the acoustic 
oscillations in this regime is displayed which turns out to be useful for the corresponding computation at 
second order in Sec. 17.51 



Symbol 


Definition 


Equation 


Mr . 1 


Vjrl cLV 1 LcLLlUllcLl JJU LdlLldlo 111 1 UlooOll g,cLU.g,ti 


(Til 




2nd-order vector perturbation in Poisson gauge 


(W 


AlJ 


2nd-order tensor nertnr nation in Poisson raiiffe 


(HI 


n 
i 


Conformal time 


CD) 


f 


Photon distribution function 


(1571) 


g 


Distribution function for massive particles 


flETJ & (jl53|) 


f(i) 


i-ih. order perturbation of the photon distribution function 


(1821) 


Ai) 
Jim 


A/Tnmp'n1"c; of flip Tinofori HicitT'i niitinvi fn n pfi nn 


(11 001) 


C(r>) 


( jolliGion fpTin 

V. v.v 1 1 1 1 v.v 1 1 LjC-1111 




l> 


A/T fi crn l tii H p pit nn pit r\v\ tyi nui fiTTTn m 1 ti — titj * 1 
IVld^lll L LIU." Ul JJ11CLU11 lllUlildlL Hill I — I 


(IK?1) 


n' 


ProDH rat ion direction 


(1651) 


A^(x\ n\ rj) 


First-order fractional energy photon fluctuations 


(|128|) 


A^(x\ n\ vj) 


Second-order fractional energy photon fluctuations 


(|132|) 


n e 


Electron number density 


(|163|) 


S e (6 b ) 


Electron (baryon) density perturbation 


m 


k 


Wavenumber 


fll39j) 


v m 


Baryon velocity perturbation 


(I142p & (I143p 


w (2)< 
U CDM 


Cold dark matter velocity 


(|198|) 


,.(2)< 

7 


Second-order photon velocity 


(|20l|) 


Sim 


Temperature source term 


(jl4()j) 


T 


Optical depth 


(fTBTD 


P~f(pb) 


Background photon (baryon) energy density 


(|195|) 




CMB angular-averaged bispectrum 


(135 1|) 


F 


Fisher matrix for the amplitudes of the bispectra 


(|352|) 


f con 
iNL 


contamination to primordial non-Gaussianity 


(I357|) 
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